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Abstract 



Naturally small neutrino masses can arise in some grand unified models. The 
mechanism of neutrino mass generation in these models typically requires the exis- 
tence of neutral heavy leptons. We study the low-energy phenomenology of these 
new fermions. Concentrating on loop corrections due to neutral heavy leptons, we 
examine how the flavour- conserving leptonic decays of the Z boson, universality 
breaking in these decays, and the W boson mass depend on the mass and mixings of 
the neutral heavy leptons. Working within the framework of a superstring-inspired 
SU(2) L x U(l)y model, we show that these flavour-conserving processes have some 
virtues over the traditionally considered flavour-violating decays. 
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Chapter 1 
Introduction 

1.1 The problem of small neutrino masses 

Why are macroscopic things around us massive ? We know they, as composite ob- 
jects, get their mass from the mass of their constituents and from the interactions 
among the constituents. But if the constituents are elementary particles, where 
is their mass coming from ? For theoretical physicists mass represents a rather 
unwelcome pollution of their elegant massless theories. They believe the world is 
essentially symmetric and it is the mass that makes this symmetry hard to see (sym- 
metry is hidden or spontaneously broken). Therefore they build massless theories 
possessing beautiful symmetries and then break these symmetries to give particles 
their masses. The exact mechanism of symmetry breaking is not known and may 
not be known for a long time. It is one of the most attractive problems of particle 
physics. Some argue we have to wait for the ultimate solution until physics at the 
Planck scale (10 19 GeV) is understood, others suggest optimistically that 10 5 GeV 
might do. The ultimate theory should predict the masses of all elementary particles 
in agreement with experiment. 

It is believed the origin of the weak boson masses (associated with the electroweak 
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2 



symmetry breaking) is better understood than that of the fermion masses (flavour 
symmetry breaking). Concentrating on the fermions, we list the quark and lepton 



content of the so-called standard model of electroweak interactions in Table 1.1 



along with their masses (T|]. These masses are not predicted by the standard model. 



1. 


Mass [MeV] 




2. 


Mass [MeV] 




3. 


Mass [GeV] 


u 


4.2 




c 


1 100 




t 


176 


d 


7.5 




s 


150 




b 


4.2 









v n 







v T 





e 


0.5110 






105.6 




T 


1.784 



Table 1.1: The three families of quarks (u,c,t,d, s,b) and leptons (neutrinos 
v e , z/^, v T \ charged leptons e, /i, r) and their masses according to the standard model 
of electroweak interactions 



The quark and charged lepton masses represent an experimental input, 9 parameters 
out of the total 17 present in the standard model [J The neutrinos are postulated 
as massless. This postulate, however, is based on the assumption that the neutrino 
is the only fermion without a right-handed field - an asymmetry going against the 
spirit of symmetric theories which work with both left-handed and right-handed 
fields. Therefore the massless neutrinos are not natural in the standard model and 
none of the twelve fermion masses is actually predicted. 

While the ultimate solution for mass prediction may be far from us, it is worth- 
while to think of partial steps which could bring us closer to it. A popular strategy 
aims at the reduction of the twelve independent mass parameters. From grand uni- 
fied theories (GUT's) which describe strong and electroweak interactions as different 
manifestations of a single force, there are hints the fermion masses are related to each 
other by simple formulae at very high energies (GUT scale ~ 10 16 GeV) and it is in 



1 Of the remaining 8 parameters, the four CKM mixing matrix parameters are also closely related 
to the origin of the quark masses. 

2 For example in the simplest version of 50(10) GUT, all masses in a given family are equal at 
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the transition to our low-energy world [] where masses pick up different corrections 
and end up in the array of seemingly unrelated numbers shown in Table O . 

In the past the efforts to partially explain fermion masses in grand unified models 
experienced great difficulties facing the problem of small neutrino masses. In the 
discussion above, we dismissed the way the standard model postulates the zero mass 
neutrinos. What is then the experimental basis for the claim the neutrino masses are 
small, possibly zero ? Ironically, even though neutrinos are possibly the second most 
abundant elementary species in the Universe, we do not know exactly what their 
mass is. This is partly due to the smallness of this mass and partly because it is so 
hard to detect neutrinos (a typical cross section is 1CT 43 cm 2 for a 10 MeV neutrino, 
more than twenty orders of magnitude below the cross section for an electron of 
the same energy). Experiments set the following upper limits on the masses of the 
electron, muon and tau type neutrino f2j: 





< 


7.2 eV 




< 


270 keV 




< 


31 MeV 



The neutrino masses are, even at their allowed maxima, strikingly small when com- 
pared with the masses of the charged leptons and quarks within their families. For 
example, within the first family of the standard model, u and d quarks have masses 
of a few MeV, and the electron has a mass of 0.511 MeV (see Table |1 . 1|) . The simple 
formulae relating the fermion masses we mentioned earlier, put different fermion 
masses on the same scale: it is thus conceivable in GUT's that masses of the first 
family are scattered around a few MeV, but the mass of the electron type neutrino 
is a deep mystery, being at least five orders of magnitude below this natural scale 

the GUT scale. 

3 By low energies we mean here energies up to a few hundred GeV and by very low energies 
those up to a few GeV. 
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for the first family. The difference of five orders cannot be explained by corrections 
masses pick up running from the GUT scale to low energies. Something else must 
be involved. Similar behaviour, although perhaps less pronounced, can be observed 
among the second and third family members. 

1.2 Solutions to the problem 

Before we start to discuss possible solutions to our problem, let us be more specific 
about the relation between GUT's and low-energy theories, such as the standard 
model. While the standard model [[| |], [5|, || has a natural scale of 10 2 GeV, GUT 
models describe physics actively operating at the GUT scale of 10 16 GeV; at 
the same time they should explain low-energy data at least as well as the standard 
model. GUT models, although more elegant than the standard model, are also more 
complicated: they have a richer gauge structure and symmetry breaking sector and, 
especially in the case of Eq, a fermion sector with more particles. A thorough 
discussion of neutrino masses in these models is beyond the scope of this work. 
Here we are mainly interested in the phenomenology of neutral heavy leptons at low 
energies. Fortunately, at low energies most of the extraneous baggage associated with 
the complicated structure of GUT models has very little impact - it is integrated out 
and the remaining effective theory often represents just a minimal extension of the 
standard model of electroweak interactions. This is in fact no big surprise. It simply 
reflects a very good agreement of experimental data with the standard model and 
the fact the GUT's extend the standard model rather than replace it. The physics 
active at high energies thus decouples at the scales currently accessible to us. For 
example, in the model studied in this thesis, it is just two extra neutrino fields per 
family which enlarge the standard model. Therefore we will focus our discussion 
on these minimal extensions, referring to unification models as motivation for a 
particular simple extension of the standard model. Interestingly, basic directions 
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in the theoretical treatment of neutrino masses can be followed with just minor 
extensions of the fermion sector of the standard model, leaving the gauge structure 
SU(2)l x U(1)y and symmetry breaking sector intact. 

The simplest model with massive neutrinos one can immediately think of is 
a straightforward extension of the standard model. One can introduce the right- 
handed neutrino field missing in the standard model and treat the neutrino in the 
same way as all other fermions - as a massive Dirac particle. This means the neutrino 
mass is still not predicted and the simplest model thus fails to address the problem 
of the small neutrino masses. 

A possible solution was found by Yanagida and Gell-Mann, Ramond and Slansky 
in the famous see-saw mechanism p, ||. A simple low-energy see-saw model has the 
same fermion content as the simplest model just described, however, mass terms 
violating the total lepton number L are allowed. This leads to the description of 
neutrinos as Majorana fermions rather than Dirac; neutral heavy leptons (NHL's) 
are introduced into this theory as a necessary ingredient. In Chapter 2 we describe 
the difference between Majorana and Dirac fermions more formally. Here it suffices 
to say that a Dirac neutrino is a particle like all other fermions with left-handed and 
right-handed particle and antiparticle states, while a Majorana neutrino is a particle 
which is its own antiparticle and therefore comes in just two states described by a 
left-handed and a right-handed field. The neutrino is the only particle of the fermion 
content of the standard model that can possibly be described as a Majorana particle 
because it is neutral. 

The see-saw mechanism comes with the following relation for the mass of a 
neutrino m v : 



The mass D is a typical family mass (say 1 — 2 MeV for the first family) and is 



m v = D 2 /M N . 
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the NHL mass. This relation tells us the neutrinos become very light with respect 
to D due to a very large mass M^- Their mass is not predicted (Mjy is unknown); 
nevertheless its smallness is understood since a very large mass scale required for 
the Mjv is naturally expected in GUT's. 

Whether neutrinos are Dirac or Majorana particles is an issue by itself. Experi- 
ments which try to address it rely on the fact that Majorana neutrinos break lepton 
number conservation, whereas Dirac neutrinos respect it. A clear answer would be 
provided by the observation of a neutrinoless double beta decay, a so far unobserved 
process mediated only by Majorana neutrinos [11 • 



Theories with massive neutrinos are popular since some motivation for nonzero 
neutrino masses comes also from outside particle physics. Massive neutrinos could 
explain the mystery of missing solar neutrinos [O through matter enhanced time 
dependent neutrino oscillations, the so-called MSW effect ]13|. In cosmology, massive 
neutrinos could explain at least part of the dark matter puzzle. The possibility of this 
increased after COBE data on the anisotropy of the cosmic microwave background 
radiation were analysed: there are hints some 10 — 30 % of the dark matter is hot 
14| and neutrinos with masses between 2-7 eV make a good candidate. Cosmology 



further constrains (subject to plausible assumptions) the mass of each of the three 
neutrinos to m v < 25 eV (see Ref. Sec. 15.3.1). 

These astrophysical and cosmological indications are quite appealing and some 
authors argue on their basis that an 5*0(10) GUT with three nearly degenerate 
Majorana neutrino masses of about 2 eV could be the correct unified theory |T(J. 

The see-saw mechanism is not the only one which addresses the issue of the 



4 Another possibility of proving the Majorana nature of the neutrino was suggested as a result 
of theoretical studies of the origin of the lepton symmetry breaking. This symmetry can be broken 
explicitly (it is not a symmetry of the Lagrangian at any energy or temperature); it can also be 
broken spontaneously. If the lepton number L is a global symmetry which is broken spontaneously, a 
massless, pseudoscalar Majoron arises jlC| ]. The discovery of this particle would prove the Majorana 
nature of the neutrino. 
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smallness of the neutrino mass. An alternative was worked out in the class of 
models wherein neutrinos remain massless while other neutral leptons can acquire 
a large mass. Here it is argued there can be some global symmetry present, such 
as lepton number L, which prevents neutrinos from becoming massive. One such 



model could arise as a low-energy limit of a superstring-inspired Eq GUT [17]. The 
superstring inspiration consists of yet another (left-handed) neutrino field added to 
the fermion content of the standard model. Therefore we have three neutrino fields 
per family: the standard left-handed one, the right-handed one required by GUT's 
and the one needed by superstrings. This field content and the lepton number 
symmetry [] give rise to three massless neutrinos and to three massive Dirac NHL's, 
as described in Chapter 3. In this thesis I study phenomenological implications of 
this model (henceforth called 'our' model), with emphasis put on signatures of NHL's 
in precision data from LEP collider at CERN (leptonic widths of the Z boson) and 
Tevatron collider at Fermilab (the mass of the W boson, Mw). 

Both theories with massive and massless neutrinos (our model) naturally exhibit 
some of the properties familiar from the quark sector of the standard model. For 
instance, neutrino mixing may arise with an analogue of the Cabibbo-Kobayashi- 
Maskawa (CKM) mixing matrix for the leptonic sector. Individual lepton family 
numbers are expected to be violated, as is lepton universality and CP symmetry. 
The total lepton number L has already been discussed. What is not predicted by 
our model, in contrast with see-saw models, are time dependent neutrino oscillations 
and neutrinoless double beta decay. The physics of light neutrinos is thus richer 
in see-saw models. On the other hand the physics of neutral heavy leptons looks 
more promising in our model. The NHL's contribution to low-energy observables is 



5 We note the conserved lepton number is also present in the standard model and in the simplest 
model with massive neutrinos. However, in the former the conserved L is a consequence of the 
missing right-handed neutrino field and in the latter, the two neutrino fields per family are not 
enough to keep neutrinos massless. 
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proportional to their mixings and these depend on the ratio In see-saw models 
this ratio is normally very small since the ratio -j^- = m v is very small. In our model 
with massless neutrinos there is no such restriction and, consequently, mixings can 
be relatively large. 

To show that the differences between theories with massive neutrinos and our 
model should not be taken too seriously, we note that there is a variant of our model 
where the lepton number symmetry is slightly broken and the neutrinos are given 



a small mass WA 181 and there are variants of see-saw models where the restriction 



from the ratio is avoided (assuming certain symmetries in the neutrino mass 



matrix) and possibly large mixings of NHL's arise as a result [19, 20 



1.3 The phenomenology of NHL's 

To investigate the phenomenology of NHL's more closely, let us be more specific 
about how they enter low-energy observables via their mixings. For simplicity we 
neglect interfamily mixings in this section. The key point, formally derived and 
discussed in Chapters 2 and 3, is that the neutrinos taking part in weak interactions, 
V\ (I = e, fi, t) are no longer states of definite mass but are combinations of light 
(massless in our model) neutrinos v and NHL's N: 

v x = K L J + K H N. (1.3) 

Here Kl and Kh are mixing parameters related through K\ + K H = 1. The chance 
of finding the signature of an NHL is proportional to the size of K^, which is equal 
to the ratio we discussed in the previous section, K H = This picture is model 
independent, whether NHL's come from our model, or a see-saw one. In the standard 
model we have Kh = 0, Kl = 1 and the weak eigenstates V\ become identical with 
the mass eigenstates v . 
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The eigenstates v\ interact weakly by coupling to the mediators of electroweak 
interactions, W and Z bosons. Through the mixing K Hl so do NHL's. For example, 
the simplest way to discover NHL's would be to directly produce them at the Z 
boson factories (LEP collider at CERN and SLC at SLAC). However, to be produced 
directly, their mass would have to be smaller than the Z boson mass - this is the 
energy available at the Z factories. Since there has been so far no evidence for 
NHL production, we conclude that either they are more massive than the Z, or 
their mixing is small enough to suppress their production to such a degree that they 
escape detection. 

Another method to probe mixings of NHL's is an indirect one, through the 
measurement of the mixing parameter Kl of the light neutrinos in, for example, 
pion and beta decays. NHL's are not produced in the decays, nevertheless their 
existence could be revealed if K L is far enough from 1 to reduce the decay rates 
beyond the experimental uncertainties. These so-called universality constraints give 
us the best limits on mixings. 

While the direct production is sensitive to NHL masses only up to the mass of 
the Z boson Mz, indirect methods are sensitive only to mixings, not masses. What 
if the NHL is heavier than the Z boson ? Can we obtain some information on its 
mass ? In cases when particle physicists face a problem of probing masses of hy- 
pothetical particles which are larger than the energy currently available, they study 
contributions of these particles in radiative corrections (loops) to some observables. 
Since loops are higher order terms of perturbation theory, they represent just a 
small correction to the lowest order (tree-level) calculation and often their size is 
not greater than experimental uncertainties. For example, only recently have preci- 
sion experiments seen some evidence for the genuine electroweak loop corrections of 
the standard model. Despite their smallness, loop corrections in combination with 
precision data can impose important restrictions on the parameter space of various 
models. 
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The major part of my thesis studies how NHL's contribute via loops to the 
regular leptonic decays of the Z boson, a lepton universality breaking parameter 
(both studied in precision tests at LEP) and the mass of the W boson. This is a 
novel approach to the study of NHL's in loops. The previous studies concentrated 
mainly on flavour- violating decays, such as \i —>■ at very low energies or Z —>■ e + 
at the Z factories' energy. We argue that our calculations probe NHL masses and 
mixings more efficiently than these traditional studies of flavour- violating processes. 
The limits on NHL mass we obtain from the leptonic decays of the Z boson and the 
lepton universality breaking parameter are comparable to the limit derived from the 
considerations of perturbative theory breakdown, discussed in Sec. |6.7| . 

This work is organized as follows: In Chapter 2 we treat models of neutrino 
mass formally. After specifying the classical Lagrangian of the standard model, we 
investigate how the peculiar nonzero energy density of the vacuum associated with 
the existence of a fundamental scalar Higgs field accommodates fermion masses in 
the standard model. We then move on to describe neutrino masses in SU(2)l x 
U(1)y models beyond the standard model; we consider the simplest extension of 
the standard model leading to massive Dirac neutrinos and a see-saw model with 
Majorana neutrinos as examples. We also briefly describe neutrino mass in grand 
unified models and show where motivation for our model comes from. 

In Chapter 3 we discuss our model, a superstring-inspired SU(2)l x U{1)y model 
of neutrino mass, in detail. We define the fermion content and the neutrino mass 
matrix and show how massless neutrinos and NHL's arise through the diagonal- 
ization of the mass matrix in the case of one family and also in the general case 
of three families. The mixing matrix is described and phenomenologically relevant 
mixing parameters are defined. In the second part of that chapter we review existing 
constraints on NHL's. 

As a prerequisite for one-loop calculations, the standard model at the one-loop 
level is discussed formally in Chapter 4. A key ingredient of the calculations is the 
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renormalization of the standard model and we spend some time dealing with its 
salient features. 

In Chapter 5 we revisit the flavour-violating leptonic decays of the Z boson in 
our model. Our results are followed by a discussion of flavour-violating processes in 
general. 

In Chapter 6 we study the impact NHL's have through loops on flavour-conserving 
leptonic decays of the Z boson, lepton universality breaking in these decays, and the 
W boson mass. One-loop corrections are classified and calculated and the most im- 
portant diagrams are identified. Violation of the decoupling theorem by the NHL's 
and its relevance is discussed. Implicit dependence of our results on muon decay is 
clarified. 

Chapter 7 completes calculations from the previous chapter by considering the 
full set of diagrams contributing to the muon decay. 

In the last chapter we conclude by summarizing our main results. 



Chapter 2 

The models of neutrino mass 



In this chapter we treat the models of neutrino mass formally. We start with the 
standard model, then we proceed with models beyond the SM. As discussed in the 
Introduction, we concentrate mainly on minimal extensions of the standard model - 
models with SU{2) L x U{l)y gauge structure. 

In Sec. |2.1| we present the classical electroweak Lagrangian defining the standard 
model. We rather state than discuss the principles upon which it is based. The 
purpose is to specify our notation and to give a practical reference point for the 
discussions of the models of neutrino mass (Chapters 2 and 3) and for one-loop 
calculations (Chapters 5,6,7) []. 

In Sec. 12.21 we examine the fermion masses in the standard model, in Sec. 



we describe the generation of neutrino masses in simple extensions of the standard 
model (a simple model with Dirac neutrinos and a see-saw model with Majorana 
neutrinos) and finally, in Sec. |2.4] we briefly deal with neutrino masses in grand 
unified models. The superstring-inspired minimal extension of the standard model, 
motivated in this last section, is discussed in detail in Chapter 3. 



1 One- loop calculations actually require the quantization and the subsequent extension of the 
classical Lagrangian by additional terms; this is discussed in Chapter 4, which directly precedes 
one-loop calculations in Chapters 5,6 and 7. 
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2.1 Classical elect roweak Lagrangian 

The standard model of electroweak interactions (SM) || ||, f| is a gauge theory- 
based on a local SU{2) L x Z7(l)y symmetry group, which describes electromagnetic 
and weak interactions as manifestations of a single electroweak force. The SU(2)l 
part is the group of the weak isospin I and the U(l)y part is the group of the weak 
hypercharge Y. The quantum numbers I3 (the third component of I) and Y are 
related to the electric charge Q via 



The SU(2)l x U(1) Y symmetry is spontaneously broken to U(l) em (the group of the 
electric charge Q) and particle masses are generated by the nonsymmetric vacuum. 

We now list the fermion content of the SM 0. The particles, represented by 
chiral (left-handed and right-handed) fields, form three families and therefore can 
be represented as three-component vectors in family space. With the theoretical 
treatment of fermion masses in mind we differentiate between the field content of 
an unbroken electroweak theory and that of a broken one. The left-handed lepton 
fields of the unbroken theory, 



transform as doublets (I = |) under SU(2)l with the hypercharge Y = —1. In 
short, their quantum numbers are (|, — 1). 
The right-handed lepton fields (0,-2) are 



Q = h + T 



(2.1) 




(2.2) 



(e R ,fi R ,T R ). 



(2.3) 



2 In fact, we already did it in Table 1.1; here we use more formal description. 
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Note there are no right-handed neutrino fields in the SM. 

The left-handed (|, |) and the right-handed up (0, |) and down (0, — |) quark 
fields of the unbroken theory are respectively 



Ql = (Qi L ,Q2 L >Qa L ) = 
u R = (u 1r ,u 2r ,u 3r ), 

d>R = {dl R i "2 H ; ^3 R )- 



d 1 



u 2 
do 



tin 



(2.4) 



The quark fields (weak eigenstates) of the broken theory are different from the fields 
of the unbroken theory (and also from the quark mass eigenstates) : 



Ql = (qi L ,q2 L ,q3 L/ 



u 
d 



u R = (u 1r ,u 2r ,u 3r ) = (u R ,c R ,t R ) 
d R = (di R ,d2 R ,d 3R ) = (dR,SR,bR) 



(2.5) 



where d = Vckm& are weak eigenstates of the broken theory obtained from mass 
eigenstates d through the Cabibbo-Kobayashi-Maskawa (CKM) matrix Vckm (see 
also Eq. ggTp . 

The classical electroweak Lagrangian is the sum of the fermion part, the gauge 
part and the Higgs part: 



£-ew — £-g + £f + C>n- 



(2.6) 



The fermion part, which describes the fermions and their interactions, is given 
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by 



3=1 

+ ^ 7 M ^ u jR + d~j"V^ d jR } , 



(2.7) 



where 



+ fa, 



(2.8) 
(2.9) 



are covariant derivatives for left and right-handed fields respectively. These deriva- 
tives ensure the gauge invariance of the Cf by introducing a weak isospin triplet of 
gauge fields = (W^W^,WT) and a weak isospin singlet gauge field B^. The 
gauge fields interact with the fermions with the strength g 2 , the SU(2)l coupling 
constant and the strength gi, the U(l)y coupling constant. I is a 2 x 2 unit matrix in 
isospin space and ~r are generators of SU{2) L transformations in two-dimensional 
representation; r are Pauli matrices (see Appendix [A]). Weak isospin I is the eigen- 
value of the operator and J3 is the eigenvalue of the operator |t3. 

The gauge part of the Lagrangian describes the gauge fields and their self-interac- 
tions; it is given by 



G 



(2.10) 



where a = 1, 2, 3 is the 577(2) index and 



W a 



B 



[IV 



W - d u W« + g 2 e abc 
duB^, 



d^B u 



(2.11) 
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are field strength tensors for the isotriplet and the isosinglet fields respec- 
tively. 

The Higgs part of the Lagrangian, responsible for spontaneous electroweak sym- 
metry breaking, is the sum of two terms: 



£>u — £>hg + 



(2.12) 



Here Lug describes the Higgs-gauge interactions and Chf the Higgs- fermion or 
so-called Yukawa interactions. Chg has the form 



C HG = {V^)\V^)-V{<$>), 



(2.13) 



where 



$ = 



(d. + i^YB^I-i^T-W, 



V(Q) = -// 2 $ f $ + A($ t $) 2 , A>0. (2.14) 

is the covariant derivative for the Y = 1 Higgs doublet $ (|, 1) with the 
charged component + and the neutral component 0°; V(<&) is the Higgs poten- 
tial constructed so it can lead to the vacuum in which the average value (vacuum 
expectation value) of the Higgs doublet, denoted (<&), is nonzero. To keep U{l) ern 
unbroken, it is the neutral component 0° which develops the vacuum expectation 
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value: 




(2.15) 



The symmetry is broken spontaneously because the electroweak Lagrangian is sym- 
metric under SU(2) L x U(1) Y transformations but the lowest energy state, the vac- 
uum, is not (here (<£>) is not symmetric). 
The Higgs doublet can be written now as 



$ = 



f ^ ) 



(2.16) 



where ± and \ are unphysical Higgs fields and H is the physical Higgs field. 

The spontaneous symmetry breaking gives rise to massive gauge fields and 
Z^ mediators of weak charged and neutral interactions leaving the massless photon 
field Ap, the mediator of electromagnetic interactions: 



W, 



± 



z„ = 



A,, = 



+ cos 9 W W* + sin 9 W B /X , 
— sin 9 W W* + cos OwB^. 



(2.17) 



The W ± mass M w and the Z mass M z are given by 



M w = \ 2 , M z = ^JgJ+ 



gi 



(2.18) 



The Weinberg angle 9 W is defined as 



(2.19) 
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The electric charge e = yi^a can be expressed as 

0102 



01+02 



(2.20) 



or 



02 = -T-5-, 0i = 7-. (2.21) 

sm 6% cos o'vy 

The second term of the Higgs part of the Lagrangian, Chf, is discussed in the 
next section. 



2.2 Fermion masses in the SM 

The spontaneous symmetry breaking is responsible also for fermion masses. The 
starting point is Chf which describes the Yukawa interactions between fermions 
and the Higgs doublet: 



C 



HF 



»=i j=i 

3 

i=l 



+ /l.C. 



(2.22) 



where 



-r 



(2.23) 



is the Y = — 1 Higgs doublet (~, — 1) and Gy, G^-, /i, are arbitrary Yukawa couplings 
which are free parameters in the SM. The purpose of the empty space in the second 
line of Eq. |2.22| will be clarified below. 
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To generate masses, one substitutes in Eq. [2.22| the vacuum expectation value 
($) (see Eq. |2.15| ) for $. Thus the first term (plus its h.c.) in the first line of 



Eq.|2~2l gives mass to u, c, t quarks; the second term gives mass to d, s, b quarks and 



the term in the second line gives mass to charged leptons. 

Let us study the charged lepton case first. We get the electron mass m e from the 
second line of Eq. |2.22j for i = 1. The substitution of ($) yields 





1 1 _ 1 _ 

--j^hie R v e L - -j=h x e L v e R = --^hv (e R e L + e L e R ) 



= -m e (e R e L + e L e R ) = -m e ee, (2.24) 

which is the familiar form of the Dirac mass term. Without any inter-generation 
couplings in the lepton part of Eq. |2.22| (Yukawa couplings hi are simple numbers 
as opposed to matrices Gij,Gij), there are no mixings among leptons in the SM 0. 
As a result, lepton family numbers (flavours) are separately conserved and there are 
no lepton flavour-violating processes. The total lepton number L is also conserved 
since it is the sum of lepton family numbers. 

Quark masses are more involved, because inter-generation couplings are allowed 
(Gij,Gij are nondiagonal matrices in flavour space). As a consequence, Qj L ,u iR and 
d[ , the weak eigenstates of the unbroken theory, are different from the mass eigen- 



When discussing fermion masses, one cannot avoid the question of possible mixings among 
fermions. It is because we look for mass effects in various weak processes where the states of 
definite weak quantum numbers (weak interaction eigenstates) participate rather than the states 
of definite mass (mass eigenstates). Mixings then relate mass eigenstates to weak eigenstates. 
Further, mixings and masses are connected through their common origin derived from Yukawa 
couplings and the vacuum expectation value of the Higgs doublet. For these reasons we will study 
mixings along with masses. 
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states u, c, t and d, c, b. They are related through the unitary matrices Al, Ar, Bl, 
Br, 0: 



\ 



L,R 



2l,r 



A 



L,R 



I u L , R \ 

CL,R 
V *L,R J 



±l,r 



V 4,« / 



B 



L,R 



1 d L , R ^ 

SL,R 



(2.25) 



To generate quark masses we again substitute ($) for $, now in the first line of 
Eq. |2.22| . We obtain mass matrices ^Gij and ^Gij which are diagonalized by the 
matrices A and B to yield masses m u , mj, of u, b quarks: 



( m„ \ 



m c 
m t J 



V2 



B r 1 GBl 



I 



m d \ 
m s 
mj y 



. (2.26) 



Mixings arise in the charged current interactions of quarks: the quark charged cur- 
rent Lagrangian (part of Cf, Eq- p.7|) is given as 



92 

V2 



"it 
4 

v 4 y 



<^=W»{u L ,c L ,t L ) A\B Lltx 



d L ^ 



Si 

\ y 



(2.27) 



where Vckm = A\Bl is the 3x3 unitary CKM mixing matrix p2[ . It is a nondi- 
agonal matrix inducing transitions between families in charged current interactions. 
Acting on mass eigenstates d, s, b, it gives us weak eigenstates d, s, b (see Eq. ^75|) . 
There is no mixing in the neutral current Lagrangian, hence no flavour-changing 
neutral currents at the lowest order of perturbation theory, although they can arise 
at the one-loop level. 

For neutrino masses, there is an empty space in the second line of Eq. |2.22 
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because no right-handed neutrino fields vr are included. Thus, there cannot be 
nonzero neutrino masses in the SM []. 

The generation of fermion masses in the SM, as we have just described it, is 
considered to be the least satisfactory part of the SM. Each mass enters as an 
unknown parameter (Yukawa couplings are not predicted) which has to be supplied 
by experiment. The SM rather accommodates fermion masses than predicts them. 
The problem of fermion masses, and neutrino masses in particular, has been a top 
priority for particle physicists for some time now. 

2.3 Neutrino masses in SU(2)i x C/(l)y models be- 
yond the SM 

As noted in the Introduction, basic directions in the theoretical treatment of neutrino 
masses can be followed in the class of models based on the same symmetry group as 
the SM, on SU{2) L x C/(l)y. This fixes the gauge sector; the fermion content and 
the Higgs (symmetry breaking) sector offer some freedom which is used by different 
models within the class. We keep here also the symmetry breaking sector of the SM 
untouched and extend the fermion sector only. We examine two such models in this 
section and the third one, our model, in Chapter 3. 

2.3.1 A simple model of neutrino mass 

In this straightforward extension of the SM one postulates one right-handed neutrino 
field vr per family with the SU(2)l x U(1)y quantum numbers (0,0). Neutrinos are 
then treated in the same manner as all other fermions in the SM. The presence of 

4 There actually could be nonzero neutrino masses without right-handed neutrino fields if the 
Higgs sector of the SM was appropriately extended (l5) . 
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right-handed neutrino fields allows new Yukawa interactions, 

3 3 

C new = -££M£($ t ^J+/i.c.. (2.28) 

i=\ j=i 

This is the term missing from Eq. |2.22j . Neutrinos acquire Dirac mass by analogy 



with up type quarks in the SM (see Sec. j2.2|) ; the only minor difference is that here 



we do not introduce mixings among the charged leptons. Neutrino mass eigenstates 
are then different from weak eigenstates u e , u^, v T , leading to neutrino mixing and 
the violation of family lepton numbers. 

The shortcoming of this model is that it provides no answer to the problem 
of smallness of neutrino masses. We can make masses small by tuning Yukawa 
couplings hij but this is not satisfactory, as there is no good reason why the 
should themselves be small. 

2.3.2 See-saw mechanism in an SU(2)l x U(1)y model, Ma- 
jorana neutrinos and Majorana NHL's 

Charged fermions are formally described by Dirac spinors. Neutrinos are described 
in the same way in the simple extension discussed above. However, because neutrinos 
are neutral, another possibility opens up. They could be Majorana particles. To 
illustrate the difference between Dirac and Majorana neutrinos, let us decompose 
the Dirac mass term into its components, which form the so-called Majorana basis 
of a matrix representation of mass terms (see Ref. [IS], Sec. 4.5), 



muu = m {v L v R + u R u L ) = m v L v R + h.c. = -m [y L v R + v c L v c R ) + h.c. 



( 



2 \n m 



m \ I u R 

m \ v c R 



h.c. 
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v 



R 



h.c, 



\ V R 



(2.29) 



where v c = Cjqis* (C = 27 2 7°) is the charge conjugate field of v, v c R = |(1 + 7s)z/ c 
is the charge conjugate of the field v L , and v c L is the charge conjugate of the field 
vr. In the above, we used the identity v^Vr = ^l^r, proven in Appendix [A.4| . From 
Eq. [2.29| it is obvious that the Dirac mass term has a very special mass matrix in 
the Majorana basis, namely, the two diagonal terms are zero. Can we make these 
two matrix elements nonzero ? The answer is yes, if we are willing to accept the 
violation of the total lepton number L, or equivalently baryon minus lepton (B — L) 
number 0. We already broke individual lepton family numbers and there is nothing 
sacred about B — L symmetry either. 

In an SU(2)l x U(1)y see-saw model we introduce the following mass matrix 
(written for the case of one family), 



& 




+ h.c, 



(2.30) 



so the fermion content is the same as that of the simple model of Sec. |2.3.1| [] , but 
here we allow Majorana mass terms breaking B — L number conservation Q 



^Mn c L n R + h.c. 



(2.31) 



5 The relevance of B — L, rather than L, is discussed in Ref. Sec. 2.4. 

6 Note that in Eq. 2.29 we use notation isl,vr for left-handed and right-handed chiral fields 
respectively; in contrast, here we use ur rather than vr for the right-handed field. The reason 
lies in the fact that for a Dirac neutrino two independent fields vl,vr combine to form a single 
particle, while in this case i>l with its partner v R form (in the limit M 3> D) a light Majorana 
neutrino and nji with its partner n c L form a Majorana NHL; hence we use a different notation for 
the field describing a different particle. 

7 For B — L to be conserved, C mass must be invariant under the following transformat ions: 

transforms as n c L nR — > e 



e "Vl; n R 



e ta n R ; n c L 



g+M* n C. n C 



The term in Eq. 2.31 



n L n R 7^ nC L n R,7 i- e -7 h breaks B — L conservation. 
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The matrix A4 



D 



now describes two massive Majorana neutrinos rather 



D M 

than a single Dirac one. To see that, we have to diagonalize M. (e.g. Ref. [[HJ 
Sec. 5.1.4), 



/ 



M = O 



T 



m% 
m 2 



-1 
1 



O. 



(2.32) 



where m li2 = |(V M 2 + AD 2 =p M) are the masses of the two Majorana neutrinos. 



In the above, 



O 



cos 9 — sin 9 
sin cos 9 



„ 2D 
tan 20 = — , 
M 



(2.33) 



is an orthogonal rotation matrix defining massive Majorana neutrinos v , N as 



v, 



Nr 



O 



n. 



-1 
1 



o 



V R 

K n R 



(2.34) 



From here we can show 



v = v L + v R = cos 9(v L - v R ) - sin 9(n c L - n R ) = -z/ c , 

N = N L + N R = sin 9(v L + v c R ) + cos 9(n c L + n R ) = N c , (2.35) 



that is, v and N are their own charge conjugates, their own antiparticles; therefore 
they are Majorana neutrinos. 

We see how this model explains the small neutrino masses when we assume that 
M 3> D. In this limit, the masses m\ of v and vri2 of iV become 



nix 



M 



m 2 = M; 



(2.36) 
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and using sin26 l = tan 29 = 29 = =M-, we find for the weak eigenstate i>l 

vl = v L + ^N L = u' L . (2.37) 



Eq. |2.36| is the famous see-saw mass relation (see also Eq. |1.2j ), whereby a weakly 



interacting neutrino, Ul = gets very light compared to the typical family fermion 
mass D thanks to the very large Majorana mass M. Assuming D ~ m T , M has 
to be greater than about 10 8 GeV in order to meet the cosmological bound (see 



Sec. y|) m v < 25 eV. 

It looks as though we have replaced the problem of the smallness of the neutrino 
mass with another one, the problem of the big mass M. Indeed, in the context of 
an SU(2)l x U(1)y model, the origin of the big mass M is a mystery. At this point 
we invoke our motivational grounds, the unification models (see the next section). 
There are in fact large scales in these models associated with the unification energies. 
The SU{2) lxU(1)y see-saw model could be a low-energy limit of some GUT theory. 

In this thesis we are specifically interested in NHL's, described in this section by 
the field iV with the mass = M. From Eqs. |2.36| , |2.37| it is obvious that NHL's 
in this model are, first, too heavy to be observed directly in the near future, and sec- 
ond, their contribution to left-handed weak eigenstates is so small that there is little 
hope to see even their indirect effects. See-saw models tend to be phenomenologi- 
cally uninteresting. There are however models with special forms of the Dirac and 
Majorana mass matrices (in the general case of n families, masses D and M in Eq. 
p.30| become Dirac and Majorana n x n mass matrices) that avoid this suppression 



19|, |20|1 . For instance, Pilaftsis |L9| finds a relation among the elements of D and 
M matrices that leads to massless neutrinos at the tree-level and small Majorana 
masses are generated radiatively. The cosmological constraint on the scale M is 
much weaker in this model and consequently, the mixing of NHL's {Kjj ~ D/M) 
is not suppressed. We shall refer to such models as see-saw models with enhanced 
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mixings. Although calculations in this work were carried out in the context of a 
superstring-inspired model, our analysis is qualitatively valid also for this class of 
see- saw models. 



2.4 Neutrino mass in grand unified models 

Here we briefly touch the question of neutrino mass in grand unified models (GUT's). 



A nice short review of the subject can be found in Ref. |15| ; the case of SO (10) is 
discussed also in Ref. @ , and that of E 6 in Ref. |23[ . 



There are 15 chiral fermion fields per generation currently known, e^, e#, v e and 
twelve u and d quark fields. In the simplest GUT model, 577(5), these 15 fields 



are assigned to {10} and {5}-dimensional representations. There is no right-handed 
neutrino postulated; therefore, one cannot generate a Dirac mass for the neutrino 



and also it is not possible to generate Majorana mass as described in Sec. [2.3.2| . One 
can still generate Majorana mass without a right-handed neutrino if an appropriate 
Higgs field is introduced. The problem is that this Higgs field is introduced ad hoc 
and, as a result, neutrino masses do not arise in SU(5) naturally. Moreover, SU(5) 
is ruled out by the proton decay measurements pi| . 

The next popular group is 50(10). This group contains left-right symmetric 
SU(2)l x SU(2)r x SU(4:)c as its subgroup, which implies automatically the right- 
handed neutrino. The number of chiral fermion fields per generation is thus 16, filling 
the fundamental {16} representation. With a right-handed neutrino in the funda- 
mental representation, neutrino masses in 5*0(10) can arise naturally via the see-saw 
mechanism (see Sec. p.3.2|) . The actual values of neutrino masses are sensitive to the 



Majorana mass matrix M (see Eq. |2.30|) , which in turn can tell us about the partic- 



ular branch of the 50(10) breaking down to low-energy SU(3)c x SU{2)l x U{1)y- 
50(10) predicts a lower rate for proton decay than does SU(5). 

Finally, a lot of attention is paid to E 6 based GUT's |23| . This is thanks to their 
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superstring connections. Green and Schwarz []25| showed that string theory in ten 
dimensions is anomaly free for the gauge group E% x E' 8 and that the compactification 
of the additional six dimensions can result in the breaking of E% down to Eq, which 
becomes an effective GUT group. 

The fundamental representation of Eq is {27}-dimensional, implying 27 chiral 
fermion fields per generation, 11 more than we had in SO (10). These eleven fields 
must be new particles, often referred to as exotics. They include a colour triplet 
weak isosinglet quark and its antiparticle and five new leptons. Of the new leptons, 
four (two charged and two neutral) form two weak isodoublets and the fifth one is a 
weak isosinglet. 

Curiously, the superstring-inspired Eq model experiences certain difficulties in 
understanding the small neutrino masses [17, : there are no appropriate Higgs 
fields to provide the large Majorana mass M for the see-saw mass matrix (see 
Eq. [2.30| ) and therefore the see-saw mechanism does not operate here. Interesting 
solutions to this problem suggest that besides the fundamental {27}-plet there exists 
an additional, E e singlet neutral fermion field Si. At low energies, Sl, along with 
the right-handed neutrino Ur, can enrich the neutral lepton spectrum of the SM. 
The other three neutral exotic leptons decouple from the low-energy spectrum. The 
mass matrix formed by uz,, Ur and Sl offers an alternative to the see-saw mechanism 
in generating naturally light (in fact massless) neutrinos. 

The phenomenological implications of such a superstring-inspired low-energy 
model, which is just a minimal extension of the SM, are studied in this thesis. 
The model itself is described in detail in the next chapter. 



Chapter 3 



A superstring-inspired 

SU(2)i x U(l)y model of neutrino 
mass 

In this thesis we study phenomenological aspects of an SU(2)l x U (l)y model, which 
extends the neutral fermion spectrum of the SM by two new fields, the right-handed 
neutrino ur and a left-handed field Sl- The model could arise as a low-energy limit 
of a superstring-inspired E 6 GUT [[H], ; it was also suggested as a low-energy limit 
of a supersymmetry-inspired 5*0(10) GUT |2~7|] . Superstring-inspired GUT's have 
an interesting problem with neutrino masses (see the discussion in Sec. |2.4|) : the see- 
saw mechanism does not apply here and unacceptably large neutrino masses arise 
as a consequence |L7|, . The existence of the field Sl was suggested as a potential 
solution to this problem. Sl is an E Q singlet which may be present in superstring 
models. At low energies it can remain in the neutral fermion spectrum along with 
the right-handed neutrino tlr and the usual left-handed neutrino Vl- These three 
fields together with imposed B — L conservation form a mass matrix leading to an 
alternative to the see-saw mechanism in addressing the problem of the smallness of 
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neutrino masses. 

In this chapter we define the model and give a detailed treatment of neutrino 
masses and mixing matrix, and the neutrino interaction Lagrangian. 



3.1 Fermion content and mass matrix 



In this superstring-inspired model we keep, in line with introductory arguments 



in Sec. the gauge sector and the Higgs sector of the SM untouched. The 

fermion content is enlarged by two neutrino fields, ur and Sl, per family. Their 
SU{2) L x U(l)y quantum numbers are (0, 0). The field tir is a right-handed neutrino, 
while Sl is an Eq singlet neutrino field. In a single family, we thus have the following 
leptons (given with their quantum numbers): 



e 



en 



n R S L 



(3.1) 



(|,-l) (0,-2) (0,0) (0,0) 



The definition of the model is completed by specifying the mass matrix Ai. In the 
Majorana basis it is given by 



( D ^ 

D T M T 
M0 



n R 
\ S nJ 



+ h.c. (3.2) 



Each UL,nR, Sl represents now a collection of three fields, one for each family, e.g. 
i>l = (^ e ,z/ M ,z/ T ) is the vector of the three SM weak eigenstate neutrinos. D and 
M are 3x3 Dirac mass matrices. The top diagonal element must vanish unless we 
extend the symmetry breaking sector of our model. A weak isotriplet Higgs field 
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could allow this term. However, we retain the symmetry breaking sector of the SM. 
The middle element is zero due to the absence of the appropriate Higgs fields that 
would provide the Majorana mass. This is enforced by imposed B — L number 
conservation, which is also responsible for all other zeros in the mass matrix M.. 
Only terms preserving the B — L number, VZ D ur + h.c. and Sl M ur + h.c, 
remain (see footnote 7, p. 23 on B — L conservation). 

To find the physical neutrino states of the model, we have to diagonalize the mass 
matrix M.. We will do it within a single family first and then we will generalize the 
procedure for the three families. 

3.1.1 Diagonalization of A4 for a single family 

In the case of a single family, vl-i tir, Sl represent each only one field and matrices 
D, M become simple numbers. We perform the following rotation, 



n. 



\S L J 



O 



n, 



O 



\ s 'rJ 



(3.3) 



where 



O = 



1 c e s e \ 

1 
y -s c e J 



D , . 

Co = cos6, so = sin#, tan^ = — . (3.4) 



The mass matrix Ai becomes 



M = (v' L n% S' L ) O 1 



D 



DOM 
M J 



O 



n R 
K S 'nJ 



h.c. 
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{y'h n l s'l 



n c L S' L 



Dc e - Ms e 

Dc e - Ms e Ds e + Mc e 

Ds e + Mc e 
















VD 2 + M 2 
\ VD 2 + M 2 



n R 
\ S 'nJ 



+ h.c, 



+ h.c. 



(3.5) 



yielding a massless neutrino v . Moreover, we recognize the submatrix 



(n\ S' L ) 



( 







VD 2 + M 2 \ / n R 



V 



VD 2 + M 2 







Q C 



(3.6) 



as the matrix representation of a Dirac mass term (see Eq. |2.29|) in the Majorana 
basis. Indeed, putting 



rf L = Nl S^ = N C R , n R = N R , S L = N L 



(3.7) 



we reproduce Eq. |2.29| and therefore, besides the massless neutrino z/', we generate a 
Dirac neutral heavy lepton iV with the mass M = y/D 2 + M 2 . The weak eigenstate 
v L is given by 



vl = cos 9 v h + sin 9 S L 
M 



D 



VD 2 + M 2 L VD 2 + M 2 
= K L v' L + K H N L , 



(3.8) 



where Kl, Kh are mixing factors (matrices in the case of three families) for massless 
neutrinos and NHL's, respectively. 

For M 3> D the mass of the NHL, M', and the weak eigenstate ui are approxi- 
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mately 



M 



»l + — Nl, 



D 

w 



(3.9) 



that is, the mixing of NHL's is Kh = jj- 



3.1.2 Mass matrix diagonalization in case of three families 

We leave the matrix representation of the mass matrix observing that 



r 



\M = \ {uZnlS L 



( D \ 

D T M T 
MO 



+ h.c. 



X -{viDn R + n c L D T u% + n\ M T S C R + S L M n R ) + h.c. 
- {vZ D n R + VZ D n R + S L M n R + S L M n R ^ + h.c. 



vl D n R + Sl M n R + h.c. . 



(3.10) 



In the above we used the identity (for which the proof is almost identical with that 
for VjV R = VjjURi see Appendix [A.4Q 



n c L D T v c R = v L Dn R 



(3.11) 



Performing the following rotation 



G 



Sl 



Ut U 2 \l v L 
U 3 Ui I \ S L 



(3.12) 
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where G is a unitary matrix, we get 

- £ m ass = (v'l U 1 + S 'l U ?) Dn R + {^L U 2 + S' L U^ Mn R + tl.C. 

= ^(U 1 D+U 2 M)n R + h.c.+~S^(U 3 D + U i M)n R + h.c. .(3.13) 
We put (compare with D cos 9 — M sin 6 = for a single generation case) 

U 1 D + U 2 M = } (3.14) 

hence 

-£ maS s = S[(U 3 D + U 4 M)n R + h.c. = S[M'n R + h.c. (3.15) 

There is no mass term for v' L , therefore v' L is a massless neutrino. M', unlike in the 
one family case, has to be further diagonalized with rotations (Z, T unitary matrices) 
in the NHL basis: 



S l = TS l> n R = Zn R , (3.16) 

yielding 

-C mass = Sl(TU 3 DZ^ + TU 4 MZ^n' R + h.c. 

= SlM"n" R + h.c., (3.17) 

such that M" is diagonal. Identifying (see also Eq. |3~7|) 



S L = N L , 

n R = N R , (3.18) 
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we arrive at 

-C mass = NM"N = M N JT A N 4 + M N5 W 5 N 5 + M N6 W 6 N 6 . (3.19) 

Here M is diagonal with elements M^ 4 , M^ 5 , M^ 6 being masses of three Dirac 
NHL's A/4, N 5 , N 6 . The weak eigenstate vector, v L , is given by 

v L = v x = {v^v T ) = Ulu L + Uls' L = U\u L + UlT^S" L 

= K L v L + K H S" L = K L v' L + K H N L . (3.20) 



3.1.3 Discussion of mass eigenstates 



The diagonalization of the mass matrix yielded three massless neutrinos v L = 
{y'i L i u 2 L i v 'i,i) a l° n g with three Dirac NHL's iV = (JV4, N 5 , Nq) with mass ~ M. 
The masslessness of the neutrinos is the consequence of the assumed B — L symme- 
try. This symmetry also prevents neutrinos from acquiring small masses in radiative 
corrections. The neutrinos are massless due to B — L symmetry also in the SM, but 
the difference is that in the SM the B — L symmetry is an automatic consequence of 
the missing right-handed neutrino fields, while here this symmetry is imposed with 
right-handed neutrinos present. Note that massless neutrinos imply there are no 
time dependent neutrino oscillations and no neutrinoless double beta decays. 

In the light of arguments for massive neutrinos (see Sec. |1.2j ) it may seem sur- 
prising that this model yields massless neutrinos. However, small neutrino masses 
can be generated in a variant of our model by introducing a small Majorana mass 
term fi jlj], |HJ in the mass matrix M. in Eq. 



- [v L n% S L 



I D \ 

D T M T 
M (i j 



( v<L \ 



R 



\ S nJ 



+ h.c. . (3.21) 
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Besides, now that we have a superstring motivation for the field content of our model, 
there is nothing unusual about massless neutrinos either. The unnaturalness of the 
SM treatment of neutrino masses is removed, dark matter has other candidates and 
the solar neutrino puzzle has alternative explanations. Whether neutrinos have any 
mass at all may actually be of secondary interest for a theorist who is trying to come 
up with some plausible explanation of the low experimental limits on this mass. 

The weak eigenstates in our model (ui) are dominated by massless neutrinos v L 
with a small admixture (Ku ~ D/M, see Eq. |3.9| ) of NHL's N. In see-saw models, 
the NHL mixing is generally suppressed due to the scale M (see Eq. |2.36| and the 
discussion afterwards) 1 , which has to be very large to explain small neutrino masses, 
m u ~ ^jj, dictated by experiments and cosmological arguments. The NHL mixing 
in our model is not, however, restricted by the dependence of neutrino masses mj 
on scales D and M (there is no such dependence as m v i = 0). Therefore, the scale 
M can be much lower than in the case of see-saw models and hence rates for many 
interesting phenomena can be large. This means that signatures of NHL's might 
be found even at current accelerator energies and luminosities. Our model is thus 
attractive not only conceptually, but also practically. 



3.2 Properties of the mixing matrix 



The weak interaction eigenstates V\ are related to six mass eigenstates v , iV via a 
3x6 mixing matrix K with components Ki a ; I = e, /i, r and a — u x , u 2 , u 3 , N±, N 5 , N 6 
(see Eq. ^201) 



V V r ) 



( 



K 



K , 


K , 


K eNi 


KeN B 


K e N 6 ^ 


K , 


K , 








K , 


K , 

TU 3 


K T N 4 


K TNb 


K tN 6 j 



( 



Nr. 



1 An exception are the see-saw models with enhanced mixings discussed at the end of Sec. 2.3.2. 
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[K L K 



11) 



Pi 



Nr 



\ U 3 J 



N 



\N 6 j 



(3.22) 



Alternatively, we can write Q 



E {K L ) u v' iL + J2 (K H ) la N aL = (K L ) u u' iL + (K H ) la N aL . (3.23) 



i=l,2,3 



a=4,5,6 



A quick inspection tells us the matrix K has 3x6 complex parameters = 36 degrees 
of freedom. Unitarity implies an important property often used throughout this 
work, 



K L K[ + K H K ] H = l. 



(3.24) 



This property reduces the number of degrees of freedom by 9 to 27. Further elimina- 
tion of unphysical parameters via redefinition (rephasing) | of physical mass eigen- 
states leaves us 3 2 angles and (3 — l) 2 phases [p9| . This allows for possible lepton 
flavour violation, universality violation and CP violation. 

The mixing factor which typically governs flavour-conserving processes, Umix, is 
given by 



al ' 



I = e, /I, T 



(3.25) 



a=4,5,6 



and the flavour- violating mixing factor it mix is defined as 



It 



E (k*UX*h 

a=4,5,6 



al ' 



Z,/' = e,/i,r, l^t. (3.26) 



2 Where not indicated in this work, indices k run through 1, 2, 3 and a, 6, c through 4, 5, 6. 
3 This operation is also done in the SM when one parametrizes the CKM matrix. 
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Further, the following important inequality holds 

\H mix\ — ttmix I I raixi I ^ I . (3.27) 



This implies that one might observe nonstandard effects in flavour-conserving pro- 
cesses even if they are absent in flavour-violating processes. 



3.3 Interaction Lagrangians 

The charged and neutral current Lagrangians are obtained from the corresponding 
terms in the SM Lagrangian substituting for v\ from Eq. |3.23| . The charged current 
Lagrangian is given by 

C cc = -^-giW £ {EMl-75)(^ + EMl-75) 

x {K H ) la N a } + h.c. (3.28) 



and the neutral current Lagrangian as 



92 
92 



Z^UK{K H ) ial ,{l - l5 )N a + h.c. 



+ -T^Z»Y,Na{V H K H ) ahlli {l- lB )N h 
A °W a ,b 

+ ^-^E^(^)«7m(1 (3-29) 

We will also need Lagrangians with neutrinos and NHL's interacting with the 
Higgs H and with unphysical Higgs <ft + , <p~ and \. The starting point is the Yukawa 
Lagrangian describing the interactions of neutrino fields vi,nn with the Higgs dou- 
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blet $, 



C = -{vi l L )h$n R + h.c, (3.30) 

where h is a matrix (in family space) of Yukawa couplings. The Sl field does not 
couple to <3>, but rather might couple to a new SU(2) L x U(l)y singlet Higgs field 



(responsible also for mass M) present in some superstring models Pq . We do not 
introduce such a field. 

The physical Higgs part, derived in Appendix [B[ is given by 

C h = -^-'N{K^Kh)M n NH 

-^V(KlK H )M N ^^NH 



2M W y ^ n/ JV 2 
^J7(rf H Kz)M N ±^vH, (3.31) 



the unphysical neutral Higgs \ part, 



£ x = +i^N(KljK H )M N ^Nx 

+ i-^-V(K\K h )M n ^^N X 
2M W v L ' 2 A 

~ %1 2^ N{K]hKl)Mn1 ~2^ V ' X ' (3 ' 32) 

and the unphysical charged Higgs 4> + , 4>~ parts 

Cr = +^B—eZK H M N N R <p-+0{^) + h.c. (3.33) 
y 2M W M w 

Feynman rules corresponding to these Lagrangians are listed in Appendix |C|. 
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3.4 Review of existing constraints on NHL's 

Constraints on neutral heavy lepton masses and mixings come from three different 
sources. 

i) First, there is the possibility of direct production of NHL's. At e + e~ colliders 
such as LEP I or SLC, they could be produced in Z decays: [] 



(3.34) 



and subsequently decay via neutral or charged currents: 




N 




W 




/ + v . du 



The rate for Z decays into an NHL and a light neutrino has been given previously 

as 



30 



T(Z ^N a + u) 



<w(i - + ^|4) r ( z - " + ")> ( 3 - 35 ) 



where 



l=e,ii,T 



(3.36) 



In the rest of this work we drop the prime from v 
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The subsequent NHL decay rate (for < My/) is then given by 

.Mn s5 



N ^mix \ 



m. 



^r M , (3.37) 



where V ^ is the muon decay rate and is the effective number of decay channels 
available to the NHL [|31|]. LEP data effectively (better than indirect constraints, 
see below) probe NHL mixings for NHL mass up to 80 GeV [3Q . 



NHL production at pp supercolliders was studied in Ref. p2| . It was concluded 
that the CERN Large Hadron Collider (LHC) has the potential to push the limits 
on a m i X below the LEP constraints for NHL mass up to 110 GeV. 

ii) Second, there are constraints on NHL mixing parameters from a variety of 
low energy experiments and from experiments at LEP I where neutral heavy leptons 
are not directly present. NHL's however, do affect observables indirectly: due to 
unitarity properties of the mixing matrix K, a nonzero NHL mixing slightly reduces 
the couplings of light neutrinos from their SM values f\ thus affecting rates for 
nuclear f3 decays, r and ir decays, and for Z decays. The following upper limits are 
consistent with experiment [R3| 

ee mix < 0.0071 
Wmix < 0.0014 

TT mix < 0.033 (3.38) 

The limit on TT mix is improved to < 0.024 if the invisible width of the Z boson is 
included in the analysis The limits in Eq. |3.38j are model independent and 



hold for any value of the NHL mass. They arise from a global analysis of results 
including lepton universality experiments, CKM matrix unitarity tests, W mass 



5 For example in case of Wev vertex, the mixing is changed from SM value = 1 to Kl, see 
Eq. p8. 
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measurements and neutral current data from LEP I experiments. Note that the 
LEP I neutral current data analysis did not include NHL loop effects but, rather, 
only coupling constant modifications due to mixing. We consider NHL loop effects 
in this work. 

Since the limit on the parameter TT m i X plays (as the least stringent one) the most 
important role in our analysis, we will pay further attention to its source. It comes 
from the \l — t universality test based on the r leptonic decays compared to the \i 
leptonic decays. The result of the test is given as the ratio of the couplings of r and 
H to the W boson, g T j g^ (in the SM we have g T = g^ = g-i). The ratio is found from 



Setting \x\x r 



T(r^euu)/T SM (T^euu) 
T(/i — > euu)/r SM (fi — > evv) 



0, we get 



9t\ 2 
9 J 



1 



(3.39) 



with 13 



(3.40) 



© 2 



0.989 ±0.016. 



(3.41) 



iii) Finally, the NHL masses and mixings can be constrained via their contri- 
bution in loops to various processes. The calculation to the one-loop level of the 
perturbation theory is naturally more involved than the mostly tree-level considera- 
tions required for direct and indirect constraints. In return we can probe regions in 
the mixings vs NHL mass parameter space currently inaccessible to the direct and 
indirect methods. For example, as we will see, we can place upper limits on the NHL 
mass. We caution though that these limits depend on the mixings and they will be 
relaxed should tighter bounds on mixings be achieved. Still this is an improvement 
over the direct and indirect methods which are blind to NHL masses larger than 
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Figure 3.1: A one- loop diagram leading to fi — > cy decay 

M z . 

There are two classes of these processes, lepton flavour-violating and lepton 
flavour- conserving. Lepton flavour- violating decays have been considered a hot can- 
didate for a new physics manifestation in general for many decades. They include 
so far unobserved, so-called rare decays of [i and r leptons and \i — e conversion in 
nuclei (A,Z): 

fi-^&y, t -»• ej, t -»• /ry, 
/i,r^ee + e~, r — ► /ie + e~, r^efi + fi~, r^[ifi + fi~, (3.42) 

li-(A,Z)^e-{A,Z), 

and the Z boson decays 

Z-^eV^ Z^e ± r T , Z^^t t . (3.43) 

We will discuss these decays in Chapter 5; here we at least mention the decay \x —>■ &y, 
which underwent an intensive experimental scrutiny and its stringent upper limit 
places a tough constraint on the mixing parameter e\i m i X . For illustration, one of 
the diagrams contributing to fi — > ej is shown in Fig. B7TL In Sec. p33 we will show 
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that /i — > e7 gives the following upper limit on the mixing parameter ejj, m i X : 

|e/w| < 0.00024. (3.44) 

By combining the indirect constraints obtained from the global analysis (see Eq. |3.38| ) 
with the inequality relations of Eq. |3.27| one obtains the following upper limits on 
the mixing factors 

\efimix\ < 0.0032 
\/iT mix \ < 0.0068 

\er mix \ < 0.015. (3.45) 

For the mixings fiT m i X and er m i X , these are the strongest available constraints. 

The second class consists of lepton flavour-conserving processes with NHL's in 
loops. The main part of this thesis (Chapters 6 and 7) is devoted to two of these 
processes, Z — > with partial leptonic width Tu and universality breaking pa- 
rameter Ubr as observables ; and \x — > ez/eZ/^ with the W mass Mw as observable. We 
will argue that the flavour-conserving processes can be competitive with and even 
have some advantages over the flavour-violating ones. 



Chapter 4 

Standard model at the one-loop 
level 



As a prerequisite for one- loop calculations in Chapters 5,6 and 7, we discuss here 
the standard model of electroweak interactions at the one-loop level. The classical 



electroweak Lagrangian was specified in Sec. |2.1| . One-loop corrections require treat- 
ment within the framework of the quantum field theory: the classical Lagrangian 
has to be quantized and extended to include some new terms. We present these 



quantum field theoretical 'amendments' to the classical Lagrangian in Sec. [LI . 

One-loop corrections calculated from the full Lagrangian typically suffer from 
divergences. A systematic way of removing these divergences, the renormalization 



of the SM, is discussed in Sec. (4.21. There are a lot of different schemes used to 



renormalize the SM. We opted for the on-shell renormalization scheme of W. Hollik 



34] , [35fl , introduced in Sec. 
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4.1 Quantization 

The quantum field theory requires two more terms to be added to the classical 
Lagrangian, C gfix and C ghost : 

C-EW = l~-G + l~-F + l~-H + £gfix + Cghost- (4.1) 

The gauge fixing term J~- g fi x is required in order to define meaningful propagators 



of the gauge fields which are otherwise singular p6[ . The linear gauge fixing of the 
't Hooft type is given by |34| 



= ~\ (*? + F 2 z + 2F + F„) , (4.2) 



where 



Itw v 



Fz = ^{d^Z,-M z i z X 



F, = -^=d»A„ (4.3) 



and £ w ,£ z ,£ 7 are gauge fixing parameters. In the 't Hooft type gauge the vector 
boson propagators have the form (V = W, Z) 



(l _ £V\ W 



k 2 - + ie \ k 2 - £ V M 2 + ie 



and propagators of unphysical Higgs particles 4> ± , x are given by 
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The unitary gauge is denned by £ y — > oo. We can see that in this gauge unphysical 
Higgs freeze out (their propagators vanish) and only physical particles appear in 
Feynman diagrams. 

In this work we use the Feynman gauge defined by £ = £ 7 = 1. In this gauge 
there are unphysical Higgs present, but the positive trade-off is the particularly 



simple form of the gauge boson propagators of Eq. \L4 . 

The Cghost term [33], [57| is specific to nonabelian theories where the one- 
loop self-energies of the gauge bosons computed from Cg + C-f + C-u + Cgfix do not 
satisfy gauge invariance and unitarity. Cghost removes this difficulty with scalar an- 
ticommuting ghost fields m ± ,« z ,« 7 (Fadeev-Popov ghosts) which appear naturally 
in Fadeev-Popov quantization based on the path-integral method f38| . u^, u z prop- 
agators are the same as the propagators of unphysical Higgs, Eq. |47|, while the u 1 
propagator is given by 



(4.6) 



k 2 + ie 

4.2 Renormalization 

In Lew there are five independent parameters (showing only lepton Yukawa cou- 
plings hi and counting them as one): 

92, 9i, A, v, hi. 

After symmetry breaking we can replace them by an equivalent set (counting mj as 
one) 



e, M w , M z , M H , m f , 



(4.7) 
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where e 2 /47r = a and masses are those of W, Z and Higgs bosons and of a fermion, 
respectively. 

Originally these parameters were identified with their physical values (a = 
1/137, Mz = 91.137 GeV, etc.). However, loop corrections calculated in terms of 
the physical values of these parameters diverge and the parameters themselves are 
modified - their physical values are changed by an infinite amount. Renormalization 
takes care of these (so-called ultraviolet) infinities through the reexamination of the 
meaning of the Lagrangian parameters in Eq. [4.7| . We will illustrate the process on 
the electric charge e. In an effort to get to the core of the one-loop renormalization 
of the SM, we present at times simplified versions of the SM formulae. The reader 
is made aware of the simplifications in a series of footnotes. 

4.2.1 Electric charge renormalization 

The piece of Lew defining the electric charge is the interaction term of the QED 
Lagrangian 

L em = ehM*- ( 4 - 8 ) 

We identify a = e 2 /4ir with its physical value 1/137 [] measured in low-energy 
(Thomson limit k 2 = (p — q) 2 — > 0, see diagram below) electron scattering. The 
photon - charged lepton vertex corresponding to L em is given by 




= ^i377 M . 



this section we enforce a = 1/137.036 using notation (2137, e\zi 
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One- loop corrections to this vertex, calculated in terms of e^, are Q 




+ 




-*ei377/,[Ay(0) + F y (A; 2 )] + .., 



where ellipses ... represent terms with Lorentz structure different from 7 M 0, Fy is 
the form factor split off so that .Fy (0) = and Ay (0) is given by [] 



Ay(0) 



«137 /2 



4tt 



+ finite constants ) 



(4.9) 



where 2/e is the ultraviolet divergence (e — > 0) regularized by dimensional regu- 
larization (see Appendix on dimensional regularization). The one-loop corrected 
vertex F is thus given by 



r = r° + r 1 



Ay(0) — F v (k 2 



+ ..., 



(4.10) 



At low energies (Thomson limit k 2 — > 0) 



T(A; 2 -.0) = ze 137 7M[l-Ay(0)] + ..., 



(4.11) 



2 Strictly speaking, lepton self-energies, 7 self-energy and 7 — Z mixing also contribute to this 
vertex at the one-loop level. They are, however, each renormalized independently (in the on-shcll 
scheme, for example, they vanish in the Thomson limit) and they will not change the essence of 
our arguments. 

3 7/j75 & n d (p + q)^, renormalized independently 

4 This form is exact only for the diagram with 7 in the loop, the other two diagrams have the 
divergence multiplied by some combinations of . There is also an infrared divergence present in 
the former diagram that we do not deal with here. 
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hence the charge is changed from its tree- level value in T by an infinite amount: 



ei37 



ei 37 [1 - A v 



(4.12) 



To explain this difficulty, we note the quantity measured by the experiment as 



ei37 is a loop corrected charge. But the quantity on the right-hand side of Eq. |4.12| is 
a loop (one-loop) corrected charge, therefore the quantity multiplying the 1 — Ay(0) 
factor, denoted e^-?, cannot be eisj. This implies the charge in the electromagnetic 
Lagrangian cannot be identified with its physical value. The correct approach is to 
admit the independent parameters appearing in Lew are in fact 'bare', unrenormal- 
ized quantities 



e\ M w , M|, M b H , m 



f 



(4.13) 



different from the physical values. The vertex T is then given by (compare with 
Eq. [POD 



r o + r i 



1 - A y (0) - F v (k 2 ) 



+ 



1 + — ( — h finite constants ) — Fy(k 2 ^ 
47r \e 



+ 



(4.14) 



Bare parameters are unambiguously fixed by the requirement that they lead to 
correct physical values. For electric charge we demand that (in view of the discussion 
above) 



T(fc 2 0) = ie% [1 - A v (0)] + 



(4.15) 



or equivalently 



e fe [l-A y (0)] = e 137 . 



(4.16) 
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,6 



b 

e = ei 37 



"137 /2 . \ 

1 — h finite constants 

4rr Ve / 



From here (after plugging in Ay(0) calculated in terms of a b ) we can solve for e b 

(4.17) 

Using Eqs. |4.16| , Q7[ the Eq. |4.14j can be written as 

r = ie 137 ^[l-F v (k 2 )] + ..., (4.18) 

where Fv{k 2 ) is of the order 0{a\ 3 j). The infinity is thus removed from the vertex 
T , i.e., the electromagnetic vertex (and the charge) is renormalized. 

4.2.2 Renormalization schemes 



The loop calculations are rarely carried out in terms of bare parameters. A widely 

4 



used technique is to split the bare charge e b , the bare fermion mass m b f and the bare 



boson mass M b as 

e b = e + 8e, 

m b = rhf + 5m f, 

Ml = M 2 + 5M 2 , (4.19) 

where e,nif,M are renormalized (finite) charge, fermion, and gauge boson masses 
and 5e,5mf,5M 2 are infinite corrections, so-called counterterms. This split intro- 
duces a degree of freedom, as there is no unique way to perform it. Renormalized 
charge and mass can take on different finite values including the physical ones. This 
freedom leads in practice to many different ways of splitting the bare parameters, 
i.e., to many different renormalization schemes (RS). 

The difference between two renormalized charges coming from two different RS 



3 We work to order 0(ei37ai37) in Sec. 4.2. Higher order terms are neglected 
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is small, of the order 0(01137), as every renormalized charge is chosen to be equal to 
ei37 in the lowest order: 



e = ei37 [1 + 0(01137)] • 



(4.20) 



With the substitution of Eq. |4 . 1 9| , the C em becomes 



C em = e b I^IA* = eh^ + SeJj^, 



(4.21) 



where the second term is called the counterterm Lagrangian. The calculation of the 
vertex T now leads to 



f = r° + r 1 



l-A v (0) + 5e/e-F v (k 2 ) 



+ ... 



(4.22) 



The conditions of Eqs. [4.15| , |4.16| are now given by 



f(k 2 0) = ie^ [1 - A y (0) + 6e/e] 



e[l-A v (0) + Se/e] = e 137 . 



(4.23) 
(4.24) 



At this point we can illustrate two different approaches to the choice of renormal- 
ization scheme. 

i) If we prefer to use some particular value of e in the calculation, say 



e = ei 37 [1 + b a 



137J 



(4.25) 



6 Note Se/e is of the order 0(a) = 0(0137). 

7 From now on we will use T = ie^^ [l — A\/(0) — fV(fc 2 )] for the unrenormalized vertex, and T 
for the renormalized vertex, T — T+ counterterm. The counterterm contains besides Se also wave 



function renormalization factors, not considered until Sec. 4.2.4 
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the counterterm SZ = 5e/e is consequently fixed by (see Eq. [4.24]) 



5Z = 5e/e = ^I+A v (0)-1. (4.26) 



In fact, the most popular scheme in electroweak calculations is an on-shell scheme 
(OS) where 

e = e os = e 137 (6 = 0), 



5Z Ub = 5e/e Ub = A y (0) = J- + finite constants j , (4.27) 

and all masses assume their physical, on-shell values {M^ s = Mz = 91.1884 GeV, 
etc.). The Eq. |4.23| with e = ei 37 , 



f (A; 2 ^ 0) = ie 137ltl , (4.28) 

is called an on-shell renormalization condition. 

ii) A different approach is to start with fixing the counterterm Se instead of e. 
For instance we may require that 5e (5rrif,5M 2 likewise) only contain infinities (no 
finite terms) []. This is the essence of minimal subtraction (MS) scheme. One chooses 

SZ MS = Se/e MS = - = (4.29) 

' 4tt e 4tt e 

and the charge is consequently given by 

e = i MS = = e J^L (4 30) 

[1 - A y (0) + 5Z MS ] 1 + ^fx finite terms' v ' ; 

hence the electric charge e MS differs from e\^ and likewise masses mf MS ,M MS do 
not assume their on-shell values. The MS scheme is frequently used in quantum 



^Compare this with OS where 5e/e os = Ay(0) contains both infinite and finite terms. 
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chromodynamics where on-shell quark masses are not well defined anyway. 



4.2.3 Mass renormalization 



The analysis performed above for the electric charge is essentially valid also for 
masses. The difference is in the form of renormalization condition. Masses can 
be defined as the poles of the propagators. For instance, gauge boson propagators 



V = W, Z have poles at bare mass M, 
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= P C 



k 2 -M\ r +ie 



The one-loop correction to P° is given by 



po p° 




P 1 



k 2 -M y +ie 



-i S y g 



a(3j 



k 2 -M°r +ie 



For a close-up of the blob — iY^vdap (unrenormalized vector boson self-energy tensor), 
see Figs. |6.4| , |6.5| and the relevant discussion in Chapter 6. The one-loop corrected 
renormalized 9 propagator is thus (compare with Eq. f4.22| ) 



p = p° + p 1 



-ig^ v 



k 2 - M^ 2 + ie 



-*v 



k 2 - My 2 + ie. 



9 Here, we mean renormalized as far as mass is concerned. There is still one divergence remaining 
in which will be removed only after the field renormalization, see Sec. 4.2.4. Therefore we will 
withhold the notation P until then. 



54 



-ig"" 1 -ig» v 



V - M b v 2 + + \ k 2 -M v 2 + Z v (k 2 ) + 



If. 



9 (4.31) 



k 2 — My — 8My + Ey(/c 2 ) + ie 

We demand that the poles (masses) of renormalized propagators remain at their 
physical values regardless of the choice of My or 6M V (compare with Eq. 4.23 ): 



P(k 2 -> My) = s ^ = . (4.32) 

Vl k 2 -My-5My + Y, V {My)+lf W ~ M^ + if 

To put the mass on shell we have to take 

M V = M V ' S = M v , (4.33) 

or equivalent ly 

S v (M y ) -8M V = 0. (4.34) 
4.2.4 Field renormalization 

For physical S-matrix elements, the renormalization of the five parameters in Eq. [4.7| 
is all that is required. However, if one wishes to also have finite Green functions, 
then the renormalization of the fields is also required (see, e.g., Ref. f34|). 

While the particle masses are given by the poles of the propagators, the normal- 
ization of the fields is given by the residues of the propagators. For gauge boson 
fields, for example, we have: 

W^WTT, < 435 > 

with the residue equal to one (ignoring —ig ttv '). The residue (field normalization) is 
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changed by the loop corrections. To show that, we expand Ey(fc 2 ) about k 2 = M v : 



S y (A; 2 



£ y (M 2 ) + (£; 2 -M 2 )£ y (M 2 ) + ... 

t! v = d^ v /dk\ 



(4.36) 



and substitute it into the propagator of Eq. |4.31| : 



k 2 -M 2 - 5M V + £y(M 2 ) + (k 2 - My)Yly[M 2 ) + ... + ze 
Applying the on-shell condition Eqs. [4.33| , [4.34| we get 



(4.37) 



///' 



fc 2 - My + - My-)S'y(My) + ... + It 



(P-M 2 [1 + S y (M 2 ) + ...]• 



(4.38) 



with the (divergent) residue 1/[1 + £'y(My)] at /c 2 — > My. The problem is fixed by 
the field counterterms generated by the substitution 



7 1/2 V 



;i + -6z v )v,. 



(4.39) 



The field counterterm 5Z V modifies the propagator as follows 



1 



(4.40) 



1 + 5Z y (1 + 5Z V ) (k 2 — M v + ze) [1 + X'y(M 2 )] ' 

To renormalize the fields (enforce their normalization) we demand that the residues 



MyV^V 
(□ + My^ 



{l + SZv)MyV tl V, 
(l + 6Z v ){a + M v )V^. 
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of the field propagators be equal to one at the poles. This implies the following 
condition (higher order term neglected): 



t! v (A4) + 6Z v 



0. 



(4.41) 



It is easy to see from Eqs. f4.31| , |4.40| that before being put on shell, P in terms of 
the renormalized self energy Ey is given by 



XV(A; 2 



k 2 -M v + E v {k 2 )+ie 

T, v {k 2 ) - 5M V + 5Z v (k 2 - My), 



(4.42) 
(4.43) 



so that the on-shell renormalization conditions Eqs. |4.34| , |4.41| become 



Zv(My) 

d±y(M v ) 
dk 2 



0. 



(4.44) 
(4.45) 



Before we go further, one remark is in order. So far we have been discussing 
a simplified renormalization of some parameters at the one-loop level. In the next 
section we will stay at one-loop level, however, we will present the full set of coun- 
terterms and OS renormalization conditions required for the processes studied in 
this thesis. 

To prove the renormalizability of the SM, it has to be shown that the infinities 
one encounters in loop calculations to any order can be removed by the finite number 
of counterterms. This was done by 't Hooft in Refs. |BP|, 



for a general case of 



non-abelian theories with spontaneous symmetry breaking. 
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4.3 The on-shell scheme of W. Hollik 

There are many renormalization schemes used in the calculation of loop corrections 
by different authors. They are distinguished in the first place by the choice of 
independent input parameters. The choice e, M\y, Mz, Mh, ffif that we are 
using is only one of several possible. Given the set of input parameters, there are 
still infinitely many possibilities for choosing the renormalized quantities e, m. The 
OS scheme is the most popular and natural in the standard model of electroweak 
interactions. 

Even then, within the OS itself, there are many different approaches to renor- 
malization. For instance, some opt for field renormalization, others do not and those 
who do, do it with different numbers of field renormalization constants. 

In this thesis we follow the OS scheme (e, Mw, Mz, Mh, rrif) of Wolfgang Hollik 



34j |35[| . We introduce multiplicative renormalization constants for each free param- 
eter and each symmetry multiplet of fields 11 at the level of the unbroken theory: 



B„ — > ( Zo ) 2 B,, 



A* J M 

3 

9i ~* (^f ) 2 </i 



1 Multiplicative renormalization and only one constant per multiplet guarantees the gauge in- 



variance of the counterterm Lagrangian. To make a connection with Eq. 4.19, note that 

e — > Z e = (1 + SZ) e = e + <5e. 
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v -> (Z*) 5 (v-<fo) 

a -> z x (z*y 2 \ 



-> (Z*) ^'fy, (4.46) 

ten constants in all (counting Yukawa couplings as one). Five of them are associated 
with fields and five with coupling constants. 

To generate the counterterm Lagrangian SC, the renormalization constants are 
expanded as 

Z = 1 + 5Z, (4.47) 



and Eqs. |4.46| - 14.47| are applied to Lew- The counterterms added to the unrenor- 



malized quantities then yield the renormalized self-energies given in Appendix [E], 
Eq. [E.2| ; and the renormalized electromagnetic, weak neutral and charged current 
vertices given in Eq. E.3| . These renormalized expressions can be compared with 
Eq. fO| and Eq. WM - 



The ten independent counterterm constants are fixed by the nine on-shell renor- 
malization conditions Q. The first set of conditions puts the masses on-shell (com- 
pare with Eq. f!44| ) Q: 

± W (M^) = ± Z (M 2 Z ) = ± H (M 2 H ) = ± f (m)) = 0, (4.48) 

where S z , Y> H and are the W, Z, Higgs and fermion renormalized self-energies 
respectively; the second set of conditions is the generalization of the QED electric 



12 The condition on & (k) in Eq. 4.49| below fixes both 8Zl and 5Zr constants. 



13 



Only real parts of self-energies enter these conditions. The imaginary parts are finite. 
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charge renormalization: 



r ee {k 2 — o) = f{k 2 -> o) 



S^(A; 2 ^0) = 



dk? 
1 



— my 
d 



£ 7 (jfe 



fc 2 =0 



E H (k 2 " 



dk< 



k 2 =M 2 H 



(4.49) 



where £ 7 and E 7 ^ are renormalized photon self-energy and photon- Z mixing respec- 
tively. The conditions involving , f 7ee and S 7 come directly from QED. Derivative 
conditions can be compared with Eq. [4.45| derived for W and Z bosons. 

When writing down renormalization conditions, one has to be careful not to vio- 



late Ward (Slavnov - Taylor) identities jy]]. These consequences of gauge symmetry 
also relate renormalization constants to one another and can be used as a cross-check 
of the consistency of the renormalization conditions. In the set above, for example, 
Ward identities make the axial part of f 7ee vanish in the Thomson limit . 

The renormalization constants calculated from Eqs. [4.48 , [4.49| are given in Ap- 
pendix [E], Eq. E.4 . 



Chapter 5 

Lepton flavour- violating processes 



Among the processes with NHL's in the loops, the lepton flavour-violating decays 



have so far received a lot more attention [^0| ^g, [£| [44], |J, |6|, |47|, @ than 
the flavour-conserving processes [^, |50[. One of the reasons could be a certain 
preconception that the experimental signature of the flavour violation is 'much more 
dramatic'. It is our intention to show in this and the next chapter that in many 
cases this expectation is rather naive. 

Another probable reason (this time justified) is that the calculation of the flavour- 
violating processes is simpler, with the smaller number of contributing diagrams and 
without having to actually renormalize. We will demonstrate this in Sec. |5.1|in case 



of the flavour- violating decays of the Z boson. These rare processes were studied in 
the context of our model previously |M |46l ; however, the limit of large NHL mass 



was not fully investigated. This was pointed out in Ref. |47[], where the branching 
ratios for Z — > l^l^ (e // t t , e t t ) were derived in the see-saw model of Ref. 



[[19] . We therefore reexamine the flavour- violating leptonic decays of the Z boson in 
our model, carefully treating the case of a large NHL mass. The diagrams are very 
similar to the flavour- conserving leptonic decays of the Z boson discussed in Chapter 
6 and we will borrow some results from there; our intention here is to focus on the 
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typical features of the flavour- violating processes rather than on the calculational 
details. 



In Sec. |5.2| we continue with the discussion of the sensitivity of the flavour- 
violating processes in general to the presence of NHL's. We will refer to this dis- 
cussion later, when in the main part of this work - the calculation of the flavour- 
conserving processes in Chapters 6 and 7 - we will confront our results with those 
for flavour- violating processes. 



5.1 Flavour- violating leptonic decays of the Z bo- 
son 

In the SM, the CKM matrix gives rise to the flavour- violating hadronic decays of 
the Z boson at the one-loop level. In our model, by analogy, the mixing matrix K 
(see Sec. |3]2|) induces the flavour- violating decays Z — > l^l^ at the one-loop level Q. 
One-loop Feynman diagrams generating these decays are given in Fig. p.l| . There 
is no tree-level contribution since there is no mixing between the charged leptons in 
the neutral current Lagrangian. We will be studying how these graphs contribute 
to the observable, the width in particular the dependence of the width on 

parameters from the neutrino sector of our model - mixings and NHL masses. The 
analysis can be simplified (without sacrificing the salient features) by assuming the 
three NHL's are degenerate, with mass M^. 



lr This feature is not exclusive to our model. Lepton flavour-violating processes are typical for 
many other nonstandard models with mixing in the lepton sector. 
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Figure 5.1: One-loop diagrams for flavour- violating leptonic decays of the Z boson. 

5.1.1 The amplitude and the width for Z — > /{" 

The total amplitude .M is given by the sum of partial amplitudes corresponding to 
the graphs of Fig. |5TT|a-j (subscripts in M.wn--- refer to particles in the loop) 0: 



a r 

M = +iee^(l-^ 5 ) — [k 1 M WN -k 1 M Wu + k 1 M^ N + hM uuW 

+ k 3 M u NW + hMNuW + k 2 MNNW + hMwWN - kxMwWv 

+ kiM^N + faMjwN + k 2 M N N<i>}, (5.1) 

where ki,k 2 ,k 3 , k 4 are mixing factors to be derived shortly and e M is a polarization 
four-vector of the Z boson. Further, functions Mwn, ■ •■ depend on masses and 
momenta of internal and external particles. M.wn is the sum of diagrams |5.1| a, |5.1| b 
with N in the loop, Aiwu is the sum of |5.1|a, [5Tl] b with v in the loop, M.^ is the sum 



'We work in the Feynman gauge, see Sec. 4.1 



63 



of diagrams |5TT| c, [5TT| d and Ad^wN is the sum of equal contributions from Fig. |5J]h 
and |5.1| i. Besides diagrams [5.1| a, p.l| b also |5.1| f comes in both with massless z/s and 



NHL's. Diagram |5.1| e comes in with four combinations of neutral lepton types. A 
sample calculation of one function, M.NN<f>, will be given in Sec. |6.4.1| . Here we 
simply state results for all functions terms neglected) : 



M.WN 

M Wv 

■M-abW 



M 



WWa 



■M.(f>rf>N 
■M<f>WN 
■MNN<f> 



- + 4^1 



Z^4-A, + h «S, + /(*)]f, 

--^—{2M 2 z [(C 23 (M a , M w , M b ) + C n (M a , M w , M b )\ 
os w cw 

-X . ^ , fill. / ■% JT H JT It r \ a1 1 JIT Tl ,T ^ 



2 - 4Ctt(M a , M w , M b ) — A M |, a,b = N,v ; M u = 

C U (M W , M a , M w ) - C 23 (M W} M a , M w ) 



3c w r2 



C Q {M w ,M a ,M, 



Wj 



112$r x \ C ff( Mwi Mat, M, 

M 2 _ 



+ 4C£T(M W , M a , M w )-^ + A,}, 



. ivl w ) + -A M 



M 2 

-— i^C (M w ,M w ,M w ), 
M 2 



(5.2) 



where 



X 

Am 
/(*) 



sin# 



7 + In 47r + In fi , 

e 

X 2 \ogX X 
(X-l) 2 + 1-X' 



c w 



cos 8 



ir- 



(5.3) 
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Our results are written in terms of the 't Hooft-Veltman integrals Co,C 2 4,,Cn, 
C23 0| defined in the Appendix [F.3| . They contain both finite and infinite parts 
regularized by dimensional regularization [40| (see Appendix 0). Infinite parts are 



parametrized by A M , where e — > is not to be confused with the polarization four- 
vector e M . 

We now evaluate the mixing parameters. The parameter fa comes from diagrams 
with one NHL in the loop, k 2 comes from diagrams with two NHL's, fa from diagrams 
with one NHL and one massless neutrino and fa from the graph with two massless 
neutrinos. The case of one massless neutrino is trivial (- fa) and is not shown below. 
Starting with terms that come directly from the Feynman rules of Appendix 0, we 
work our way through to the final form using the properties of the mixing matrix K 
from Sec. |3.2| (if not explicitly shown, repeated indices are summed over): 



k 2 = (K H ) la (K^ H K H ) ab {K^ H ) bl , = £ (K H ) la (^ H ) am (K H ) mb (K^ H 



h = (K H ) la (K ] H ) al , = U mix , 

\ V- / ts \ / tA \ / \ /tA\ 

bl' 

m=e,fi,T 
^ m mix m ^ m ixi 

m=e,fj,,T 

fa = {K L ) u {K{K H ) ia {K ] H ) al , = £ {K L ) u {Kl) im {K H ) ma {tf H ) al , 

m=e,fi,T 

= [ 6 im - ( K H) lh (K j H ) bm ] (K H ) ma (K ] H ) al , = U' mix -Y,lm mix ml'mix 

m ' m 

= fa- fa, 

k 4 = (KLUKlK^iKi)^ = -2k 1 + k 2 . (5.4) 

For ki we show only the initial and final step. 

To address the question of infinities, we note that we do not have to actually 
renormalize. Indeed, we easily observe the mass independent divergences (in fact 
any terms independent of mass) are cancelled in the sums 



M + M WWN, 
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M NuW + M vNW + M NNW + M 

Mwu + M W n- (5.5) 

The origin of this, so-called GIM cancellation can be traced back to the unitarity of 
the mixing matrix K []. The remaining divergent amplitudes have their divergence 
multiplied by the mass term X and therefore GIM cancellation does not apply here. 
However, these divergences vanish in the sum of mass dependent diagrams, 

M HN + M^n. (5.6) 

Using Eq. [5]4] it can be shown that the width for the flavour- violating decays of 
the Z boson to lil 2 is given in terms of k\ and k 2 as 

= ^j^MzlhM. + hM^, (5.7) 

where 

Mi = M (f>WN + M M WWv + M WWN + M NuW + M vNW 

- Mv/u + MfftN + MwNi 

M2 = M]S[N(j) — MjsTuW — M„nw + M w w + Mnnw- 
The amplitude squared can be written as 

\hM 1 + k 2 M2\ 2 = \k 1 \ 2 \M 1 \ 2 + \k 2 \ 2 \M 2 \ 2 + 2Re(k 1 k* 2 M 1 M* 2 ). (5.9) 

The mixing factors ki,k 2 are process dependent and the following relations hold 

3 This cancellation is referred to as the GIM cancellation since it has a similar origin as the 
cancellations due to the CKM matrix in K° — > which lead to the postulation of the c quark 

by Glashow, Iliopoulos and Maiani (GIM) ||. 



— 2M vu w 

(5.8) 
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between CP conjugate final states: 

h,2 = K 2 {i-i+) = ki 2 (ip~), (5.10) 

implying [] 

Re {fa^lt) k^lilt) MiM* 2 ) + Re {fci^fZj) k* 2 {l+l 2 ) MiM* 2 ) 

= 2Re{k x k* 2 )Re{MxM* 2 ) , (5.11) 



giving the total rate for Z — > l±l 2 + l x I 



2 



4 (T 3 r 

r^ 2++ ^- = -^M z {|A; 1 | 2 |>l 1 | 2 + |A; 2 | 2 |A< 2 | 2 + 2 J Re(A; 1 ^) 

x i2e(A^i^)}. (5.12) 



5.1.2 Approximate relations in the limit of large NHL mass 

While we can easily see how the width F ri + +l + r depends on mixing factors, the 
dependence on the NHL mass Mjq is obscured by the algebraic complexity of the 
't Hooft - Veltman integrals. Fortunately, in the most interesting case, which is that 
of a large Mjy, the amplitudes become particularly simple. It is the most interesting 
case since the signal is the largest due to quadratic nondecoupling effects. This means 



Note the difference 



i?e{fci(^Z+) k* 2 {l^l+) MiM* 2 ) - Re{k l {l+l 2 )k* 2 {l+l 2 ) MiM* 2 ) 

= -2Im{kik^)Im(MiM* 2 ) 



may lead to a CP violating asymmetry 

r,- ,+ — r,+ § q,3 m 

1 z o (4:ir) z f z 

We found that the maximum value allowed, r\ < 2.2 x 10~ 14 for er mode at = 5 TeV, is very 



small (see experimental limits in Eq. 5.22) 
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that some of the diagrams give rise to terms 0{Mj 4 ). These effects, as well as the 
question of how high we can go with the mass Mn without disturbing perturbation 



theory, will be discussed in Sees. |6.6| and |6.7j For now, we just state that by large 
NHL mass we mean Mz < M^ < 5 TeV. 

In the limit of a large NHL mass Mat, the amplitudes exhibit the following 
behaviour: 



M-WWN 
■M-abW 



3c w 



w 



Ss w cw 



{4Cit(M w , M N , M w )-- + A,}, 

{2 - 4C£?(M a , M w , M b ) - A At }, (5.13) 



a, b 



N,u ; v,N ; N,N ; M v = 0. 



These formulae differ by less than one percent from the exact ones in Eq. |5.2| , at 
Mm = 500 GeV and the difference decreases with rising Mm to less than 0.1 percent 
at Mjy = 5000 GeV. C functions in the same limit behave as 



C {M W , M N , M w ) 

C (M Nl M w ,M N ) 
C£l n (M w ,M N ,M w ) 



M% 



1 



ln^ + 2 A /4c^- 1 [6 



9 = arctan \/4c^ 



1 + (x- 1 ) 



M% 

\-\\nM% + o(x^ 



(5.14) 



and also C2T function of any other combination of arguments involving M^ varies 
slowly as InMjy. 

With the help of Eqs. 5.13| , 5.14j , we can see there are three amplitudes in 
Eq. [O] with nondecoupling behaviour, namely the quadratic dependence on NHL 
mass. They are M.nn</>, M-^n and M^n- However, as numerical calculations show, 
■M^n ~^ —-M-^n for large M N , leaving us A4nn<p as the only amplitude with the 
nondecoupling behaviour. M.NN<j> gives the dominant contribution to M. 2 and, more- 
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over, it ensures that for large Mjv 

\k 2 \ 2 \M 2 \ 2 > N'l-Mil 2 , (5.15) 

despite the fact that the \k 2 \ is quadratically small compared to the linear \k\\. In 
Refs. p6| , |46| the authors neglected terms proportional to \k 2 \, therefore their results 
do not apply in the large Mjv limit. 



5.1.3 Numerical results 



In the numerical calculations, the term Re (kik^) from the Eq. |5.12| was treated as 



follows. We only have limits on \k±\, \k 2 \ (input parameters for our calculations), not 
on real and imaginary parts of ki,k 2 . Thus for given \ki\, \k 2 \, the real part of k\k\ 
can vary as 

-NN < Re{k x k*) < \k x \\k 2 \. (5.16) 
In our calculations we set 

Re(k 1 k*) = (5.17) 

and 5 is varied between —1 and +1 as an independent input parameter. To find a 
numerical value of the parameter |A; 2 | 2 , we express it in terms of ll m %x and hl 2m i X , 

\k 2 \ 2 = |&2| 2 (^i^) = (hhmix + hhmix) \hhmix\ 2 + \hhmix\ 2 \hh m ix\ 2 

+ 2 (hhmix + hhmix) Re {hhmix hhmix hhmix} ■ (5.18) 
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The smallness of e/j, m i X effectively removes some of the terms. For the e/i final state, 
the first and the third terms above are negligible, leaving 

\h\ 2 (efi) = \eT mix \ 2 \Tfi mix \ 2 , (5.19) 

while for er and /it sector we have 

|^2| (^t) (cCmix TTmix) \&Tmix\ 5 

|&2| 2 M = (Hfl mix +TT mix f \^T mix \ 2 . (5.20) 



The maximally allowed mixings (Eqs. |3.38| , |3.45| , 3.44) imply \k\\ = (0.00024, 



0.015,0.0068) and \k 2 \ = (0.0001, 0.0006,0.00023) for e/i, er, /it modes respectively. 
As noted before, we assume degenerate NHL's with mass M^. Gauge boson masses 



used in the numerical calculations are M z = 91.1884 GeV |52[ and M w = 80.410 GeV 
[j5"3| . The total decay width of the Z boson is taken as Y z = 2.4963 GeV . 

The results are shown in Fig. |5.2j a,b. They show how the branching ratio 
BR (ifl^j = F l + l - +l - l +/T z varies with the NHL mass. In Fig. [5~2"| a we set 5 = —1, 



in Fig. |5.2[ b S = +1. The graphs start at = 100 GeV. For NHL masses less 
than M\y, the rates are negligibly small. A sudden rise in branching ratio just 
above = 1 TeV for er and fir modes in Fig. p.2| a signals that at this point 
the |fc 2 | 2 |A / i2| 2 term overtakes the IfciplA^il 2 term and the nondecoupling behaviour 
(generated by the Feynman graph Fig. |5.1|j ) becomes dominant. We predict the 



following branching ratio limits for = 5 TeV and 5 = +1: 



BR th (Z -> eV T ) < 3.3 x 10" 8 , 
BR th (Z -> e ± r T ) < 1.4 x 10~ 6 , 

BRt h (Z -> ^t t ) < 2.2 x 10~ 7 . (5.21) 
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These results are similar to those of Ref. HSU, where the calculation was done in the 



context of a see-saw model with enhanced mixings. 

Current experimental upper limits on branching ratios are || 

BR exp (Z -> e ± /i T ) < 6 x 10~ 6 , 95% C.L., 

BR exp (Z -> e ± r T ) < 1.3 x 10~ 5 , 95% C.L., 

BR exp (Z -> ^r T ) < 1.9 x 10" 5 , 95% C.L. . (5.22) 

Our prediction for er mode is thus at least one order of magnitude below the exper- 
imental limit. 



5.2 Flavour- violating processes at very low ener- 
gies 



How does this compare with flavour-violating processes at very low energies ? The 
rare decay /i — > (see Fig. |3j]) is very well measured and supplies us with a 
stringent limit on \efi m i X \ (see Eq. 3.44 ), which we use as an input parameter for 
our calculations. We now derive this limit. The decay n — > ej was studied in the 
context of our model and see-saw models with enhanced mixings by several authors 



2D, [42|. |43]. [44], [45]]. In our model, with mass degenerate NHL's, the \x — > e-f branching 



ratio is 43 



BR(fi -> erf) 



3a 
327 



^l^mix | | Rj ( X ) | j 



(5.23) 
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Figure 5.2: The branching ratio Z ift^ as a function of M^r for (a) 5 = —1, (b) 
5 = +1 . 
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where 



is an NHL mass dependent form factor. For NHL masses Mjy > 500 GeV, which we 
ultimately consider, the formfactor becomes independent of mass, 

F 7 (X) -»■ -2. (5.25) 

This is another example of the nondecoupling behaviour (F^(X) does not vanish) 
of NHL's 0. It is the mildest nondecoupling case; we encountered in the previous 
section amplitudes with quadratic dependence on NHL mass. Given the current 
experimental limit on the \x — > ej branching ratio 0, 

BR(n -> &y) < 4.9 x 10~ u 90% C.L., (5.26) 

Eqs. |5.23| , |5.25| yield an upper limit on the mixing of \efi m i X \ < 0.00024, given 
previously as Eq. 



5 It is instructive to see how this result arises from the dimensional analysis argument. The 
effective Lagrangian for the fi — > is given by ]54| 

£ eff = T(X)e L <r\umF Xu ', 
where the field operators e^, fin, F Xl/ have mass dimensions 3/2, 3/2 and 2, respectively; hence 



T(X) has to have dimension -1. For Fig. 3.1, the large mass Mm dominance thus suggests 



T{X) ~ jj^- on dimensional grounds. However, there is also a possibility of T(X) ~ jjfi-; this is 
indeed what happens since it is (or m e ) which gives the right helicity flip to yield the required 
C e ff (these points can be best understood after writing down the amplitude for the graph). The 



amplitude for Fig. 3.1 thus decouples quadratically. 

There is another graph where two internal W's are replaced by the unphysical Higgs <p. The 
large Mm behaviour is in this case boosted by Mm dependent couplings of NHL's to (p's, so this 



dominant graph yields T(X) ~ F~(X) ~ ^M- ~ const, in agreement with Eq. 5 . 2£ 
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On the other hand, experimental limits on r — > cy and r — > ^7 



BR exp {r erf) < 1.2 x 10" 
BR exp {r -> /x 7 ) < 4.2 x 10" 



C.L. 



90% C.L. 



(5.27) 



are much weaker. The predicted rate for both ey and fxy modes is f43 



BRt h {r ->e7,^7) = 7 x 10" 7 , for Myv > 500 GeV. 



(5.28) 



However, the limits on mixing parameters used in Ref. [|^] are out of date now. 
With the current limits the predicted rate would be smaller by at least one order 
of magnitude, implying that the theoretical result is two orders of magnitude below 
the experimental upper limit for /ry mode and about three orders for cj mode. This 
explains why we had to use indirect limits of Eq. |3.45| for fiT mix and er m i X . 
Another well-measured muon decay mode is [i — > e~e~e + , with H 



BR exp (n — > e e e + ) < 1.0 x 10 



-12 



90%C.L. . 



(5.29) 



This process was considered by Refs. p0| , |42| , [P] , [45|| . The calculation shows the 
quadratic nondecoupling we encountered in the lepton flavour-violating decays of 
the Z boson. Ref. ||5| gives (with an assumption discussed therein) the following 
constraint for the parameters of the superstring-inspired (our) model: 



C^mix I &H>mix I — 0.93 X 10 



ITeV' 



(5.30) 



Mj^TeV 2 )' 

which for > 3 TeV is competitive with a constraint implied by Eqs. |3.38| , ^.45 



3.44: 



< 0.17 xlO" 5 . 



(5.31) 
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Also considered in Refs. [42], |45j is /x — e conversion in nuclei, /i (A, Z) 



e _ (A, Z). The constraint on the product ee m i X \efj, m i X \ |^5[ is similar to the one above. 

For the flavour- violating decays of the tau into three leptons (r — > e _ e _ e + , 
e~/i~/i + , etc) there is to my knowledge no calculation studying the large (TeV) 
NHL mass limit in the context of our model. Within the see-saw model of Ref . flUJl . 



Pilaftsis predicts with the current limits on mixings and for Mn = 3 TeV [flSfl : 

BRt h {r -> e~e~e + ) = 5 x 10~ 7 , 

Si2tfc(r -> e»+) = 3 x 10~ 7 . (5.32) 
The current experimental limits are [0] 



BR exp (r -»■ e- e -e+) < 1.4 x HT 5 , 90% C.L. , 
BR exp (r -> e»+) < 1.4 x 10" 5 , 90% C.L. . (5.33) 



Finally, hadronic decay modes of the r lepton, r — > /7r° are disfavoured 
by loose limits, e.g. Bi^r -> ^vr ) < 4.4 x 10~ 5 §. 

In conclusion, to probe large NHL masses, we would have to push experimental 
upper limits by at least one order of magnitude for flavour- violating leptonic decays 
of the Z boson, and by one to two orders of magnitude for flavour- violating decays of 
the t lepton. This most likely requires increased high luminosity running at LEP I 
energy and r factory f)5| . Note the dominant contribution to the total rate for 
Z — > Ifl^ + l^lt, I ^2 1 2 \-M-2 1 2 1 depends quartically on small mixings (see Eqs. |5.19 



5.2C| ) and also quartically on the NHL mass Mat. Further mass independent limits 



on mixings will therefore suppress this dominant contribution rather quickly, unless 
M N is very large. 



Chapter 6 

Lepton flavour-conserving 
processes 



In this and the following chapter we will examine two lepton flavour-conserving 
processes: i) Z — > l + l~ (I = e, /i, r) with observables Tu (the width) and U\, r (uni- 
versality breaking parameter); and ii) fi —>■ evv with observable Mw (the W boson 
mass) . We will show that these observables probe the mixings vs NHL mass param- 
eter space of our model in many respects more efficiently than the flavour-violating 
decays discussed in the previous chapter. We work to the one-loop (0(a)) level of 
perturbation theory. 

In Sec. |6.1| we classify, closely following the SM case of Ref. |34 |, one-loop cor- 



rections to Z — > into three groups - oblique, vertex and QED corrections. Each 
group is then individually studied in Sees. [T2]- We note that a large number of 
contributing diagrams comes directly from the SM without being modified by NHL's. 



In such cases we use the SM results of Ref. ||34j| . As far as non-SM contributions are 
concerned, we present a detailed calculation of two Feynman diagrams (one oblique 
and one vertex) and a summary of results for the remaining ones. Divergent results 
are then renormalized as discussed in Chapter 4. 
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Sees. |T5] - [T7| are less technical and hopefully more intriguing. We show the 
impact of loop corrections to a /z-decay on the Z width, calculate the W mass 
Mw, discuss the violation of the decoupling theorem and the quadratic dependence 
of the loop corrections on the NHL mass, and the limitations of the perturbative 
calculations. 

Our numerical results are presented in Sec. [6l| and discussed in Sec. |6i| We 
investigate here only a part of the mixings vs NHL mass parameter space by setting 
eemix = HHmix = 0. The full space is studied in Chapter 7. 



6.1 Z — ► 1 + 1 : the tree- level and the corrections 

The tree-level leptonic width of the Z boson in the SM is given by 

To = ^M z (vf + af); (6.1) 

with v\ = (—1 + Asyy) / (4:Sw c w) an d ai = —1/{AswCw) being, respectively, the 
vector and axial vector couplings of the charged leptons to Z. We neglected terms 
proportional to mf/M^. In this approximation, as a consequence of the lepton 
universality of the SM, the partial widths for all three modes (ee, fifi, rr) are equal. 

One-loop corrected leptonic decays of the Z boson in the SM were thoroughly 
discussed by W. Hollik in Ref. El. He parametrizes the leptonic width as 



i + n z (M|) 

The one-loop electroweak corrections include Z boson propagator (so-called oblique) 
corrections Ilz] vertex corrections STu and QED corrections Sqed- To give the 
reader some feeling for the numbers involved, we note that the SM prediction with 
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M z = 91.1884 GeV, m t = 176 GeV and M H = 200 GeV is 

T = 81.45 MeV, 

r H = 84.03 MeV, (6.3) 
i.e., loops account for Fu — T = 2.5 MeV. The current experimental value under the 



assumption of lepton universality is [[52 



T™ p = 83.93 ±0.14 MeV. (6.4) 
Without assuming universality 



r^P = 83.92 ±0.17 MeV, 
r«f = 83.92 ±0.23 MeV, 

T exp = 83.85 ±0.29 MeV. (6.5) 



In our model, Eqs. |6.1| , |6.2| keep the same form. It is 11^ and STu which are 
modified by the contribution of NHL's. Also, T is modified (via syy) in an indirect 
way (see Sec. |6.5| ); the QED parameter 5qed is not affected by NHL's. We now 
address these corrections one by one, starting with Sqed- 



6.2 QED corrections 



QED corrections (Fig. |6.1|) form a gauge invariant subset and therefore can be 



treated independently of the genuine electroweak corrections |34l . The graphs of 



Fig. |6.1| were calculated in Ref. |56 where the results were shown to modify the 
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a) 







Figure 6.1: QED corrections. 



Z width by a factor 5qed (see Eq. |T 



jqed 



3a 
Aix' 



Our inclusion of NHL's has no impact on this SM result. 



(6.6) 



6.3 Z-propagator corrections Ylz 

Z-propagator corrections liz are related to the real part of the renormalized Z self- 
energy E z via 



f[ z (Mf) 



dp 2 Wzh 



(6.7) 



where p is the 4-momentum of the Z boson. E z includes, besides the unrenormalized 
Z self-energy E^, through the renormalization constant 5Z^ (see Eq. E.4|) , also all 
other unrenormalized gauge boson self-energies E^,S 7 and E 7 ^. The diagrams 
contributing to these self-energies are in Figs. |6.2| 



The photon self-energy (Fig. |6.2| ) and the photon - Z mixing energy (Fig. |6.3| ) are 
not modified by the NHL's (the sum of the fermion loops runs over all fermions except 
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neutrinos) and therefore we will use the SM analytical formulae of Refs. [33], |57] given 
in Eqs. |E3J, |F^6. 



The Z self-energy (Fig. |6.4| ) and the W self-energy (Fig. |6.5|) , are modified in 



our model as NHL's enter the fermion loops. The non-SM graphs from Figs. |6.4j , pT5 
are shown explicitly in Fig. |6.6| . They include the graphs with massless neutrinos 
(no NHL's), since these differ from the SM in the mixing factors. We will calculate 
these graphs and the resulting amplitudes will replace the SM neutrino contribution 
in Eq. [E/7|. 

Y,z is associated with the transverse (<? M ") part of the Z self-energy tensor E^": 

= g^Ez+p^p'Ez. (6.8) 



The longitudinal part E^ does not contribute to S-matrix elements |58| and we will 
not consider it here. The non-SM part of the Z (unrenormalized) self-energy tensor 
is the sum of four terms corresponding to Figs. p.6| a — d : 



Ef = Ef (M N ,M N ) + E% (M w ,0) + Ef (0,M^) + Ef (0,0). (6.9) 

We evaluate the contribution of one of these terms, E^(Mjv, Mw) in detail below 
as an example. Divergent integrals are regularized and evaluated in n dimensions 
using the technique of dimensional regularization due to 't Hooft and Veltman [40 
(see Appendix 0). 

Following the Feynman rules of Appendix [C], the contribution of Fig. |6.6| a to the 
Z self-energy tensor is 

-iY%{M N ,M N ) = - E /^^{iTV^^U^a-Tfe) 



a, 6=4, 5 



z ze 




Figure 6.3: Photon-Z mixing energy. 




Figure 6.5: W boson self-energy. 
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z 



N 



a) 




b) 




c) 




d) 




e) 




f) 



Figure 6.6: Non-SM loops with NHL's and massless neutrinos. 



ia q 2 r 2 



Y,{ K H K H) ah {K\ I K H ) ba J 



d n q 
(2vr) n 



x 



Tr{ 7M (l - 75 )( k~ t> + M^t^I - 75) ( A + M N )} 



\(q-pY-M%}[ q *-M%] 



(6.10) 



In the above we sum over NHL's of all three families (a, b = 4,5,6). Using the 



relations and theorems of Appendix A.l we now evaluate the trace: 



Tr{...} = Tr{^(l- l5 )(A-^ + M N ) lu (l- l5 )(/i + M N )} 
= 2Tr{ 7M (l- 75 )( A- £>)!»{ A + M N )} 

= 2 Tr { 7m ( jq- £>) lv /ft - 2 Tr {7^ jq- §) lv A}. (6.11) 
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The trace with 75 does not contribute: 



Tr {7 M 7 5 ( h- i>)i v = U e mu/3 {q-p) a q P = 0, 



(6.12) 



usint 



ing 6n a v8 q a q 13 = and Lq a ... = p a --. ■ The original trace is thus given by 



Tr{...} = 2Tr { 7m ( jq- §) lv 

= 8[(g - v)»qv - g» v (q - p)q + q^(q - P)u}- 



(6.13) 



We plug this result back in Eq. |6.10 
-iV%(M N ,M N ) = 



0,2 2 



X 



d n g (q-p)f,qu- g^u(q-p)q + q^(q-p) l 



(27T) 



iVj 



Z7T r 



2^c 2 w ™(2tt) 
- p„B f ,(p;M N ,M N ) + B a (p;M N ,M N ) 

+ 2B^(p-M N} M N )-g^g aP B af3 (p-M N} M N )}, (6.14) 

where = S {K^ H KH) ab {K^ H KH) ba and functions B^B^ are 't Hooft scalar 

n-dimensional integrals defined in Appendix E.3| . Using (see Eq. |E. 14j ) 



B„ 



-p^Bi, B^ = PnP v B 2 \ - g^B-n, 



(6.15) 



the n-dimensional space-time relation (see Eq. P-15|) 



g^g^u = n, 



(6.16) 



and recollecting that only the transverse (terms with g^) part contributes to 
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S-matrix elements, we get 

-iE^(M N ,M N ) = - 7 r^ r k_ 



+ nB 2 2 ■ 



p 2 B 1 - p 2 B 21 - 2B : 



22 



(6.17) 



With the help of formulae from Appendix |E.3| , we arrive at the following result for 
the self-energy £z(Mat, Mn): 



i T, Z (M N , M 



n 



e 2 iir 2 1 , / 

_U It O /\/f2 i „2 



2s 2 w c 2 w " HH \2nr3 
2At(M N ) + 2M 2 N -^+ (p 2 - M 2 N )Bt(p; M N , M N )] } 



-\ ( — 3Mff+p )A 



+ 



1 k HH ^l{(-3M 2 N+P 2 )A 



+ 



2s 2 w c 2 w aa (2ir)»3 

2 

2M 2 N In M 2 N - |- + f F(p; M w , Mjv) - In M 2 



N 



x (p 2 -M^)]} 



(6.18) 



where the function F is related to the function Bq by Eq. [E.17| and AQ in (m) = 
-m 2 (-lnm 2 + 1). The divergence is displayed as a pole at n = 4 (see Appendix 



2 2 
7 — hi7r = 7 — ln7r. 

4 — n e 



(6.19) 



In n dimensions, a becomes a dimensional quantity and we should do the the fol- 
lowing replacement: 



e 2 / e 

a = — -> au e = all + -ln/i 2 + 0(e 2 ) + ... 



(6.20) 
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where \i is an arbitrary mass. Together with the expansion 
1 1 1 



(2tt)™ (2vr) 4 - e (2tt)< 

this yields 



(1 + e In 2n + 0(e 2 ) + ...), (6.21) 



- — — ct A = - — —a ( 7 + In 47r + In u 2 ) = - — —a A, 

(2tt)™ (2tt) 4 V e ' p y (2tt) 4 a 



,2 



1 2 zr 

a(A m + lnm 2 ), A m = - - 7 + ln4vr + In (6.22) 



(27r) 4 e m 2 

The self-energy thus becomes 



Z Z (M N ,M N ) = ^-^^k HH -{(-3M 2 N +p 2 )(A MN + \nM 2 N ) 
8vr >f r 'f r 6^ 

+ ■2\ri\n.\ri | 

a 1 r /p 2 2 \ p 



2M 2 N hi M 2 - ^- + (F{p; M n , M n ) - In M 2 N ) {f - M 2 N )] } 



+ ^M^M^-M 2 /)}. (6.23) 
For the other three contributions we get, following the same steps, 

+ f(p; M N ,0>(f-^-^)}, 



2 2 M| 



,p 2 M 2 N Mj 
3 6 6p 2 

E z {0,M N ) = S z (Mjv,0), 



Sz(0,0) = ^^^k LL ^(A m + F(p;m,m)-\), (6.24) 

where m 2 <C p 2 , otherwise m can be arbitrary since F(p;m,m) = 1 — ln(— p 2 /m 2 
— it) and therefore A m + F(p; m, m) is independent of m. The mixing factors 
kHH, kuL and kn can be cast into a more convenient form by converting Kl matrices 
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into Kh matrices with the help of Eq. |3.24j : 



Z HH = Y.( K HK H ) ah {K H K H ) ba = E ( K H)al( K H) lb (K H ) bj (K H ) ja 
a,b &,b,l,j 

= E {K* H ) la {K H ).{K H ) lb {K* H ) jh = ee 2 mix + \e^ mix \ 2 + \er n 2 



I ja\ M nb\ H ) jb 
a,b,l,j 



' mix 



~^~\&fJ"mix\ ^f^mix \^^~mix\ \(^Tmix\ \^Tmix\ ^"^"mixi (6.25) 



a,i a,i,j,k 

= E i K H)aM K »)ka ~ E (K H ) aj (K H ) jb (K H ) bk (KE) ka 
a,j,k a,b,j,k 



E (Ki) ak (K H ) ka - E (K H K H ) ab (K H K H ) ba 

a,k a,b 

&&mix f^f^mix 7~T m i x ^HHi (6.26) 



C LL = E ( K L K L)ji(KL K L)ij = E { K L)jk( K L)ki( K L)il{ K L)lj 

= ... = 3 -2(ee mix + fifi mix + TT mix ) + k HH . (6.27) 



The total Z self-energy in our model is obtained by cutting out the neutrino contri- 
bution from the total Z self-energy in the SM (the first line of Eq. E.7|) and replacing 
it with the sum S z (Mjv, M N ) + TL Z {M N , 0) + E z (0, 0). 

The W self-energy calculation goes along the same lines yielding 



ot r r /\% , 5 

E w (M N ,mi) = — — E llmix[—~- [P 2 ~ 7} M n - 



m 



Y1txs 1 w y , ^ "*~ L 2 V 2 " 2 



2 2p 
+ V p 2 A 1 M- m 2 m m J 3 J / 
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W l=e,[i,r " 



+ ( p - ^-^) F ^°' m ') + ^ T"]}' (6 ' 28) 

for the diagrams of Figs. |6.6| e,f respectively. The total W self-energy in our model 
is obtained by cutting out the lepton contribution from the total W self-energy in 
the SM (the first two lines of Eq. |E.8| ) and replacing it with the sum ~Ew{Mn, mi) + 
Svy(0,m;). The self-energies are then renormalized using (see Eq. [E.2| ) 

± z (f) = S z {f) - 6M 2 Z + SZi {f - M 2 Z ) , 
t w (f) = E w (f) - 6M 2 W + 8Zf (f - M&) , 



with renormalization constants given by Eq. E.4| . Note the form of the equations 
above is the same as in the SM. 

In order to better see the dependence of II ^ on M N for M N 3> M w (this is the 
limit we are ultimately interested in, see Sec. WM), we split Ilz as 



fi z = nf M + rr™ L , (6.29) 

where ITf M is the SM limit of U z and ITf HL are corrections due to NHL's. Expanding 
the F functions in powers of M^, we obtain in the limit of 3> Mw 

K N Z HL = -{^^^k HH + 0(lnM N /M^ w ) + ...}. (6.30) 

Although this formula looks very simple, one should be aware of one important fact. 
The leading term is suppressed by the mixing parameter. Indeed, the kuH mixing 
is quadratic in TT mix , while some of the O (In Mj^/M^) terms are only linear. As a 
result, a few of them are comparable in size to the leading term in M% expansion 



up to ~ 1 TeV NHL mass. We illustrate this point in Table 5.1. Here we show 
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numerical predictions for the (exact) oblique parameter tlz and compare them with 
the approximate parameter flz apP x = n| A/ + flzappx where flzappx * s the first term 
in Eq. |6.30| []• The contribution of higher order terms from Eq. |6.30 is given by 
the difference d = Tlz — Tlz ap px- Input numbers used are Mz = 91.1884 GeV, 
M H = 200 GeV, m t = 176 GeV, TT mix = 0.033 and ee mix = fifi mix = 0. 



M N 


0.5 TeV 


1 TeV 


3 TeV 


5 TeV 




frSM 
lL Z 

tlz 

TT _ ftSM , frNHL 
HZappx — i-i-z 1L Zappx 
frNHL 

Zappx 

d = Tlz — Tlzappx 


- 4.299 

- 4.313 

- 4.292 
0.007 

- 0.021 


- 4.297 

- 4.298 

- 4.270 
0.027 

- 0.028 


- 4.265 

- 4.054 

- 4.013 
0.252 

- 0.041 


- 4.199 

- 3.526 

- 3.479 
0.720 

- 0.047 


xlO" 2 
xl0~ 2 
xlO" 2 
xlO" 2 
xlO" 2 



Table 6.1: Comparison of t\ Za ppx with fl z 

For M N = 0.5 TeV, fjf^ contributes 0.007x 10~ 2 of the total difference between 
flz and flf A/ (flz - fl% M = d + flzappx)- The higher order terms contribute more 
(with opposite sign), d = -0.021 x 10~ 2 . At 1 TeV, fl^px = 0.027 x 10" 2 is 
comparable with \d\ = 0.028 x 10~ 2 and at 3 TeV it already dominates. 

Overall, flz apP x differs from flz by approximately 1% in the considered range of 
NHL masses. 



6.4 Vertex factor 5T 



u 



The vertex factor 



ST, 



-aM z \vi 



TZeF v (Ml) 



+ ai neF A {M z )}, (6.31) 



x The dependence of n| M on the NHL mass ha s its origin in a different value of the input 
parameter Mw (as calculated from G M , see Sec. 6.5) for different NHL masses. Thus the formula 
for n| M comes from the SM, but the choice of Mw comes from our model, not the SM. 
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is the source of the lepton universality breaking in r«. It includes, besides irreducible 
vertex corrections, also lepton wave function renormalization and the renormalized 
mixing energy £ 7 z- The irreducible vertices (Fig. |6.7| ) and lepton wave function 
corrections (Fig. |6.8|) are absorbed in the renormalized formfactors Fy, Fa while the 



mixing energy (Fig. 6J3) comes in as 



IMP 2 ) = ^1^1. (6.32) 



As in the case of Y*z, £ 7 z (given in Appendix |F.1[ ) includes besides £ 7 z also S 



and Yjz (through the renormalization constants bZ\ Z \bZl Z \ This is the source of 
the M N dependence of £ 7 z, in spite of E 7 ^ being a purely SM quantity. 

With T: w , Xz calculated in Sec. |T3] we get in the limit of M N 3> M w in the 
leading order an expression similar to Eq. |6.30| : 

n 7Z ,Ml, = U jzi M^-^^ kHH + 0(^) + ...}. (6.33, 



6.4.1 Irreducible vertices 

Here we examine irreducible vertices. As an example we calculate the contribution 



to the unrenormalized form factors Fy,F^ of the diagram of Fig. |6.7| f, which we 
redraw in Fig. |6.9| to show our convention of momenta flow {pi,P2,P = Pi + P2,<l) 
and charge flow (arrows on internal lines). 




Figure 6.7: Irreducible vertex corrections. 




Figure 6.8: Lepton self-energies. 
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q- p,- p 2 




Figure 6.9: Momenta and charge flow convention. 

Following the Feynman rules of Appendix |C|, the vertex V^ N( p is given by (terms 
mf/Mw neglected) 



E 

a,b 



d n q 



(2tt)" (q - vi) 2 ~ 



+m IK \ M 1 + 75 

{^H)i a M N 



+ie 
As w c w 
1 - 75 



(K^Kh)^! - 75 ) 



2 ^ — Mtv 

(6.34) 



Introducing a shorthand h = i.^H)i a {.K\ 1 K H ) ab {K ]i H ) w collecting numerical fac- 
tors at the front and merging 1 ± 75 factors we get 



X 



e 3 M% 
32s 3 w c w 2 
d n q 



4M^(l- 75 ) 



(2tt)« (g 2 - M&) [(<? - V i) 2 - M 2 ,] [(q - Pl - p 2 ) 



M 2 N ] 



(6.35) 



Now we can identify the integral as a finite C function (see Eq. |E.9| ): 



^NN<f> 



e 3 M% 



-VK 



8s 3 w c w M^ z (4 



7T r7T° 



r C (Mjv, Mw, Mat)7 m (1 - 75 
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= +te^(l-7 5 )^l 2 ^-^C (M N ,M w ,M N ) 

47T Ss^rCw M^r 
ry 

= + ie^{l-^)—l 2 M NN<t> . (6.36) 

The total contribution of the irreducible vertex diagrams and the definition of the 
unrenormalized form factors Fy , Fa is given by 



= +iej fl F v - ie^^F A 

O, r 

= +iej<*— [hM<f,wN + hM N N<t> + hM^ N - (1 - h)M W Wu 

+ hM WWN NuW vNW + UM-uvW + l 2 M NNW — 

- ie^js—^hM^wN + hM N N4> + hM^N - (1 - k)M W wu + 1\M W wn 
+ IsMnuw + hMuNw + hM V uW + hM NNW - (c| - c 3 R )M uz }, (6.37) 

where 

1 

_ _i_ 2 2 



M 



uz 



1 - • r 2 



o Q 3 „3 



2Mi(C 2 3(mi, M z , mi) + C n {m h M z , mi)) + 2 



- 4C£i n (mi,M z , mi )-A^ (6.38) 

A M is given in Eq. p.9| and M^wn, ••• were defined before (see Eq. |5T2|) . The mixing 
factors are obtained from the flavour- violating ones (see Eq. [O]) by setting I — t: 

h = {K H ) la (K ] H ) al = ll mix , 

h = {K H ) la (K ] H K H ) ab (K ] H ) bl = \le mix \ 2 + \lfi mix \ 2 + \lT mix \ 2 , 

h = {K L ) H {K{K H ) ia {K^ H ) al = h - 1 2 , 

h = (K L ) h (KlK L ) tj (Kl) 3l = l-2/ 1 + / 2 . (6.39) 
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Renormalized form factors Fy, Fa are defined (using Eq. [E.3| ) as 

= +ie^F v - iej^ 5 F A 

= +iej»{F v + Vl (5Z? - 5Zl) - (5Zf - 5Zf) + (v, 6Z l v + a t 5Z l A )} 
ie^j 5 {F A + ai (5Zf -SZ^ + ivtSZ^ + atSZy)}. (6.40) 

The check of the cancellation of divergences in shows that the divergence of 



Fig. |6.7| c is cancelled by those of Figs. |6.8| c,d and the divergence of Fig. |6.7^ by 



those of Figs. |6.8| e,f. The sum of the remaining divergences of Figs. |6.7| d,e is 
cancelled by the sum of the divergences of Figs. |6.8| a,b plus those associated with 
the counterterms 8Zf — bZ^ and 5Zj z — 5Z^ Z (these come from the bosonic loops 
of the photon-Z mixing energy). 

6.4.2 Lepton self-energies 

We define lepton self-energies S ? at the one- loop level as shown in Fig. 3.10] , i.e., 



+ = I + 3 « ^ l h = i - i Z l h, 

V V V V V V 



Figure 6.10: The definition of the lepton self-energy E' 
we call the blob i rather than — i £' (Hollik's convention). T} has the form 



= t^^-^W) + t^A^^W) + m^ l s(P% (6-41) 
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where T} VA SL R (p 2 ) are vector, axial vector, scalar, left-handed and right-handed 
parts respectively. Renormalized charged lepton self-energies T> 1 (see Eq. |E.2| ) do not 
contribute to the renormalized vector and axial vector formfactors Fv,Fa. Indeed, 
the graphs of Fig. |6.8| give zero after the on-shell renormalization conditions (see 
Eq. 4.48 ) are applied. However, the unrenormalized lepton self-energies S z do con- 
tribute to the formfactors Fv,Fa through the renormalization constants 5Zy,5Z l A . 
The relation between these renormalization constants and £' can be found in Eq. [E.4 



The charged lepton self-energy is the sum of four parts corresponding to loops 
with W — N,W — — N and Z — I respectively (see Fig. |6.8[) : 



S'(p) = + + 4 + T} ZU (6.42) 



where 



V WN = -— ^/i[l-2A M + 21nM^ + 2/(A')] i>(l- 75 ), 

6ZTTS W 

£U = -^2-(l-^i)[l-2A M + 21nM^] ^(1-75), 



32ns 



2 
u 



Z l zi = ~ 1R ^ 2 ci[l-2A^ + 21nMl] #1 - 75) 



o 



4[1-2A M + 2 In Mf] ^(1 + 75), 



= (j^^+T^ * = S" (6 ' 43) 

These expressions are evaluated at p 2 = mf as required by 5Z l v ,5Z l A (see Eq. |E.4|) , 

2 

and terms -^2- were neglected. Note the Y} zl part is a pure SM result. For complete- 
ness, we also give SL: 



= £[A M - 1 + 2fi /m (p; m x , m,) + 2£?f >; m A , m,)] ^ 
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- -[A fl -- + Bt n (p;m x ,m l )}m l . (6.44) 

7T Z 

Here m\ is the regularized photon mass and Bq,Bi (evaluated at p 2 = mf) are 
given in Eq. |E.19| . is, in this case, part of the QED subset treated independently 



of genuine electroweak corrections (see Sec. |6.2|) . However, in Chapter 7, this 
photonic correction will be included as a part of the total lepton self-energy and the 
counterterm SZ l L together with the genuine electroweak corrections. 

6.4.3 Form factors Fv,Fa in the limit of large NHL mass 

Studying the behaviour of F v , Fa in the limit of M N 3> M w , we observe in agreement 
with Sec. |5.1.2| the three sources of quadratic nondecoupling: M.NN4>-,M-<t>4>N an d 
Ai^N (which contains S^y). The quadratic nondecoupling of these amplitudes is 
also expected from dimensional analysis considerations. For illustration, for M.nn4> 
we get from Eq. |6.35| 

(q 2 - M Z N ) [(q - Pl ) 2 - Mjy] [{q - p 1 - p 2 ) 2 - M Z N \ 
Setting n=4 and neglecting all masses and momenta except Mat, we obtain 

M "» ~ M> " I w^Sw^wr (6 - 46) 

The integral has to be of the form (Mjv) fc ; the power counting yields k = —2, so 

M NN4> ~ .1/;.U V - = Ml,. (6.47) 
The exact result is (see Eqs. fT2"| , |5.14| ) 
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As mentioned in Sec. |5TT| , the contribution of J^A^n cancels with that of M^n, 
leaving A4nn<p as the only amplitude with the nondecoupling behaviour. We will 
now shed more light on this curious cancellation. The way to go is to replace in 
the diagrams corresponding to these amplitudes the Z boson with the photon and 
to use a Ward-identity [[34], which relates the vertex formfactors Fy A evaluated at 
(pi + P2) 2 = to electron self-energies represented by the counterterms 5Zv,a' 

where E 7 ,z(0) is the term originating in the bosonic loops of the 7-Z mixing. At zero 
M/v the graphs with unphysical Higgs <fi are negligible, however, with rising the 
two graphs dominate the left-hand side of Eq. |6.49 : the vertex M.^ N and the self- 
energy M1 N . Since the right-hand side of Eq. |6.49| is not affected by the NHL's, it 



remains constant and (very) small with respect to M-J^ or AiT N at M N = O(TeV) 

{~4>N — ~ JV1 4>4>N 



Hence the only way to meet the above formula is to have M^n = ~^ / Q>an i n the 



limit of large M^. 

If we now return from the photon to the Z boson, it suffices to check how the 
Feynman rules for the vertices change. It turns out that 

M+n = Ml N = cM} N , (6.50) 
M HN = M Z UN = cMlw ( 6 - 51 ) 

with the same constant c (a simple function of Weinberg angle) in both eqs., so we 
conclude that 

M^ N = -M^n- (6.52) 

Note that for the photon case the form factors Fy A are evaluated at {jp\ + P2) 2 = 0, 
while for the Z case Fy,A are calculated at (pi +P2) 2 = M\\ this is acceptable, since 
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in the limit of a large Mn the Z mass can be neglected. 

This cancellation implies that the formfactors are dominated by the single am- 
plitude M.NN<f>'- 



a 



47T 8%c w 



(6.53) 



As in the case of the oblique parameter II z, we have to carefully examine the reli- 
ability of this approximation. We note again the leading term is suppressed by the 
mixing parameter Z 2 = TT^ ix On the other hand, some of the higher order terms 
(found in many of the irreducible vertex diagrams, not just M.NN<j>) are only linear 
in TT mix , which implies they can be larger than expected on the basis of an ex- 



pansion. In Table |6.2| we show numerical predictions for the (exact) formfactor F\ 



and compare them with the approximate parameter F Vappx — ± v 

fiNHL 
Vappx 



+ £™ where 



is the leading term in Eq. |6.53| . The dependence of Fy M on the NHL mass 
has its origin in a different value of the input parameter Mw (as calculated from G^, 



see Sec. |6.5| ) for different NHL masses. Input numbers used are Mz = 91.1884 GeV, 

= 0. 



M H = 200 GeV, m t = 176 GeV, rr, 



• < •mix — 0.033 and 6& m ix 



M N 


0.5 TeV 


1 TeV 


3 TeV 


5 TeV 




fpSM 

Fv 

P _ pSM , pNHL 
r Vappx r V ' r Vappx 
pNHL 
Vappx 

d Fy Fvappx 


1.938 
2.056 
1.971 
0.033 
0.085 


1.939 
2.247 
2.071 
0.132 
0.176 


1.949 
3.525 
3.150 
1.201 
0.375 


1.971 
5.903 
5.345 
3.374 
0.558 


xl0~ 3 
xl0~ 3 
xl0~ 3 
xl0~ 3 
xl0~ 3 



Table 6.2: Comparison of Fv apP x with F\, 



For M 



N 

hSM 



0.5 TeV, F^ x contributes 0.033 x 



10 3 of the total difference 



between F v and Fy. It is less than the contribution of the higher order terms, 



2 Assuming ee m i X — [i[i m ix — and Z — > tt mode 
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d = 0.085 x 10 -3 . However, at 3 TeV the roles are switched, with FyJ^p X dominating 
as expected. 

Overall, Fy appx differs from Fy by approximately 4% at 0.5 TeV and by 9.5% 
at 5 TeV. This is worse than the 1% found for the flz parameter and it signals 
that the terms linear in mixing are more important in this case. Indeed, the second 
largest non-SM contribution to Fy comes from the graph of Fig. |6.7| d (associated 
with amplitude M.wwn)- It is given by (see Eqs. |5.13| , |5.14| ) 



a 3cw /5 
4^ 1 4s^\6 

and even at 5 TeV it is as large as 2.02 x 10~ 3 compared to F$^£ x = 3.374 x 10~ 3 . 
Other terms are also relatively large, partly cancelling the M.y/WN effect to produce 
the comparatively small difference of 9.5%. 



6.5 Imprecise Mw, precise G 



Our on-shell renormalization scheme takes a, Mz and My/ as input parameters (see 
Sec. |4.3|) . However, the direct measurement of My/ |53|| . 



My/ = 80.410 ± 0.180 GeV, (6.55) 
as opposed to that of M z |5^| . 

M z = 91.1884 ±0.0022 GeV, (6.56) 



is not yet precise enough for its use as an input parameter. To appreciate this point, 
let us examine the sensitivity of Fu to My/. It mainly comes from the tree- level 
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formula (loops are suppressed by factors of a) 



-M z ( Vl + a,) = -M z —^-2 , (6.57) 



where c w = M w jM\. It turns out that <7m w = 0.18 GeV induces <7r = or u 
= 1 MeV. This is very large compared to the experimental value (see Eq. |6. 1) 



T e * p = 83.93 ±0.14 MeV. 



As a consequence Mw, even though convenient as an input parameter for one-loop 
calculations, is usually replaced by the more precisely measured muon decay constant 

G, I : 

G M = 1.16637 (2) x 1(T 5 GeV" 2 . (6.58) 

The replacement is done in such a way that while Mw is still kept in our formulae 
for one-loop self-energies and vertices, its actual value is no longer taken from the 
direct measurement, but is rather calculated from G M . To calculate Mw from G M , 
one first computes the /i-decay rate to one-loop level in the Fermi model (this defines 
Gfj) and then equates it to the one-loop calculation of the same quantity in the SM 
34| . The result is a formula relating Mw to G^. 



VlGJl-Ai 



with Ar SM (the notation Ar is reserved for our model) representing loop effects in 
/i-decay, 



Ar SM = Ar SM (a)Mw)Mz)Mfl)m(; 



100 



a (a , 7 ~ 4s w i 2 



M, 2 



^ eE ^ ( °UC (6.60) 



Mp^ is found from Eq. |6.59| by iterative procedure. 

The reason we discuss this in some detail is that our model modifies this relation 
between and Mw both at tree and one- loop level and thus NHL's, besides con- 
tributing directly to Tu through Z-decay loops, contribute also indirectly via M w . 
This complication in the calculation of Tu is something we can actually benefit from 
since in Mw we obtained another observable sensitive to NHL's. In our model, 
Eq. |6.59| is modified as 



M w s w = ^p2Q (l Ar) _ 2 e6mix ~ 2^™^' (6.61) 

where 1 — \ee m ix — \^mix represents the tree-level modification The loop quantity 
Ar is calculated from the diagrams depicted in Figs. |T.1| - [7.4| and the diagrams with 
the corrected W propagator. We devote the next chapter to the detailed calculation 
of the Ar for the arbitrary values of the mixings U m ix- Here we make a specific 
choice of mixings which will help us to avoid most of the /z-decay loop diagrams: 
we put ee mix = fifi m ix = 0. This choice also implies e^ mix = er mix = \ir mix = 
and leaves rr mix as the only nonzero mixing parameter. Of the non-SM /i-decay 
loops, the vertex corrections, neutrino self-energy corrections and the boxes are all 
proportional to either ee mix or /i/i m j X and therefore only £^ ( see Sec. 6.3| ), the W 
propagator correction remains to modify Ar : 



Ar = ^#+*v = + (6.62) 



3 T M is modified by 1 — ee m i X ~ [i^mix and we take the square root of V ^ to get Eq. 6.61 
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In the limit of Mn ^> My/ we obtain an expression similar to Eq. |6.30| , 

Ar = Ar SM -^-^k HH ^ + 0(lnM^ w ) + ...}. (6.63) 

71 1V1 W 

The neglect of terms beyond the 0{M% j 'M^f) results in an error less than 2.5% for 
M N <5 TeV. 



6.6 Violation of the decoupling theorem 

The reader may wonder whether NHL loop effects can have any noticeable effect 
at all. After all, these loops are suppressed by factors such as ll m %x or even ll^ix 
(see for example Eq. |6.63j , where Uhh = TT mix) compared to their SM counterparts. 
Although this reasoning is correct, it is not complete. In spite of the smallness of the 
mixing parameters, the NHL loop effects can actually be larger than SM ones. The 
reason is the (possibly) large NHL mass and a violation of the decoupling theorem. 

The decoupling theorem was established and proven by Appelquist and Caraz- 
zone in Ref. ||59|| . It describes how the heavy particles of a renormalizable theory 
A enter into the low-energy theory B [^0|: All effects of the heavy particle in the 
low-energy theory B appear either as a renormalization of the coupling constants or 
else are suppressed by powers of the heavy particle mass. For instance, heavy W and 
Z bosons of the SM (theory A) decouple from the low-energy QED (theory B). 

There are, however, cases in spontaneously broken theories when the decoupling 
theorem is violated. A well-known example can be found in the SM itself with 



respect to top quark behaviour. There are two diagrams (Fig. |6.11|) which exhibit a 



quadratic dependence on the top quark mass and therefore do not vanish as m t — > oo. 



The nondecoupling effects are easily visible in this case - the diagrams of Fig. |6.11 
led to indirect bounds on the top quark mass from LEP I observables. The bounds, 
m t = 170 ± 10^x9 GeV J52|, are actually competitive with the value obtained so far 
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Figure 6.11: Diagrams with the top quark nondecoupling. 

from the direct observation at CDF at Fermilab, m t = 176 ± 8 ± 10 GeV PH . 
The NHL's also exhibit nondecoupling as can be seen in Eqs. |6.30| , |6.33| , |6.63 



where the dominant terms are quadratically dependent on M^. We note however, 
that the nondecoupling of NHL's is the consequence of our treatment of the mixings 
as the parameters independent of M N . If we replaced TT mix with TT mix ~ ^rr 

M N 

(Kh ~ the dominant terms would change from ~ rr^M^ (see Eqs. |6.30 , 

PI, |H) to 

ttLM - = g, (6.64) 

and the decoupling would be recovered. 

The analogy with the top quark can be useful in another aspect. We can make 
a naive estimate of how large an NHL mass Mn should be to produce an effect 
comparable to that of the top quark. For example, oblique corrections rise with the 
top quark mass as fif ~ m t/^-wy while in our model, the dependence on Mjv 



4 According to Okun |62[ , this point is subtle: a large positive contribution from the top quark 
is cancelled by a large negative contribution from all other virtual particles. As a result, genuine 
electroweak corrections are negligible and the bulk of ~ 2.5 MeV loop corrections is associated 
with the running coupling constant a. 
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(besides the dependence on m t ) is 

M N 2 M N / \ 

Uz ~ kHH M^ = TT ^M^> (6 ' 65) 

yielding 

m 2 t ~ rrLMj. (6.66) 

As a result, for rr mix = 0.033, = 30 m t = 5 TeV is required to produce 'the 
top size' effect. We note this is not the necessary condition for observing the NHL 
given the current mixings. As our numerical results indicate, even NHL's with the 
mass Mat ~ 3 TeV could lead to non-SM effects in IV In the same sense we do not 
have to change the top mass by 175 GeV to see conflict with the data. 

6.7 Heavy NHL's and perturbation theory break- 
down 

With one- loop corrections rising as M^/M^, we have to know at what value of 
Mtv these corrections are comparable in size to the tree-level contribution. At this 
point the theory becomes strongly interacting and the perturbative treatment fails. 
Alternatively, we can study the transition to the strongly interacting regime through 
the tree-level width of the NHL. We investigate this latter approach here and show 
the analogy to the the strongly interacting Higgs. 

There are three decay modes of an NHL, N a — > + l T ,N a — > Z + Vi and 
N a — > H + b>i, open for M N > M w , M z , M H . They are represented by the diagrams 



of Fig. |6.12| . In the limit Mm 3> Mw-, Mz, Mh, the partial decay widths are given 
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v 



Figure 6.12: Decay modes of an NHL. 



by0 

T(N a ^W ± + n 

i=l,2,3 
i=l,2,3 

where the following relations were used: 

'AKh\ I 2 



a 






16s^ 




a 




Ml 


32s^/ 


M 2 W 


a 




M% 


32s^/ 


M 2 W 



E 

i 



I 



(Kh) 1o 



The total tree-level width of the NHL is the sum of the partial widths: 

„ 3a 



N 



Mk 

-, r o "'mix „ ■ 
16% Mfr 



(6.67) 



(6.68) 



(6.69) 



5 The partial width for N a -> W ± + l T is the same as in the see-saw model of Ref . , while 
the widths for N a — > Z + v% and N a — > H + i/$ arc half of their see-saw counterparts. 
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For a comparison, the tree-level width of a very heavy Higgs is given by 



r - = ikw (670) 



Both Tn and Th rise swiftly with the particle mass and at some critical point 
they become larger than the masses Mn, Mh themselves - a clear indication of a 
perturbative breakdown. The tree-level formulae are no longer appropriate. We can 
get a safe estimate on the critical mass by demanding that 

?n,h < \m n>h . (6.71) 

With a mix < ee mix + ^ mix + rr mix = TT mix = 0.033, we obtain 

M% lt ~ 4TeV, (6.72) 

and for the Higgs we get the well-known bound of Mff lt ~ 1 TeV Q. Note that M^ lt 
is not to be interpreted as the upper bound on the NHL mass. Its sole purpose is 
to make us aware of the limitations of the perturbative treatment. 



6.8 Results 

In this section we present our numerical results. 

Our FORTRAN program is written as a modification of the routines from the 
CERN electroweak library |57| used for the SM predictions of LEP I parameters. 



This applies namely to the oblique parameter, Hz, with most of the contributing 
diagrams (see Fig. |6]^ - 16.5| ) being SM. The non-SM contributions were implemented 
as specified in Sec. jOj The vertex factor STu was derived independently of the 



3 A similar value is obtained from considerations of perturbative unitarity violation. 
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CERN electroweak library, including the SM diagrams. 

As a 'standard set' of input parameters, we used Mz = 91.1884 GeV, a -1 = 
137.036, A = = 37.281 GeV. Also part of this set are M H = 200 GeV and 

m t = 176 GeV. Below we only make it explicit when values different from the 
standard set are used, e.g., M#,m t in Figs. |6.13| , |6.17| a. To keep the discussion 
simple (without losing any qualitative features), we have reduced the free parameter 
space by assuming degenerate masses for the three NHL's. In this chapter, we 
have also imposed restrictions on the mixing parameters. We assume that ee m ix 
and /ifimix are very small relative to TT m i X . The model and NHL mass independent 
limits quoted in Eq. |3.38| are more stringent for e and /i than for r. In addition, our 
assumption is also partially [] supported by the smallness of efi m i X (see Eq. |3.44|) , as 
determined from \x — > cy, in combination with the inequality Eq. p.27| . This neglect 
of ee mix and n^mix proves useful practically in that many of the muon decay loops 
(boxes and vertex corrections, but not W oblique correction) are eliminated as a 
result. The general case of arbitrary ee m i X , ^fi m ix is considered in the next chapter. 
We present results mainly for the NHL mass range 0.5 TeV < < 5 TeV, as 
motivated by the non-decoupling arguments given in Sec. |6.6| . 

Given the relative complexity of the formulae involved, we looked for alternative 
ways of checking them, other than the special care taken in their derivation. 

First, a logical step to take is to run our program with all mixing parameters 
zero, to see if it reduces to the SM, as predicted by the CERN electroweak library. 
For the standard set of input parameters we get precise agreement, see Table |6.3| . 

Second, we checked that infinities cancelled out in all renormalized quantities. 

Third, throughout this chapter we tried to separate the dominating contributions 
in the limit of large NHL mass and represent them by simple formulae. Here, we 
collect these approximations and see, if we can understand the behaviour of partial 



7 Partially, because small e\x m i X does not necessarily imply small ee m ix and 
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frSM 




If [MeV] 


CERN library 
Our program 


-4.29769 xl0~ 2 
-4.29769 xlO" 2 


-1.16976 xlO" 3 
-1.16976 xlO" 3 


84.0297 
84.0297 



Table 6.3: SM limit of our model 
leptonic widths on their basis. Thus, collecting Eqs. |6.30| , |6.33| and p. 53] , we get 

r + 8t n 



x 



'1 + 5, 



QED) 



-pappx 
1 U 



i + n z (Mf) 

T + 5T Z M + laM z (viF$H px + aiF^^ x - viU^zJj^x) 



[1 + 6, 



QED) 



1 _L T\SM , tjNHL 
1 -|- ll z -|- tLzappx 



(6.73) 



where 



frNHL 

Zappx 



a c 



w 



' w 



M 



A' 



vr 16s^ M w 



n 



NHL 

■yZappx 



a c w m n i 



pNHL 

Vappx 



a 



1 M 2 



47T Ss-^cw M, 



N j 

14' 



(6.74) 



The third, fourth and the fifth term in the numerator in r^ pp3: are all increasingly 
negative with rising M^, therefore they make the partial width smaller. The de- 
nominator rises with Mn and also makes the partial width smaller. On the other 
hand, the tree- level result, T , whose indirect dependence on NHL's via M\y was 
discussed in Sec. |B~5|, rises with M N (see Table 



). As a result, depending on the 
particular values of the mixing parameters I2 and k^H, the partial width can either 
increase with NHL mass or decrease. An example of the latter case can be found 
in Table |6.4j which shows (for the standard set of input parameters, TT mix = 0.033 
and the above assumptions on ee mix , fiu, mix ) the Z — > rr mode partial width T rT for 
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different NHL masses. Also shown are Tq and V^P X , the prediction of Eq. |6.73 



M N 


0.5 TeV 


1 TeV 


3 TeV 


5 TeV 


T [MeV] 


81.43 


81.45 


81.82 


82.60 


T rr [MeV] 


83.99 


83.94 


83.62 


83.02 


T appx [ Me y] 


84.00 


83.97 


83.70 


83.18 



Table 6.4: To and T TT function of Mn 
Once we can predict partial leptonic widths, we can take advantage of that and 



study the universality breaking parameter defined as 49 



U, 



br 



r — r 

x tt 1 ee 




5T TT - 5T ee 


r + r 

1 tt ~ 1 ee 




2T + 5T TT + 5T ee 



(6.75) 



The new feature here is that the universality breaking parameter depends only on the 
vertex factor 6Tu] the indirect influence of Mw via T Q is suppressed and the direct 
oblique corrections cancell out in the ratio. Thus Ub r represents an independent and 
complementary quantity to the partial leptonic widths for the study of the NHL's. 
The experimental limit on the universality breaking parameter can be derived from 
the limits on partial leptonic widths as follows. From the limits on partial leptonic 
widths (see Eq. |6]^) we get 



0.9991 ±0.0040. 



(6.76) 



This ratio implies 



jjexp 
°br 



r — r 

^ tt L ee 




r TT i 


r 4- r 

1 ee T i ee 







0.00045 ±0.00200, 



leading to the upper limit (at 2a level) 



(6.77) 



U^ p < 0.00445. 



(6.78) 
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The third experimental parameter sensitive to NHL's is the W mass, Mw Under 
the assumptions made above it is only sensitive to W oblique corrections and thus 
complements the former two quantities. 

The Z leptonic widths are given as a function of NHL mass in Figs. |6.13| - 16.15 



We mostly show r rr , as the most NHL-sensitive mode under the stated assumptions 
on mixing parameters. Y ee = r w is also plotted in Figs. |6.13| a, |6.14| b . The 
dashed lines represent the la variation about the current experimental results for 
the individual Z leptonic widths V^P or Tl* p , see Eq. 

In Fig. |6.13| a, the widths T ee , T TT are shown for three values of the top quark 
mass, m t = 158, 176 and 194 GeV. The mixing was fixed at TT mix = 0.033. The 
striking difference between the two modes can be seen to arise from the competition 
between the r on one side and the 5Tz and Ilz on the other side (see the discussion 
above). For the rr mode, the mixing parameter l 2 = kuu = TT mixi while f° r the 
ee mode, kuu = T^mix an d ^2 = (see Eq. |6.39|) . Zero parameter l 2 eliminates the 
third and the fourth terms in the numerator in the second line of Eq. |6.73| and the 
rising T dominates as a result. In Fig. |6.13| b, we again set TT mix = 0.033 and show 
the Z width to t+t~ for three values of the Higgs mass, M H = 100, 200, 800 GeV. 
The dependence on the top quark mass and the weak dependence on the Higgs mass 
got transferred from the SM as expected. 



In Fig. |6.14| we vary the mixings from rr m i X = 0.02 to TT m i X = 0.07. Fig. |6.14| a 
shows r TT , Fig. |6. 14f b shows T ee . The quadratic dependence on the mixing parame- 
ter, anticipated from Eq. |6.73| , is confirmed. 

Fig. |6.15| is the only figure where we look at a 'lower energy' scale. It is a 
continuation of Fig. |6.14| a to the lower NHL mass range. As expected, the partial 
width is well within the la band and the NHL effects are negligible. 



The universality breaking factor U{, r is plotted in Fig. |6.16| as a function of Mjy, 
with the mixing parameter varied about rr m j X = 0.033. The la experimental limit 
on Ubr is indicated as the dashed line. 
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Finally, we present the NHL mass dependence of the W mass in Figs. |6.17| a,b. 



The top quark mass is varied about m t = 176 GeV in Fig. |6.17| a, while the mixing 
is held constant at rr m i X = 0.033. In Fig. p. 17 b, the mixing is varied about rr m j X = 
0.033 for a fixed top quark mass. 

6.9 Discussion and Conclusions 

Our primary consideration here has been the inclusion of neutral heavy leptons in the 
calculation of the flavour-conserving Z decays to charged leptons at one-loop level. 
The dependence of the Z leptonic widths on the NHL mass, M^, and on the mixing 
parameter rr m j X was given in Figs. |6.13j - 16.15| . We see for the experimentally allowed 



upper limit of TT mix = 0.033, and assuming ee mix = \x\i mix = 0, m t = 176 GeV, the 
Z decay width to r leptons becomes sensitive to NHL masses of about 4.3 TeV at 
the 2a level. Curiously, this is an upper limit, 

M N < 4.3 TeV, (6.79) 

rather than the lower one. The cause is the nondecoupling of the NHL's. 

Apart from this comparison of each leptonic width prediction with experiment 
we can also exploit the lepton flavour universality violation which takes place in the 
model. The universality breaking ratio, Ub r (see Fig. |6.16| ), leads to a yet better 
upper limit, 

M N < 3.8 TeV, (6.80) 

at the 2a level. The importance of Ubr is underlined by the fact that it is sensitive 
only to the vertex parameter 5Tu, unlike Tu, which besides the vertex parameter also 
depends on the Z oblique corrections and indirectly, via the W boson mass, on the 
W oblique corrections. Thus the universality breaking complements the Z leptonic 
partial widths as far as sensitivity to NHL's is concerned. 
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Figure 6.13: Z leptonic width as a function of M N in the interval 0.5 TeV < M N < 
5 TeV for (a) fixed mixing parameter TT mix , fixed Higgs mass and different values 
of m t ; both Z — > rr and Z — > ee modes shown, (b) fixed mixing parameter TT mix , 
fixed m t and different values of the Higgs mass M#; only Z — > rr mode shown. 
The dashed lines represent la band about the current experimental value = 
83.85 ± 0.29 MeV 



112 



r(M 

84.2 

84 
83.8 
83.6 
83.4 
83.2 

83 
82.8 
82.6 
82.4 
82.2 

r(M 

84.6 
84.5 
84.4 
84.3 
84.2 
84.1 
84 
83.9 
83.8 
83.7 
83.6 



eV) 




eV) 



10 

M N (TeV) 







L eV, TT mi , = 0.033 




- e + e~, TT mi 




= e + e~, TT mi 


= 0.01 


L m, = 176 GeV ^ ^ 


) 







10 

M N (TeV) 



Figure 6.14: Z leptonic width as a function of Mn in the interval 0.5 TeV < Mjv < 
5 TeV for fixed m t , fixed Higgs mass and different values of the mixing parameter, 
(a) Z — > tt mode, (b) Z — > ee mode. The dashed lines represent la band about the 
current experimental value (a)r^ = 83.85 ±0.29 MeV, (b)T e e x e p = 83.92 ±0.17 MeV. 
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Figure 6.15: Z — > rr width as a function of Mjy in the interval 50 GeV < Mn < 
500 GeV for a fixed top mass m t = 176 GeV and different values of the mixing 
parameter TT mix . The dashed lines represent the la band about the current experi- 
mental value TZ P = 83.85 ± 0.29 MeV. 
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Figure 6.16: Universality breaking parameter \J hr as a function of M N for fixed top 
quark mass (m t = 176 GeV) and different values of the mixing parameter. The 
dashed line represents la experimental limit (< 0.00245). 
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Figure 6.17: W mass as a function of Mjy for (a) fixed mixing parameter {jT mix = 
0.033) and different values of m t , (b) fixed top quark mass (m t = 176 GeV) and 
different values of the mixing parameter. The dashed lines represent la band about 
the current experimental value M w = 80.410 ± 0.180 GeV. 
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The W boson mass also exhibits some sensitivity to NHL parameters arising from 
the mixing factor modifications and the presence of one-loop diagrams containing 



NHL's, as described in Sec. |6.5| . From Figs. |6.17| a,b we see that the sensitivity of the 
W mass, currently measured as M w = 80.410 ± 0.180 GeV f>3|, to NHL's depends 
to a large degree on the top mass. The experimental error on Mw might be expected 
to come down to about 0.05 GeV once LEP II measures W pair production |63| . 



In conclusion to this chapter, the Z boson flavour-conserving leptonic widths 
along with the lepton universality breaking parameter and the W mass represent 
some of the best quantities sensitive to the NHL mass. This applies especially to 
NHL's that are far too heavy to be produced directly at present colliders. The 
only way to probe their mass in this case is via their loop contributions. Much of 
the previous effort on NHL studies has been so far concentrated on flavour-violating 
processes, either at very low energies (//, r decays, see Eq. |3.42| ), or at Z-peak energy 



[Z decays, see Eq. |3.43| and Sec. |5.1|) . We feel there are at least two reasons which 



give the processes studied in this thesis a distinct advantage over the flavour- violating 
decays. 

First, we were able to actually reduce the allowed region in the mixings - NHL 
mass parameter space here (see Eqs. |6.79 , |6.80|) using the current experimental 



data. The only flavour- violating process that competes with the limits of Eqs. |6.79| , 
3.8C is the decay fi — > eee, which sets an upper limit on NHL mass (see Eqs. |5.29 



5.30| , assuming ee m i X = 0.0071, \e\i m ix\ = 0.00024) at 2a level 

M N < 3 TeV. (6.81) 

The flavour-violating decay rates for r and Z are below the current experimental 
sensitivity (see Sees. p.l.3| , |5T2| ). 

Second, the inequality Eq. |3.27| can further suppress the flavour-violating pro- 



cesses against the flavour-conserving ones via the 'conspiracy of the phases' in the 
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sum of complex terms making up the flavour- violating parameters ll' mix - In this case, 
flavour- conserving processes are the only way to probe very heavy NHL's. 



Chapter 7 



Muon decay and W mass 



Here we generalize our analysis from the previous chapter by considering the case of 
arbitrary mixings ee m i X , fifimix and TT mix . In this case non-SM box, vertex and self- 



energy diagrams contributing to the muon decay (see Figs. - |7.4j ) may become 
important for the calculation of M^y. This is in contrast to the previous chapter 
where, as a result of the assumption ee m ix = H^mix = 0, only oblique corrections 
(corrections to the W propagator) had to be considered (see Sec. |6.5|) . Still, we 
assume in this chapter efi mix = and fiT mix = 0. The last parameter, \iT mix does 
not contribute to the muon decay. The neglect of e\i m ix is supported by experiment 
(see Eq. j3.44|) and the neglect of fir mix is motivated by our intention to keep the 



discussion as straightforward as possible. 



7.1 Box diagrams 



The set of box diagrams contributing to the muon decay is depicted in Fig. |77 
As an example, we will calculate the diagrams of Fig. |7.1| e. These two were 

chosen since, as we will show, they exhibit quadratic nondecoupling. 

The amplitude for the diagram with the Higgs boson H is given by (we sum over 
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Figure 7.1: Box diagrams for muon decay 

NHL's iV a , iV;, with M Na = M Nb = M N and we neglect external momenta in the 
internal propagators) 

M * NHN ~ jfj J2^ U »—M^ {KlKh) - 2 ^-M n ^M^ {Kh) ^ 

X (1- 75 K -e^^(^Ul+75)^^^^(^) b , 



k 2 -M^k 2 - Mfj 



9 > ^Ei^KHU^jK^JlClK^f ^ 



x 



128M^f^ v L >iay Hla ^ >eby H Jbj J (2tt) 4 
[u^(l + 7s)(^ + M jV )(l - TsMgl + 7s) (/: + M jV )(l - 75KJ 

g\M% r dk 4 [iz^(l + 75) ^][^(1 + 75) MJ n 



32 (2tt) 4 (jfe2 _ M 2 ) 2 (k 2 - M^)(k 2 - M 2 H ) ' 



where 



k mix = (KlK H )jKi) afi (K H )jK^ H K L ) bj = (Kl)^(K L ) ej ee mix ^ r , 



Using / k € k p = jg eu J k 2 ... we get 

o 4 M 4 

M^nhn = ^j^k mix [u^(l + 75)7^] [^(1 + 75)7^. 

where 2s (to) is the integral 

(4tt) 2 /■ c/ 4 fc k 2 



J 2 (m) = — — / 



i ^ (2vr) 4 (P-M2 r ) 2 (P-M2 / )(P-m 2 )' 

In the next step we use the tree-level amplitude 

■ 2 

and (?2 = e/% to obtain 

■M-cbNHN = —^—ee m i x fX fJ, m i x „ — ^- J2 (Mg)M tree . 

The integral X 2 (M H ) is calculated in Appendix 0. The result is 

± 2 {M H ) - — — T < -5- + 



M| It, M i 



The amplitude M^nxN is obtained from M.^nhn by replacing M# with M/. 



121 



The total contribution of the box diagrams (Figs. [O] a-g) is 



M 



box 



= -M-ZeWfi + M-WuZu + M-WvZN + M-WNZu + M-WNZN 
+ -M-^NZN + M-WNHN + A^jVxJV + M-^NHN + -i^<f>NxN 



1 



-Ml 



1 2 

4 (~2 + S w) lo+M 1 ~ Wmix)( l ~ ee mix) 



47T [ 4sfyCfy 

X 1 {M Z ){1 - nn mix )ee mix + T x {Mz)\x\i mix {\ - ee mix ) + 2 2 (M z )ee 



1 



-^1 3 (M Z )+I 3 (M H )+1 3 (M Z ) - —^I^Mh 



"W 



4M 2 , 



4M 2 , 



X 2 (M 2 



C&mixH'H'mix ~r 



2( 2 + s wy „ 



„2 2 



-M, 



-^o(l &&mix) 



X l (M z )ee mix + J (l - fc) + T x (M z )nn mi ^ 



(7.8) 



where the integrals T Q ,Ii(m),T 3 (m) are 



2n 



Ji(m) 



(47T) 



In 



i J (27ry k*(k* - M^ik* - Ml) Ml- Ml Ml 



(47T) 



1 



i ^ (27r) 4 (^-M 2 r )(A ; 2 -M^)(fc 2 -m 2 ) 
1 



m 2 - Af^ 



+ 



M 2 M 2 
In 



2 M 2 , - M 2 Ml m 2 -Ml Ml 



N 



N 1V± N 



(7.9) 



(7.10) 



I 3 (m) 



(47T) 



z 7 (2vr) 4 (A; 2 -M 2 r ) 2 (A; 2 -M 2 ,)(A; 2 -m 2 ) 



+ 



M^ln^l 



m 2 -M 2 v \M 2 N -M 2 v (Ml-M 2 ,) 2 M 2 N - m 2 
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m 2 In 

(M% - w?) 



(7.11) 



For the calculation of these integrals see Appendix [B|. In the limit of 3> Mw, 
the leading contribution to X 2 (m) is very simple, 

ir px (m) = -± (7.12) 



implying the Mjy dependence of the J^A^nhn and M. 



<j>N X N, 



(7.13) 

47T lbSw Jllm 



The leading contribution to 1^{m) goes as ~ V-^Kzv and to T\{m) as ~ 1/Mjy-, 
therefore no other box depends quadratically on the NHL mass. 
The last box diagram, Fig. |7.1| h can be written as [ 3"4" |: 



M ieW » = A^tree— (in— + ln— -21n^-21n^ + 77) + ... 7.14 

47r v m e A A 2 7 

The ellipses denote additional terms discussed in Ref. 



7.2 Neutrino self-energy and its renormalization 

One half of the neutrino self-energy diagrams contributing to muon decay is shown in 
Fig. [7.2| . The corresponding self-energy is denoted as £ Vm . The other half consists of 
the same loops sitting on the bottom neutrino leg with the corresponding self-energy 
Tj Ve . In all these diagrams, we sum over the internal massless neutrinos z/^, k = 1, 2, 3. 
In principle, the graphs with z/^ replaced by N a are also present, however, they are 
suppressed by the large mass Mjy. 

Without derivation, we present results for the unrenormalized neutrino self- 
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Figure 7.2: Neutrino self-energy diagrams for muon decay 



energy H Vl (I = e, //). It has the form 



57* = -E2 ^(1-75), 



(7.15) 



where E L ' receives the following contributions from the diagrams of Fig. [T2 



E£ = Ef(p)+E^(p)+Ef^) + Ef'»+E 



+ 



+ 



+ . n n II 



M\ 



^ n (p;M H ,M N ) + B{ in {p;M H , 



Mi 



N 



■N 



2 A m + 5 o M 



1 ri 1 

/i 2" 2 v ""mix) ~ 



1 



2 ^ 



i + £ /in (p; M w , m, 0) + >; Af w , m t 0)J |. (7.16) 



(p;M z ,M N ) + B{ m (p;M z ,M N 
+ B{ in (p;M Z} M 1 
~A„ - i + 2Jjf*(p; M z , 0) + 5f >; M z , 0) 



Here s = = < M 2 H) M 2 Z) M^, M 2 Nl therefore Eqs. |QB| apply for 5^ m , fif 1 ". 
The amplitude for the diagrams of Fig. [7.2| in terms of E^ can be shown to be 
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equal to 



M se lf = -Mtree^i (7.17) 



where the one half comes from our dealing with the external wave function rather 
than the neutrino propagator. 

7.2.1 Renormalization of the neutrino self-energies 

Let us now investigate the question of the renormalization of the diagrams of Fig. \l.2\ . 
It is the only case in this work when the counterterms are modified from their SM 
form []. The problem is how to renormalize a part of a theory where interaction 
eigenstates are different from mass eigenstates. Curiously, this also happens in the 
SM quark sector. The difference is that in the SM the problem is circumvented by 
arguing the off-diagonal quark mixings are too small (with the required accuracy in 
mind) to have any effect in the loops and the renormalization procedure is effectively 
simplified to that of mass eigenstates being also the flavour eigenstates. In our model, 
we cannot neglect the off-diagonal mixings (U m i X ), since they (in combination with 
TeV NHL masses) lead to the dominant terms in the predicted signals Q. 

To derive the required counterterm, we start with the counterterm Lagrangian, 
which has the same form in both SM and our model: 

i 5Z e L ir e jdv e + i SZ^rjdu^ + i 5Z T L ir jdv T , (7.18) 

where SZ l L is the sum of 5Z l v and 5Z l A renormalization constants which are given 



1 So far we have used the SM form of the counterterms, see Appendix E.l. The actual value of 
the counterterms was, of course, different from SM. 

2 Just a few days before the submission of this thesis an interesting paper appeared on the 
preprint bulletin board |64| ], where the arguments of this paragraph are also made. The authors 
then continue to develop the first formal framework ever for the renormalization of theories with 
nonncgligiblc mixings between mass and interaction eigenstates. Their treatment is more general 
than ours below. 
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in Eq. |K.4| . Weak eigenstates v\ can be expressed in our model in terms of mass 
eigenstates v^N a (see Eq. |3.23[ ) 



E {K L ) u Ui + £ (K H ) la N a 



£ K lna n a , 

n a =vi,...,Ne 



vi = n a K ] nal . 



(7.19) 



This gives us for the product viv\ 



v\Vi = nJ\'J\n 



K^Kn = £ Vi(Kl) il (K L ) lh v k + ...(ViN t Nv k ,NN), (7.20) 



fe,i=l,2,3 



and Eq. |7.18| thus contributes the following terms as the massless neutrino coun- 
terterm Lagrangians: 



E 

fc,i=l,2,3 



{SZl(Kl) ie (K L ) ek + 5Z£(Kl) iti (K L )^ + 5Zl{K{) iT {K L ) Tk }v-^v k . (7.21) 



In our case, however, we sum over internal u k but not over external z/j. The graphic 
representation of the relevant counterterm (embedded in muon decay) is in Fig. \l.3[ 
The amplitude for this diagram is 



Figure 7.3: Counterterm diagram for neutrino self-energy in muon decay 
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Mc = -\M™ £ E «M)aWA- (7-22) 

l=e,fi,r fc=l,2,3 

Again, the factor | comes from our dealing with the external wave function rather 
than the internal propagator and the mixing factor {K\) k originates at the fiWfk 
vertex. The amplitude M.q can be further simplified, 

Mc = -~M™ £ 5Zi{Kl) a £ {K L ) lk {K\) hll 

" l=e,fi,r fc=l,2,3 

= -\m£Z £ SZ^KlU^-l^) = -\M™JZ£(l-w mix )(Kl) ifl 

= -]-5Z£(l - HHmix)Mtree- (7.23) 

The factor (K^)- was absorbed by M tree = Mf r f e {K.\).. 

Now we can write down the final expressions for the renormalized amplitude 
M. se if and the renormalized neutrino self-energy: 

E" ! 5Z l 

M self = M self + M C = ^Mtree- ^(l-Um^Mtree, (7.24) 

t v l = ^ + 5Z l L (l-ll mix ). (7.25) 
The constant 5Z\ is found from 



6Z l L = -s^K)-mn^K 2 ) + s^K) + 2S^K 2 )], 

dp 2 



KnA m i) = -^?K 2 ). (7-26) 



where E^, E# and E5 are respectively the left-handed, right-handed and the scalar 
part of the lepton self-energy given in Eqs. |6.41| - |6.44| []. The only graph with 



significant contribution to the term with derivatives, is the photonic loop (Eq. |6.44j ). 



3 Now we include the photonic loop , in the total lepton self-energy, see Eq. 6.44 and the 



J 7* 

following comments. 
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The other diagrams lead to derivatives of Bq IU , B( in functions which, when multiplied 
by mf, are of the order mj/Myy and therefore negligible. This can be easily seen from 
Eqs. E.18| , E.19 , wherein also derivatives of the Bq IH , B( m functions corresponding 



to the photonic loop are given. The final answer is 

5Z[ = -Ei(mf) + ^(21n^-l). (7.27) 

To prove the cancellation of the infinities, we note that the infinite part of 5Z l L 
is given by 



and the infinite part of the neutrino self-energy by 

47T S]y 4 2 4c% 

+ -i^(l-// m „) 2 }A M . (7.29) 
From the formulae above it can be easily seen that infinities cancel out in Eq. [7.25 



7.2.2 Limit M N > M w , M z , M H 

We will investigate the large behaviour of the renormalized neutrino self-energy 
E^. For large Mn we get from Eq. E.18 



B (p;M H , ZiW ,M N ) = l-21nM JV) 
Bi(p;M HtZtW ,M N ) = -0.25 + \nM N . 



(7.30) 
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This implies the quadratic nondecoupling for £|f(p) and see Eq. |7.16| : 

a 1 



Sf(p) + Sj(p) 



27r4^ 



m 2 n n 



IV 



-A,, + 0.75- In 



(7.31) 



and for the , a left-handed part of E^y (see Eq. |6.43|) , which contributes to the 
fl 1 ^ via the counterterm SZ l L , see Eq. |7.27| : 



a 



11 , 



Ml 

'M 2 W 



A„ + --2\nM N 



(7.32) 



From here we can see the not only cancels out infinities in the £^ (p) and (p) i 
but, in the limit investigated, it also cancels out the finite parts. As a result, there 
is no quadratic nondecoupling in the renormalized neutrino self-energy. 



7.3 Vertex diagrams 

Diagrams modifying the WfiUi vertex are depicted in Fig. |7.4j. Another set, one 



that modifies Wevj vertex, is not shown. 



The sum over the depicted set of diagrams gives the muon vertex amplitude 



vertex • 



M 

vertex M^z + M^ NZ + M 

+ M WZN + M <t>ZN + M-WHN + M-<t>HN + M<pxN- (7.33) 



The computation of the diagrams yields 
M vertex - M tree —< (A A / Z - -J(l 



2^ - 1 , 1 _J^_ ln ^lx 



fit/ j' 




Figure 7.4: Vertex diagrams for muon decay 
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(7.34) 
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where 

Ml Ml Ml M% m 2 

Vim) = 111 77 — — T~p\ TTT^-^TTT 5 T7T ^ n TTT- (7.35) 

m 2 Ml - Ml Ml m? - M 2 N Ml v ; 

The vertices are renormalized as (see Eq. |E.3| ) 

A" = A" + 5ZY - 5Zf + 5Z£, (7.36) 

where 



= MZ ertex /M tree , (7.37) 

SZ^-dZ™ = -J?L-A„ W . (7.38) 

Looking for the dominant graphs in the limit M^ 3> Mw, Mz, Mh, we note the 
leading contribution to the function v(m) is 

B(mr . = _L( m M„^-<ln^), (7.39) 



which implies the graphs of Fig. |7.4j f have quadratic nondecoupling However, as 



in the case of the neutrino self-energy, both infinite and finite terms of these graphs 

f 2 N 



are cancelled by the term in the counterterm 5Z l L . Therefore there are no M| 



dependent terms in the renormalized vertex diagrams either. 



4 Another interesting point the authors of Ref. 
that the diagrams for the muon decay (tt decay in Ref. 64 ), Figs 
out as dominant in an elegant way by using the Goldstone boson equivalence theorem. 



make (see also footnote 2 on page 122) is 
|7l] e, f^2]c, [M|f can be singled 
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7.4 Results 

The muon decay loops modify the Ar quantity in the implicit relation between My/ 
and G M , see Eq. |6.61| : 

2 2 7ra / I 1 \ 

M w s w = - - — (1 - -ee mix - -^mfa j , (7-40) 

and Ar can be written as (see Eq. |6.62| ) 



TZet w (0) . 
Ar = M 2 + 5 v- (7.41) 

Here the parameter 5y is the sum of the loop diagrams calculated in this section: 

Sy = ^W t+ ^ + ^e_lf,u e _ij : u, + Mbo^ (?42) 
Ai tree 2 2 Alt ree 



Numerical results are shown in Table |7.1|. As input data we used masses from 



the standard set, mixings ee mix = 0.0071, fifi m ix = 0.0014 and TT mix = 0. We 
suppressed the TT mix , which at its maximal value currently allowed (0.033) would 
make corrections to £pp(0) / 'M^ much larger than the corrections to 5y (these only 
depend on ee mix , fifi m i X ) and we would like to show the case when also the latter are 
important. In the first three lines of the table we show how much the self-energy, 
vertex and box diagrams contribute to 5y (line 4) for NHL masses of up to 30 
TeV. Also shown (lines 5,6) are Sp^(0)/M^ and Ar. Ultimately we are interested 
in NHL effects in the observable Mw (line 7). 

The results confirm expectations from the previous sections. There is no nonde- 
coupling for self-energies and vertices and there is a quadratic dependence on M/v 
for the boxes. The boxes are becoming important at very high masses. Still, they 
are small compared to the change in £jy(0)/M^. This is due to the fact the dom- 
inant boxes depend on the product ee m i X fifi m i X (see Eq. |7.13|) , while the correction 



132 





SM 


0.5 TeV 


5 TeV 


15 TeV 


30 TeV 






- 4.995 


- 4.972 


- 4.982 


- 4.988 


- 4.992 


xl0~ 2 




- 1.441 


- 1.442 


- 1.444 


- 1.444 


-1.445 


xlO" 2 


Mbox/Mtree 


4.273 


4.300 


4.315 


4.457 


4.950 


xlO" 3 


5v 


6.670 


6.539 


6.525 


6.652 


7.133 


xlO^ 3 


±w(0)/M^ 


2.396 


2.346 


2.301 


1.872 


0.329 


xl0~ 2 


Ar 


3.063 


3.000 


2.954 


2.537 


1.043 


xl0~ 2 


M w 


80.459 


80.537 


80.545 


80.612 


80.846 


xl 



Table 7.1: Contribution of the muon decay loops to 5y and A r 



to the W propagator is proportional to knu = ee 2 ^ + ^^ ix (see Eq. |6.25| ), which 
is allowed to be larger given the current bounds on the mixings. 

The W mass jumps from M^ M = 80.459 GeV to M w = 80.537 GeV at M N = 0.5 
TeV mainly as a result of the tree-level correction ^1 — \ee m i X — ^nHmix^ > see Eq. [7.40 



After that it rises very slowly until the dependent amplitudes become dominant 
above 5 TeV. 

In conclusion, the numerical analysis of Chapter 6 turns out to be basically valid 
even after the restriction ee mix = \i\i m ix = is relaxed. It can be improved by the 
inclusion of the tree-level correction (l — \ee m ix — \^mix)i while the largest loop 
corrections, the box diagrams of Fig. |7.1|e, are only marginally important. 



Chapter 8 



Conclusions 



Two theoretical schemes were discussed here as possible solutions of the problem 
of the small neutrino masses: the see-saw mechanism of Yanagida and Gell-Mann, 
Ramond and Slansky g, [| (Sec. |2.3.2| ); and the superstring- inspired low-energy 
model of neutrino masses suggested in Refs. [[17], |27j (Chapter 3). Both these 
solutions introduce NHL's into the theory as a necessary ingredient. Our main 
consideration has been the phenomenology of NHL's in the superstring-inspired low- 
energy model. The qualitative features of our analysis are applicable also in the 
context of see-saw models with enhanced mixings. 

The superstring-inspired low-energy model (Chapter 3) is a simple extension 
of the SM. It enriches only the neutral fermion spectrum of the SM, leaving the 
gauge and Higgs sectors intact. We found that among the new parameters of 
the model, seven are especially important: 'flavour- conserving' mixing parameters 
ee mix , Wmix, TT mix ; 'flavour- violating' mixing parameters e/i mix , er mix , [iT mix and the 
mass scale of NHL's (we assumed all three NHL's have mass Mjy). The bounds 
on the six mixing parameters (Sec. |3.4| ) are largely independent of the mass Mat- 
On the other hand bounds on always depend on the mixings. 

The mass Mn, if larger than Mz, can presently only be probed in radiative 



133 



134 



corrections (loops). A traditional approach was mostly limited to hypothetical lepton 
flavour- violating processes such as \i — > ej; fj,,r — > ee + e~; Z — > e ± /i =F etc |20], [4^, ^6], 



[43] , |44] , |45| , [46l , [47| , |48| . We reviewed constraints from these processes in Chapter 5. 

Besides flavour- violating processes NHL's could also induce (via radiative cor- 
rections) deviations from the SM in currently observed processes. We calculated 
radiative corrections due to NHL's to three such observables: leptonic widths of the 
Z boson Tn, lepton universality breaking parameter Ub r and the mass of the W boson 
M w (Chapters 6 and 7). We found that these observables form three complemen- 
tary quantities as far as sensitivity to NHL masses and mixings is concerned. Yu 
depends on three kinds of radiative corrections: the vertex corrections (Sec. 



the Z oblique corrections and the W oblique corrections (Sec. 6.3); the universal- 
ity breaking parameter depends only on the vertex corrections; the W mass Mw 
depends to a large degree only on the W oblique corrections. 

The effect of the NHL mass Mn in radiative corrections is, on the one hand, 
suppressed by the small mixings; on the other hand it is enhanced due to the violation 
of the Appelquist-Carazzone theorem (Sec. |6.6p . These competing tendencies are 
reflected by the typical behaviour of the dominant terms in radiative corrections due 
to NHL's (see Eqs. |[3(| gjjp , 

To make up for the small mixings, only NHL's with masses in the TeV range can 
lead to significant deviations from the SM. 

Assuming TT mix = 0.033, ee m i X = unmix = and m t = 176 GeV, we derived the 
following limit on the NHL mass at the 2a level from the Z decay width to r leptons: 

M N < 4.3 TeV. (8.2) 
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The universality breaking ratio, Ub r , leads to a yet better upper limit, 

M N < 3.8 TeV, (8.3) 

at the 2a level. We can use Eq. |871] to display the approximate dependence of the 
above limits on rr m i X : 

D33 

M N < 4.3 x TeV (8.4) 



from Eq. |S.2| and 



fl D33 

M N < 3.8 x TeV (8.5) 



from Eq. |3.3|. Note the limits of Eqs. |3.2|, 8.3 are comparable to the limit 



M N < 4 TeV, (8.6) 
derived from the considerations of perturbation theory breakdown in Sec. 1677 . 



We also found some sensitivity of the W mass My/ to the NHL mass and mixings, 
which depends to a large degree on the top quark mass (Figs. |6.17| a,b). 



In Chapter 7 we generalized our analysis of Chapter 6 by relaxing the restriction 
ee m j X = nHmix = 0. We found that while the numerical results of Chapter 6 remain 
basically valid, they can be improved by the inclusion of the tree-level correction to 
the muon decay, (l - \ee mix - \^ m i^j . 

As already noted in Chapter 6, we feel there are at least two reasons which 
give the (flavour-conserving) processes studied in this thesis a distinct advantage 
over the flavour- violating ones. First, the limits on which we derived are only 
matched by those from /i — > eee. The flavour-violating decay rates for r and Z are 
below the current experimental sensitivity (see Sec. |5.1.3| and |5.2p . Moreover, the 
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[i — > eee decay depends only on ee™,^, e/j, m i x , the two of the six mixing parameters 
(see Eq. |5.30| ), and may be unobservable if ee m ix or e\i m ix are very small. Second, 
the inequality Eq. [3.27 can further suppress the flavour- violating processes against 
the flavour-conserving ones via the 'conspiracy of the phases' in the sum of complex 
terms making up the flavour-violating parameters. 

For these two reasons, first signatures of neutral heavy leptons could come from 
flavour-conserving observables. At this time, LEP has stopped its runs at the Z-peak 
energy and is running at 130 — 140 GeV. It will eventually be producing W pairs 



which will allow the mass My/ to be measured with a precision of 0.05 GeV ||63 



(currently My/ = 80.410 ± 0.180 [j53"|). Combined with more precise measurements 
of the top quark mass we might be in a position to place even more stringent limits 
on NHL masses and mixings from our prediction of My/ (see Figs. |6.17j a,b). 

The observation of neutral heavy leptons is essential for our understanding of the 
small neutrino masses. It would provide us with significant hints on grand unified 
theories and possibly superstring theories. 



Appendix A 

Dirac algebra and trace theorems 



In this Appendix we define the Dirac gamma matrices and collect their properties. 
We also show spin sums, some Fierz identities and the proof of the identity v c L v R = 
vlv R . 

A.l Gamma matrices 

The so-called Dirac representation of the gamma matrices is given by 




where I is 2 x 2 identity matrix and the Pauli matrices o % are defined as 



1 




1 







\ 



(A.2) 



V 



1 



-1 



/ 



Some properties of the gamma matrices follow: 



7o = 7°, 



7i = -7, 
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7o 


= 7o, 


7o7o 


= 1, 


707^70 


7/i? 


T 

7o 


= 7o, 


7 5 T 


= 75, 


ll 


= 1, 


v\Yv 2 


= -v 2 T Y T v x 



where vi, u 2 are spinors, 1 is 4 x 4 unit matrix, and T represents product of gamma 
matrices. 

Trace theorems: 



Tr(l) 

Tr(Yl") 
Tr(Yl"l X l a ) 

Tr( l5 YY) 

Tr( l5 Yl u l X l a ) 
Tr(odd number of gamma matrices) 

Product rules: 



= 4, 

= 4<r, 

= 4(^V CT -^V <T + ^V A ), 

= o, 

= o, 

= 0. (A.4) 



7Y + 7Y = 2<T, 

7^7 M = 4, 

7 „yy = -2y, 

7 yyy = 4 «f\ 
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lifl 7 7 7 = ~ 2 1 7 7 , 
757^ + 7 M 75 = 

<r<^ = 4. 

An explicit representation of the charge conjugation matrix C is 



(A.5) 



2„,0 



C = ^7 7 







1 





1 

-1 





A. 2 Spin sums 



s 

s)vp(p, s) 

s 
A 



-1 







c 1 c 



-1. 



(^-™)a/J , 



(A.6) 



(A.7) 



where u a (p, s), f a (p> s) are spinors with momentum p, spin s and mass m; e M (p, A) 
is a polarization vector of a weak boson V = W, Z with momentum p, spin A and 
mass My. 

A. 3 Fierz identities used for the calculation of the 
boxes 



[wy a (l - 75)7e7 7 M ] x [^e7 Q (! - TsbV^J = 16[m7 m (1 - 75) it] 

x [U e y(l-75K e ], (Ai 
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[«7a(l - 1^)1 el ^U] X [W e 7 7 (l - 7 5 )7 £ 7a^J = 4[«7^(l - 75) It] 

x [tvf(l-75He]- (A.9) 

A. 4 Proof of the identity Vl v r = vlVr 

Using the definition of the charge conjugate field, 



we get 



v c = Clou*, (C = t 1 2 1 °) 



(A.10) 



L R 



1-75. 



1 +75 c ct l -75 f 1 +75 c 



= ^ ct 7o 



1 + 75 1 + 75 c = ^ 1 + 75 



T/^y 1 + 75 ^-y * *T 
= Z/ 6 — CloU = — V 



C 



1 + 75 



C 7 G 



-i/t 



7o^i±^^ 



1/ = -vC 



-nTnT 1 + 75 _ 1 +75 
— vC C v = v i 



(A.ll) 



Appendix B 

Couplings of v and TV to Higgs 



We begin with some useful properties of the rotation matrix G (see Eq. |3.12| ). From 
GGt = G^G = 1 we have 



UiU{ 


+ u 2 u\ 


= 1, 


uli/! + u\u 3 = 


1, 


u 3 ui 


+ u A u\ 


= 1, 


u\u 2 + u\u 4 = 


1, 




+ u 2 u\ 


= 0, 


u\u 2 + U\Ui = 


0. 


u 3 u\ 




= 0, 


V\U X + U\U 3 = 


0. 



(B.l) 



Further, from (see Sec. |3.1.2|) 



U X D + U 2 M = 0, 
U 3 D + U A M = M' 



(B.2) 



we get 



UlUxD + UlU 2 M = 0, 



U\U 3 D + U\U A M 



u\m'. 



(B.3) 
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Adding these two equations and using Eq. [B.l| we find the following relation between 
M and M'\ 



M = uIm'. 



(B.4) 



To derive the Lagrangian describing the couplings of v and iV to H, Ch, we 
rewrite Eq. [3.30 



as 



C 



92 



w 



VlIl) D$n R + h.c. , 



(B.5) 



where D is a 3 x 3 matrix in family space and, (see Eqs. [2.16| , [2.23| ) 



i,0* 



V 



(B.6) 



Selecting the H part we get 



C 



92 



2M, 



-vlDhrH + /i.e. 



(B.7) 



In the next step we add and subtract a term: 



H 



92 



2M, 



-vlDtirH + /i.e. 



IV 



,/2 S L Mn R H + h.c. 



+ 



2M W 
2Mi 



S L Mn R H + h.c. . 



w 



(B-8) 



The first two lines of this relation can be compared with Eq. [3.10| . We can now use 



the results of Sec. 3.1.2, which give 



H 



' J - S' L M' n R H + h.c. 



+ 



2M W 

92 f— TT , TT 



2M, 



w 



(y' L U 2 + S' L U 4 ) Mn R H + h.c. 
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■J?-SlTM'Z^n R H + h.c. 
2M W 



+ (v' L U 2 MZ ] + SlTU±MZ ] ) n' R H + h.c. . (B.9) 



Now we use M = U\m' (see Eq. ^§ and M" = TM'Z^ (see Eq. |07|) : 



£ H = --^—S" L TM 'Z^n' R H + h.c. 



2M W 
92 (— 



2M, 



n R H + /i. 



^ 2 YlM"n' R H + h.c. 



2M W 

+ (^"[^It+TMV + SlTU A U\T^TM' Z^) n" R H + h.c. 

2Mw v / 

^ 2 SlM"n" R H + h.c. 



2M W 

+ (V L U 2 UlT ] M" + YiTUJj\T ] M") n R H + /i.e. . (B.10) 
Using Z7 2 t4 = -C/iC/J and C/ 4 £/"| = 1 - E/gL/J, see Eq. Q 

+ r§-K(-CiC3 , )TtAr + Slr(l-t/ 3 C/t)rtAf"] n ' R ff + A.c, (B.ll) 
and putting 



we get 



Kr 



K 



H 



(B.12) 



H 



U2 ~SlM" n R H + h.c. 



+ 



+ 



2M W 

#2 



2M W 

#2 77T 

2M, 



v' L {-K\K H )M n R H + /i.e. 



-5^(1 - K ] H K H )M"n R H + /i.e. 
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-^Sl(Kl I K H )M''n R H + h.c. 
^^L(KlK H )M"n R H + h.c. 
-7^T W l{K ] h K h )M" N R H + h.c. 

92 



2M, 



-v' L {K{K H )M"N R H + h.c. 



w 



-^L-N{K^ H K H )M N NH 
g 2 — It a T r N ,, 1 +75 



2M W * 
The couplings of x, + an d 0~ are found by analogy. 



2M W y u 11J 2 
■ ll ~N(K\jKl)Mn — ~—^-vH. (B.13) 



Appendix C 
Feynman rules 



We list here the Feynman rules needed for the computation of the non-SM dia- 
grams contributing to the processes studied in this thesis. They are given in the 



t Hooft-Feynman gauge (see Sec. O). The rules for the vertices correspond to the 



interaction Lagrangians of Sec. |3.3|. The SM case is obtained in the limit 



K 



H 



0, Kj 



(C.l) 



In vertices, where applicable, the arrows indicate in addition to the flow of the charge 
also the flow of momenta. 
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Asw c w 



l 3 
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+ 3^(*!*hL7m(1-7b) 




N, 



N. 



+ 



(K* H K H ) abl »{l- l5 ) 




w 



/ 

/ 



+ ^(^l)h7,(1-7 5 ) 



+ ||(^)^(l-75) 



N, 



»g2 

2M W 



m n (k[k h ). 

\ / la 



(1+75) 
2 



N 



a 



W2 
2M W 



m n (k^ h k l ) ( 



(1-75) 

2 



N 



ff2 
2M W 



M N (K{K H )^ 



N 



N. 



a 



(1-75) 
2 



W2 
2M W 



M N (K\jK H ) 



ah 



N 



N 



ff2 

2M W 



M N (K* H K H ) ablb 



N, 



+ 7^ ( R H)la 



mmi t jy- \ (1-75) 
V2M W \^ L 'U 2 



N. 



+ 



W2 



V2M W \ HJ al 



M N £=^±- mi ^ 




mm ( js\ \ (1+75) 

V2M W \ LJu 2 



3 } 1 W. (p.) 
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w 
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1 



X 
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r \Jf\J~< 



4> 
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w 

4> -(p_) 

H ( Pl ) 
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z„ 
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Appendix D 

Dimensional regularization and 
some useful integrals 



D.l Dimensional regularization 

Before one-loop amplitudes with divergent momentum integrals can be renormalized, 
they have to be regularized. Regularization defines integrals, parametrizes their 



divergences, and separates their finite parts. Dimensional regularization [[40] defines 
integrals by analytically continuing them from 4-dimensional to n-dimensional space- 
time. The computation of integrals in n dimensions typically yields (see Eq. |4.9|) 

^137 ( 2 \ 

Ay(0) = ( — V finite constants ) , e = 4 — n, 

A-k \e J 

that is, the divergence is parametrized as a simple pole at n = 4. In our calculations, 
we used the following momentum integrals in n dimensions: 



d n k 1 i(-ix) n l 2 T{l -n/2) 



(2tt)» (A; 2 + 2k ■ s + t) 1 "~ (2n) n T(l) (t - S 2 )('-"/ 2 ) 
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V(0 





-s 2 ) 


(l-n/2) ' 




r d n k 




k/j, 




J (2tt)" 


(k 2 


+ 2k-s 


+ t) 


r d n k 




k^ky 




J (2tt)" 


(k 2 


+ 2k-s 


+ t) 


1 









X 



l-n/2 - 1 
Specifically, for I = 2, 3 we have 

iHQj r(2-j) 



iV(2) 
iV(3) 

where we used 



(2vr) n T(2) 
i(-7r)tr(3- 



(27r) n T(3) 



— ^-(1 + -ln47r ( 7 

(4tt) 2 V ^ 2 [ 

2(4tt) 2 V 2 A 2 [ 



r(2-=) = r(2-izf 



2 

- -7, 



n 

K3-- 



r i + - 

' 2 



£ rf- 

2 V2 



1--7, 

2 [ 



(D.l) 
(D.2) 



(D.3) 



(D.4) 
(D.5) 



(D.6) 



where 7 = 0.5772 is Euler-Mascheroni constant. In n = 4 — e dimensions, a becomes 
a dimensional quantity: 



a = > aa € = a( 1 H — In u 2 + ... 

47T V 2 



(D.7) 



where /x is an arbitrary mass scale. The combination of Eq. |D.4| and Eq. p.7| yields 



a 2 2 
a/i e (l + -ln47r)( 7) = a( 7 + In 4vr + In /i 2 ) = aA^. (D.8) 



It is this A M rather than -, which is usually thought of as the parameterization of 
the divergence since factors 7,ln47r and ln/i 2 are always present along with 2 and 
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they all together cancel out in renormalized quantities. We use also other variants 
of the A symbol: 

2 2 

A = 7 — In 7r = 7 — In 7r, 

4 — n e 

2 

A M = --7 + ln47r + ln/i , 

2 m 2 

A m = 7 + ln47r-ln— . (D.9) 

e n 2 



To cast momentum integrals into the form of Eq. p.l , the following Feynman 
parameterization is used: 

\ rl rxi rx n -i 

= T(n + 1) I dxi / dx 2 --- I dx n 



aQa\a2...a n Jo Jo Jo 

x T 7 x —7 \ TTT- ( D - 10 ) 

[a + (ai - a )xi + ...{a n - a„_i)x n J n+i 

Higher powers in a; are obtained by differentiation with respect to this parameter. 
In this work, these specific expressions were used: 

If 1 1 

-r = / 7a Ha' ( D - u ) 

ao Jo |a + (o — ajxp 

-1- = 2[ 1 dx[ X dy r -— (D.12) 

aoc Jo Jo [a + [b — a)x + (c — oji/J' 3 

= 2 [ dx [ dy ~, (D.13) 

abc Jo Jo [axy + by(l - x) + c(l - ?/)] 

1 = Q f 1 dx f dy - . (D.14) 

abc 2 Jo Jo [a + (b — a)x + (c — 6)y] 4 

In n dimensions, the algebra of the Dirac matrices described in Appendix A is 
generalized as follows: 

YY + YY = 2g" u , 
g^Qnv = n, 
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7m7 m 


= n, 






ill/ 

= rig 1 ", 






= n{g^g x °- 






= (2-n)-f, 






= Ag uX + (n- 


4)tV, 




= -2 7 a 7 V" 


-(n-4) 7 V7 CT , 


7s7 M + 7 M 75 


= 0. 





D.2 The computation of 2o,2i(ra), 22(771) and X^(m) 
integrals 

We start with two useful integrals one often encounters during the computation of 
momentum integrals: 

r 1 D C + D 

/ dxln(Cx + D) = ln(C + D) + — In — ^ 1, 

JO / / 1 

/ cMaTn(F:c + F)] = 
Jo 

+ 
+ 

To compute Xo,X 1 (m),X 2 (m) and X 3 (m) we note these integrals (defined in 
Sec. f^lf) are related via 

X 2 (M Z ) = Ji + M 2 J 3 (M Z ), (D.17) 
and the integral Xq is obtained as a special case of Xi(m) for M^r — >■ 0. To calculate 



ln(F + F) 



2F 2 

1 F 2 F(F + F) 
4~EP F 2 
F 2 InF 



2E 2 



InF 



1 (E + FY 
2E 2 2 



ln(F + F) + 



FF 
F 2 " 



F 2 



(D.16) 
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the Ti(m) integral, we use the parameterization of Eq. p. 13 



— — = 2 f dx f dy- 



abc Jo Jo [axy + by(l — x) + c(l - y)] c 

The term in the denominator is 



axy + by(l — x) + c(l — y) = {k 2 — M^)xy + {k 2 — m 2 )(y — xy) 

+ (k 2 - M 2 N )(l-y) 

= k 2 — M^xy + m 2 xy — m 2 y 

- M 2 N (l-y). (D.18) 



Using the momentum integral Eq. |D.1| , we can write for T\ (m) 

d 4 k 1 



(47T) 



-Urn) = J 



x 



(2tt) 4 (k 2 - M%){k 2 - M^k 2 - m 2 ) 
f 1 r f 1 , f d*k 
2 Jo dX Jo dy Jj2^ 

y 

[k 2 — Myyxy + m?xy — m 2 y — Mjy(l — y)Y 
1 rl 



I dx dy 
o Jo 



x 



(m 2 - Mfr)yx - m 2 y - M 2 N (l - y) ' 



(D.19) 



First we integrate over the x parameter: 



rl rl _1 

T x (m) = J o dyyJ Q dx ^ _ - — - - ^ 

dyy I dx ~ = - / rfyyi[ln(A + 5) -InB] 

'^H , 2 l M2 , [HM^y + M 2 N - M 2 N y) 
io (m 2 - Mfc)y 

- \n(m 2 y + M 2 N - M 2 N y)], (D.20) 
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and then the integration over y: 



Mm) = Cdy 2 1 [\n(M^y + M 2 N -M 2 N y)-\n(m 2 y + M 2 N -M 2 N y)] 
jo (in iyiy^) 



(m 2 - M 2 / 



■ [ £ dy ln(Cy + D) - £ dy \n(Ey + F) 



F E + F 

- t; ln Z, + 1 
E F 



1 



(m 2 - M 2 ,) 



ln4 + 



M 2 



ln 



w 



M 2 



ln 



M^-M 2 , M 2 r m 2 -M 2 r M 2 



N 



N ~ 



.(D.21) 



The X 3 (m) integral is found using similar steps: 



X 3 (m) = 



(4tt) ; 



= 6 



/ J (27T) 4 (k 2 -M 2 N y(k 2 -M 2 v )(k 2 -m 2 ) 
d A k 



(2k) 



rl rx 

/ dx dy 
Jo Jo 



y 



dx / dy 
o Jo 



[k 2 - M 2 , + (M 2 , - m 2 )x + (m 2 - M 2 N )y\ 

y 



[(M 2 , - m 2 ):r - M 2 , + (m 2 - M 2 ^] 2 



+ 



/If 2 In 



m 2 - M^XMl - M 2 , (Ml-M 2 ,) 2 M 2 N -m 2 

^ 2 j- 



m 2 ln^ 



(M 2 , - m 2 )' 



(D.22) 



Appendix E 

Renormalization constants, 
unrenormalized self-energies and 't 
Hooft scalar integrals 



E.l Renormalization constants 

It is convenient to define the following linear combinations of the renormalization 
constants (i = 1,2): 

5Z] = s 2 w 5ZY + c 2 w SZf, 
SZf = c* w 5Z? + s 2 w 5Zf, 

szf = (szf _ 5Z1 ) , 

c w s w 

SZl = (5Z f L + 5Z f R )/2, 6Z f A = (5Z f L -5Z f R )/2. (E.l) 
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The renormalized self-energies are obtained from unrenormalized ones by adding 
appropriate counterterms: 

s 7 (p 2 ) = £ 7 {p 2 ) + sz] P \ 

t z (p 2 ) = [f) - m\ + szi (p 2 - Ml 

t w (p 2 ) = Z w (p 2 ) - 5M 2 W + 5Z? (p 2 - M 2 

t~< z (p 2 ) = (p 2 ) - 5Zfp 2 + (bZf - 8Zf) M 2 , 

±f(p) = j6(4(p 2 )+54) + fi l5 (z f A (p 2 )-5Z f A ) 



1 > 



+ m f (zi(p 2 )-5Zl- 5 -^), (E.2) 

where Hy,H A and Eg are vector, axial vector and scalar part of the fermion self- 
energy, respectively (see Eq. |6.41|) . 

The renormalized electromagnetic, weak neutral and charged current vertices are 
given by 

f iff = T yff + %txzilii ( SZ i _ §z i + 5z f y _ 5z f Al ^ 

- iew^Vf - a/75) (6ZJ Z - 5Zj z ^j , 

fzff = r z " + ie lsrfli {v f -a n5 )(6Zf-6Zi)-ie 1ST7ll (SZ? z -6Z^ z ) 

+iei37 r y»\Vf5Zl + a f 5Z A ) - ie 137 -f^ 5 (v f 5Z f A + a f 6Zf) , 
\wiu _ T wiu + i _^r_ J ^ 1 _ / x + 5z w _ 6z w + SZ A (E 3) 



The renormalization constants are obtained from the OS renormalization conditions 
Eqs. [4.48| - [4.49l (we show only constants needed for our calculations): 



bM 2 w = ReZ w (M^), 



SM 2 Z = J ReS z (M|), 
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ap 2 ch/ M| 

x 7 z _ o ^-4^ 7Z (0) ^-^ / jjg _ 5_M$A 

2 " J SwCw Mf + ^ ^Mf M*J 

1 ~ Sp» W SwCw Mf + s 2 w \m 2 Myy 

2 dv 2 s w M\ 

5Z{ = +Z f A (m 2 ). (E.4) 




E.2 Unrenormalized self-energies in the SM 



Below we present complete SM gauge boson self-energies corresponding to Figs. 5.2 



- |6.5| . They were calculated in Ref. f34| . For the definition of the function F see 
Sec. [E.3| ; for the definition of the A factors see Eq. p.9| ; s = p 2 , where p is the 
4- momentum of the gauge boson; w = Myy, z = M|, h = M\. 



4?r <-3 



3s A w - (3s + Aw)F(p; M w , M w )}, 



(E.5) 



^{ ~ 3 E Qf v f [ sA f + ( s + 2m 2 )F(p; m f , m f ) - - 



+ 



+ 



1 



Cy/Sw 
1 



2w 



A 



(34 + 

(34 + i) s + (Ac 2 w + |) w] F(p; M w , M w ) 



9cwSw J 



, (E.6) 
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+ 



+ 
+ 



£{« E 2 a ? s( A, + §-.„( ' - fc )) 



E + 4) ( sA / + ( s + 2mJ) F ^ m /' m /) " 3) 
3 



2 c 2 m f 



■(^f + F(p;m f ,m f )) 



w°w 



( 3 -^- + ^2-) s +( 4 + ^-^) M : 

(-c^(40s + 80u>) + (c^ - s 2 w ) 2 (Sw + s) + I2w) F(p; M w , M w ) 



+ f Itiz - 2h + s + ^ ^-)F(p;M H ,M z ) -2h\n- -2z\nzw 



/i 



+ (10z-2h + s)(l- 



s 

h + z 



In 



M, 



H 



h-z M ? 



In 



M H M Z 



w 



) 



w 



4) 2 - 4c 2 w ) 



10 r 2 2 



(E.7) 



+ 



+ 



W l=e,n,T 



+ _^) F(!); o, mi ) + ? s - f 



2s 

A. 3L 2 V 



+ ( 
+ ( 



m 
2 

—) + —(' 

q— doublets 

m 2 + + m 2 _ {m 2 + — m 2 _) 2 



3 

5 2 m 2 _ \ A_ / 5 , 



s - 
19 



m 2 , + mi 



2 - + 3 »( 1 -f) 



1 

2 

\± 

c w' 



2 ' 2 

)F(p;m+,m-) 
m+ + m 2 , \ s 



2 2 
mi — mi 



lnm + m_j 



A 



s 4 _ ^Lh z + 7w + 10a - 2 ( ^ 



1 / s — u>) 2 



— - u> + z 



F(p; M z , M w ) 



<? 2 



2s 

- ~-Z-(-4w-10s + —)F(p;0,M w ) 
3 v s ' 

1/ , s (/i-u>) 2 



2s 



-)F(p;M H ,M w ) 
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' w (jz + 7w + 10s - A(z - w)) - s^z + U2w 



3z z 
m — 

z — w w 



3 v v '> " 6 V - 2 

v 3 12' h — w w 3 v 3 7 

+ g^ + ^-^-^-X^ + f 8 )}- (E.8) 

E.3 't Hooft scalar integrals 

Here we define various C,B,A and F functions and reduce them to scalar integrals 
Cq, Bq and A$. For the calculation of Co and Bq we refer the reader to the original 
work of 't Hooft and Veltman, Ref. ||51|| . 

The Cq function is defined as (with finite parts indicated by the superscript): 



Co(m 1 ,m 2 ,m 3 ) = Co(pi,p 2 ;m 1 ,m 2 ,m 3 ) = C^ n (mi,m 2 ,m 3 ) 

d n q 1 



(E.9) 



where 

D = (q 2 - m\ + ie) [(q - pi) 2 - m\ + ie] [(q - p\ - p 2 ) 2 - m\ + it}. (E.10) 
The functions CV,- are defined by: 

C/j, — — [ — 7-77 = —PifiCn ~ P2^Cu, 

J 27T Z D 



C^ u 



d"q q^q u 
in 2 D 

PlnPluC 2 l + P2fiP2uC 22 + {p\^p 2v + PluP2 f i)C 2 3 - g^C 2 4- (E.H) 



The functions Cn, C 2 4, C 23 are reduced (in the limit p\ = p\ = mj <C (pi + p 2 ) 
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Mf , applicable for our considerations of the leptonic decays of the Z boson) to: 

t- 1 

B o n (pi + P2; m 1 ,m 3 ) + B^ n (pi, mi, m 2 )], 



C 11 (m 1 ,m 2 ,m 3 ) = Ci^fa^m^m^ = -j^[f 2 C Q (mi,m 2 ,m 3 ) 

z 



1 , in 

C 24 (wi, m 2 ,m 3 ) = -A + C 2 7(mi,m2,m 3 ), 
C < 24 r l (m u m 2 ,m 3 ) = [ra?C (rai, m 2 , m 3 ) + /iCn(mi, m 2 , m 3 ) 
+ Bf n (p 1 +p 2 ;m 1 ,m 3 )](-^) + ^ 

C 23 (mi, m 2 , m 3 ) = ^"(mi, m 2 , m 3 ) = --L.[Bf OT (pi + p 2 ; mi, m 3 ) 

z 

+ B^ n (p 2 ;m 2 ,m 3 ) + / 1 Cu(m ll m 2l m 3 )] 

2 

+ C|4 n (mi,m 2 ,m 3 )^2-, (E.12) 

z 



where 



/ 2 = Mf + m 2 -™ 2 , 

/1 = m?-ml. (E.13) 



The functions _B , -Bi are defined as: 

d n q 



B (p;m 1 ,m 2 ) = 

J m A [q 2 — m{ + te ) K<? — P) ~ m 2 + te \ 

= A + BQ tn (p; m 1 , m 2 ), 

Q in (p; mi,m 2 ) = — dx In [p 2 x 2 + — (p 2 + — m 2 )x], 

•/n 

2 , ^TT7 ^ 2— T = 

in 2 (q 2 — m{ + «e) [(g — p) 2 — m 2 + zej 
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1 f . 

B 1 (p;m 1 ,m 2 ) = --A + B( m (p; m u m 2 ), 



B( in (p; mi, m 2 ) = / dx In [p 2 x 2 + m 2 — (p 2 + m 2 — m%)x]x, 
Jo 



^ v J in 2 (q 2 — m\ + ie)[(q — p) 2 — m\ + ie] 

= V^VvBn - g^B 22 - (E.14) 



The functions -B 2 i and _B 22 can be reduced to 
1 

Ip 
1 



B 1 = ^{-A (m 1 ) + A (m 2 )-(p 2 + m 2 1 -m 2 2 )B }, 



B 21 = ^{-A {m 2 )-2(p 2 + ml-m 2 2 )B 1 -mlB -l/2(m 2 1 + m 2 2 -p 2 /3)}, 
1 



B 22 = ^{ + A (m 2 ) - (p 2 + mj- m 2 2 )B 1 - 2m 2 5 - {m\ + m\ - p 2 /3)}. (E.15) 
The functions A are defined as 

d n q 1 [d n q I 



A (m) [ [ 

J in 2 (q — p) 2 — m 2 J 



in 2 (q — p) 2 — m 2 J in 2 q 2 — m 2 

= -m 2 (A - lnm 2 + 1), 



A,(p,m) = f^l- % = -p,A (m). (E.16) 

J iiT 2 (q — p) 2 — m 2 

Relations between F and B functions: 

(p;mi,m 2 ) = — H 5 ~ m V Inm-i + lnm 2 + B {p; m 1 , m 2 ), 

m\ — m 2 m 2 

F(p;0,m) = -1 + lnm 2 + B o (p;0,m), 

u( \ m 2 2 -m 2 l F(p]m 1 ,m 2 ) 1 

Si(p;mi,m 2 ) = 5 -B (p;m 1 ,m 2 ). (E.17) 

2 2 
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For s = p 2 small with respect to mf, m\, m 2 , we have 



F(p;m 1 ,m 2 ) = — -3 

(mf — m^J 



m^ + wig 



22 2 
mimi , mf 

— — In — - 

2 2 2 
mf — m 2 mj 



B (p;m 1 ,m 2 ) 



m? + m , mi , ^ . , 

1 5 5- In in mi — lnm2 + Cms), 

mf — wij wi2 



B (p;0,m) = 1 - 2mm + 0(s), 
B 1 (p;m 1 ,m 2 ) 



2 t«2 — mf 



mf + 7B2 m\m\ ^ mf 
2 m\ — m\ m\ 



- ^B (p;m u m 2 ), 



Si(p;0,m) 



-- + lnm + 0(s), 



Si (p; m,0) = -- + lnm + 0(s). 



(E.18) 



Finally, in photon loops we encounter functions with regularized photon mass 
m A -> 0: 



S (p; m A ,m,) | p2=m 2 
Si(p; m A ,m ; ) | p2=m 2 
#i(p; m,,m A ) | p2=m 2 



dB 
dp 2 

ggi 



(p;m x ,mi) 
(p;m x ,mi) 



2 - 2 lnm/, 

-- + lnm /5 

3 , 

-- + lnm/, 

-l(-l-ln^), 

mf v m/ 7 



p 2 =m? 



2m 2 ' 



(E.19) 
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